Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,966,751 C→T A112V (GCG→GTG)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,7510CT98.6% 273.6 / NA 74A112V (GCG→GTG) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (33/39);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCG  >  NC_000913/3966606‑3966895
                                                                                                                                                 |                                                                                                                                                
ttCCTCCGTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTgctg                                                                                                                                              <  2:764703/150‑1 (MQ=255)
ttCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGtgctgc                                                                                                                                             >  1:92732/1‑151 (MQ=255)
tccctcctTCCCGCAAACACCCCTTCCCTCGCCGTTCCTAATGACATTTCCTTTCCCCCAAGCAAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTgct                                                                                                                                               <  1:369643/147‑1 (MQ=255)
  cctccGTTCCGCAGCAACCTCCTACCGCGCCGGCCCGGATAACCTCTATGTTCCCCCTAGCAAAATCCGCCGTTTCACCCTCCGCGCTGGTGATACCATCTCTGGAAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTgctgct                                                                                                                                            <  2:964545/150‑1 (MQ=255)
  cctccGTTCCGCAGACATCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgct                                                                                                                                            >  1:961547/1‑150 (MQ=255)
      cGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGaa                                                                                                                                         >  1:551984/1‑149 (MQ=255)
       gTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGt                                                                                                                                      >  1:738700/1‑151 (MQ=255)
       gTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGt                                                                                                                                      >  2:146421/1‑151 (MQ=255)
       gTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGt                                                                                                                                      >  2:543694/1‑151 (MQ=255)
        ttCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAg                                                                                                                                       <  2:808776/149‑1 (MQ=255)
          ccGCAGACAGCTCCTCCCTCGCGGGTCCGGATGCCATCTAGTTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTa                                                                                                                                    <  2:714519/150‑1 (MQ=255)
                aCAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGATGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGt                                                                                                                             >  1:508014/1‑151 (MQ=255)
                  aGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGtt                                                                                                                            >  1:584195/1‑150 (MQ=255)
                   gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTa                                                                                                                           >  1:371519/1‑150 (MQ=255)
                    cccctccctcCCCGTTCCTGATGACATCTACGTTTCCCCCACCCAAACCCGCCTTTTCAACCTCCGCACGGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGt                                                                                                                             <  1:566700/145‑1 (MQ=255)
                      ccTCCCCCGCCGGTCCTGATGACATCTCCTTTTCCCCTACCCAATTCCGCCGTTTCACCCTCCCCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGa                                                                                                                                <  2:904937/142‑1 (MQ=255)
                      ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACtt                                                                                                                       >  2:415481/1‑151 (MQ=255)
                           cTCGCCGGTCCTGATGACATCTAAGTTTCCCCTAGCAAAATCCGCCGTTTCAACCTCCGAACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGAGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                                                                                                   >  1:464151/1‑150 (MQ=255)
                             cGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACaa                                                                                                                 >  1:5493/1‑150 (MQ=255)
                               ccGGTCCTGATGACATCTCCGTTCCCCCTAGCCAAACCCGCCGTTTCACCCTCCGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACCTCGACaaa                                                                                                                <  2:431052/149‑1 (MQ=255)
                               ccGGTCCTGATGACATCTACTTTCCCCCTAGCCAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                                                                                                              <  2:383742/151‑1 (MQ=255)
                                cgttccgaagcccATCTCCGTTTCCCCGAGCAAAATCCGCCGTTTCACCCACCGCACTGGAGATACCAGCTCTGTAAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTc                                                                                                                      <  2:332516/130‑1 (MQ=255)
                                 ggTCCTGATGACATCTACGTTTCCCCTAGCCAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                                                                                            <  2:189443/151‑1 (MQ=255)
                                 gctCCTGATGACATCTCCTTTCCCCCTAGCCAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                                                                                             <  2:71333/148‑1 (MQ=255)
                                  gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                                                                                            >  1:207875/1‑150 (MQ=255)
                                    cttgaggaCATCTACGTTCCCCCAACCCAAACCCGCCTTTTCACCCCCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaa                                                                                                          <  1:543308/148‑1 (MQ=255)
                                         tgaCATCTCCGTTCCCCCTAGCAAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGc                                                                                                     <  2:885926/150‑1 (MQ=255)
                                         tgaCATCTACGTTCCCCCTAGCCAATTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                                                                                                        <  2:614307/147‑1 (MQ=255)
                                            cttttcctttcccccTACCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                                                                                                       <  1:795719/134‑1 (MQ=255)
                                             aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                                                                                                 >  1:355101/1‑150 (MQ=255)
                                             aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                                                                                                 >  1:239802/1‑150 (MQ=255)
                                                 actttcccccTGGCCAAATCCGCCGTTTCACCCCCCGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaac                                                                                              <  2:71508/143‑1 (MQ=255)
                                                 aCGTTCCCCCTACCCAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaa                                                                                               <  2:797288/148‑1 (MQ=255)
                                                      tCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc                                                                                       >  1:433395/1‑151 (MQ=255)
                                                        cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc                                                                                     >  2:27389/1‑151 (MQ=255)
                                                         ccTACCAAAACCCGCCGTTTCACCCCCCGCACTGTTGATACCATCTCTGGTAAAATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct                                                                                    <  2:181319/151‑1 (MQ=255)
                                                             gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg                                                                                 >  2:517073/1‑150 (MQ=255)
                                                             gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg                                                                                 >  2:441172/1‑150 (MQ=255)
                                                               cAAATCCGCCGTTTCAACCTCCGCACTGGTTATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGGGCTGCTGAAAGTTGACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg                                                                                 >  1:279566/1‑148 (MQ=255)
                                                                    ccgccgTTTCACCCTCCGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt                                                                                  <  1:931294/142‑1 (MQ=255)
                                                                       ccgTTTCACCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAAc                                                                       <  1:101300/150‑1 (MQ=255)
                                                                          tttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                                                   >  2:200509/1‑151 (MQ=255)
                                                                             caccccccGCCCTGGTAATACCATCTTGGGTAAGTTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCGCTTTGAGAACTTAACCCCGCt                                                                 <  2:507965/144‑1 (MQ=255)
                                                                              aaCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGc                                                               >  1:730970/1‑151 (MQ=255)
                                                                                    cGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCaa                                                         <  2:907777/151‑1 (MQ=255)
                                                                                    cGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCa                                                          >  2:566151/1‑150 (MQ=255)
                                                                                          ggTGATCCCATCTCTGTTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACtctc                                                   <  2:665464/151‑1 (MQ=255)
                                                                                          ggTGATACCATCTCTGTTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCaaa                                                        <  2:905168/146‑1 (MQ=255)
                                                                                          ggTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACtctc                                                   <  2:293343/151‑1 (MQ=255)
                                                                                            tGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGt                                                 >  1:236883/1‑151 (MQ=255)
                                                                                              aTACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGt                                                 <  2:614715/149‑1 (MQ=255)
                                                                                              aTACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCt                                               >  1:879386/1‑151 (MQ=255)
                                                                                                 ccATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCt                                               <  2:317999/148‑1 (MQ=255)
                                                                                                 ccATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTg                                              <  1:954219/149‑1 (MQ=255)
                                                                                                 ccATCTCGGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTc                                                <  1:449936/147‑1 (MQ=255)
                                                                                                         ggTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAAc                                    >  1:888316/1‑151 (MQ=255)
                                                                                                         ggTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAAc                                    <  1:253585/151‑1 (MQ=255)
                                                                                                            aaGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTg                                 >  2:103290/1‑151 (MQ=255)
                                                                                                              gATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGt                               >  1:61093/1‑151 (MQ=255)
                                                                                                                ttCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTa                              <  2:144737/150‑1 (MQ=255)
                                                                                                                    cccgccGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACgg                          <  2:617416/150‑1 (MQ=255)
                                                                                                                     ccgccgAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGtt                        >  1:648570/1‑151 (MQ=255)
                                                                                                                        ccgAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCt                      >  2:611325/1‑150 (MQ=255)
                                                                                                                        ccgAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCACAACAAAATCCTCTTTGAGAACTTAACCCCGCTTCACGCAAACTCTCGTCTGCGTATTAAACGTAGTAACTTTTCt                      <  1:304508/150‑1 (MQ=255)
                                                                                                                          gAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACt                   <  2:885417/151‑1 (MQ=255)
                                                                                                                          gAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACt                   >  2:693191/1‑151 (MQ=255)
                                                                                                                              gaagGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGa                 >  1:267769/1‑149 (MQ=255)
                                                                                                                               aagGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGa              <  2:330741/151‑1 (MQ=255)
                                                                                                                                agGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAt             <  1:543603/151‑1 (MQ=255)
                                                                                                                                  gTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAttt           <  2:881160/151‑1 (MQ=255)
                                                                                                                                   tGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAttt           >  2:892490/1‑150 (MQ=255)
                                                                                                                                      aCGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACt       >  2:800077/1‑151 (MQ=255)
                                                                                                                                         cTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCt    >  2:561873/1‑151 (MQ=255)
                                                                                                                                           aTTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTcg  >  2:590418/1‑151 (MQ=255)
                                                                                                                                           aTTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTc   <  2:946591/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCG  >  NC_000913/3966606‑3966895

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: