Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,101,541 G→C 7.5% P28R (CCG→CGG)  yggN ← DUF2884 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,5410GC 7.5% 549.8 / 28.7 160P28R (CCG→CGG) yggNDUF2884 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (75/72);  new base C (3/9);  total (78/82)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 171 total aligned reads displayed.

CTGATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTT  >  NC_000913/3101393‑3101689
                                                                                                                                                    |                                                                                                                                                    
     aaTCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTg                                                                                                                                                 <  3:69476/149‑1 (MQ=255)
      aTCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGa                                                                                                                                                >  1:705756/1‑149 (MQ=255)
        ccTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACg                                                                                                                                              >  1:599124/1‑149 (MQ=255)
         cTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGc                                                                                                                                             >  2:536268/1‑149 (MQ=255)
         cTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGc                                                                                                                                             >  3:386770/1‑149 (MQ=255)
          ttCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACgat                                                                                                                                            <  3:32738/149‑3 (MQ=255)
          ttCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACgat                                                                                                                                            <  1:269101/149‑3 (MQ=255)
             gCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTTCCGTCTGGCGTGATCACCAGATTGCTGTTTTCGCCCTCCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACGATgtg                                                                                                                                         <  4:1083893/149‑3 (MQ=255)
              ccTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCAc                                                                                                                                        >  1:1358143/1‑149 (MQ=255)
                 gCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTgg                                                                                                                                     <  2:700440/149‑1 (MQ=255)
                   tCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTa                                                                                                                                   <  1:1477721/149‑1 (MQ=255)
                   tCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTa                                                                                                                                   <  1:349324/149‑1 (MQ=255)
                       gcTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCg                                                                                                                               <  3:70094/149‑1 (MQ=255)
                       gcTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCg                                                                                                                               >  4:340695/1‑149 (MQ=255)
                         tGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGc                                                                                                                             >  2:90294/1‑149 (MQ=255)
                          gggCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCg                                                                                                                            >  2:340090/1‑149 (MQ=255)
                             cggcATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGa                                                                                                                         <  1:406563/149‑1 (MQ=255)
                                cATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCg                                                                                                                      <  2:1512497/149‑1 (MQ=255)
                                 aTTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCgg                                                                                                                     >  2:294372/1‑149 (MQ=255)
                                 aTTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGAGAGCgg                                                                                                                     >  4:390572/1‑149 (MQ=255)
                                     aGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCAt                                                                                                                 <  3:973778/149‑1 (MQ=255)
                                      gggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCTTCACGCGGCGTGACGCTGCACTGGTGGTCGGCGTGAGCGGTCAtt                                                                                                                >  2:878600/1‑149 (MQ=255)
                                      gggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTg                                                                                                                                      <  3:1405535/127‑1 (MQ=255)
                                      gggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTg                                                                                                                                      >  4:1405535/1‑127 (MQ=255)
                                      gggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTa                                                                                                                                   <  1:453667/130‑1 (MQ=255)
                                      gggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCAtt                                                                                                                <  1:536268/149‑1 (MQ=255)
                                      gggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCAtt                                                                                                                >  4:1392071/1‑149 (MQ=255)
                                           tatTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGt                                                                                                                           >  2:240557/1‑133 (MQ=255)
                                           tatTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGt                                                                                                                           <  1:240557/133‑1 (MQ=255)
                                             ttGCTTACCGTTATACATCGCGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCa                                                                                                         <  1:257990/149‑1 (MQ=255)
                                                 ttACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGa                                                                                                     >  4:1256734/1‑149 (MQ=255)
                                                   aCCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGaaa                                                                                                   >  3:1354728/1‑149 (MQ=255)
                                                     cGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGt                                                                                                 >  1:511939/1‑149 (MQ=255)
                                                      gTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtg                                                                                                >  1:131763/1‑149 (MQ=255)
                                                       ttATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtgc                                                                                               <  3:201/149‑1 (MQ=255)
                                                       ttATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtgc                                                                                               <  4:1179192/149‑1 (MQ=255)
                                                          taCATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtgctgc                                                                                            <  2:1472385/149‑1 (MQ=255)
                                                          taCATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtgctgc                                                                                            >  4:414492/1‑149 (MQ=255)
                                                           acacCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTgcc                                                                                           <  2:1239433/149‑1 (MQ=255)
                                                             aTCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtg                                                                                                <  2:165115/142‑1 (MQ=255)
                                                             aTCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtg                                                                                                >  1:165115/1‑142 (MQ=255)
                                                             aTCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTgccgc                                                                                         >  3:1558514/1‑149 (MQ=255)
                                                             aTCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTgccgc                                                                                         >  4:1050302/1‑149 (MQ=255)
                                                              tCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCAt                                                                                                                 <  2:513298/124‑1 (MQ=255)
                                                              tCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTgccgcc                                                                                        >  4:824724/1‑149 (MQ=255)
                                                              tCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTgccgcc                                                                                        >  4:590518/1‑149 (MQ=255)
                                                                aCGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCcag                                                                                      <  3:1038114/149‑1 (MQ=255)
                                                                aCGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGCCATTGCCGTCACTGAAAGTGCTGCCGCcag                                                                                      >  1:39197/1‑149 (MQ=255)
                                                                        gTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCAttt                                                                              >  1:516238/1‑149 (MQ=255)
                                                                        gTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCAttt                                                                              >  4:1438322/1‑149 (MQ=255)
                                                                           tGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTgc                                                                           <  3:225103/149‑1 (MQ=255)
                                                                            ggCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTgcg                                                                          >  3:1233831/1‑149 (MQ=255)
                                                                             gCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCcagca                                                                                    <  1:883550/138‑1 (MQ=255)
                                                                               gTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGcat                                                                       <  3:1160767/149‑1 (MQ=255)
                                                                                   tCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATa                                                                   >  4:1328024/1‑149 (MQ=255)
                                                                                   tCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATa                                                                   <  4:1233155/149‑1 (MQ=255)
                                                                                     aCCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCTGGCGGCAGCActttca                                                                           <  4:1436498/139‑3 (MQ=255)
                                                                                     aCCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCTGGCGGCAGCActttca                                                                           >  3:1436498/1‑137 (MQ=255)
                                                                                       cAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTgg                                                                                                                                     >  3:1183764/1‑79 (MQ=255)
                                                                                          aTTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTcc                                                            <  3:76234/149‑1 (MQ=255)
                                                                                              ccGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTTTTCCCTCa                                                        >  2:1305346/1‑149 (MQ=255)
                                                                                              ccGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCa                                                        >  3:621015/1‑149 (MQ=255)
                                                                                              ccGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCa                                                        >  1:133082/1‑149 (MQ=255)
                                                                                               cGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGcc                                                                                                             >  2:640296/1‑95 (MQ=255)
                                                                                               cGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGcc                                                                                                             <  1:640296/95‑1 (MQ=255)
                                                                                                gTTTTCGCCCTTCACCTGCACGGTTTGCGGGATGACAATCACATCGTCACGCGGCGTGACGCTGCACCGGTAGTCGGCGTGAGCGGTCATTGCCGTCATTGAAATTCCTGCAGCAAGCAGCATTTTATGCATCATTTTCTTCCCTCAAg                                                      >  2:503324/1‑149 (MQ=255)
                                                                                                    tCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaa                                                  <  4:518446/149‑1 (MQ=255)
                                                                                                     cGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaa                                                 <  2:339576/149‑1 (MQ=255)
                                                                                                      gCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaa                                                >  4:1192579/1‑149 (MQ=255)
                                                                                                        ccTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATc                                              >  1:1201252/1‑149 (MQ=255)
                                                                                                        ccTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATc                                              >  3:568315/1‑149 (MQ=255)
                                                                                                         cTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCa                                             <  1:141020/149‑1 (MQ=255)
                                                                                                            cACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATATTCTTCCCTCAAGAAAAAATCAGGc                                          >  4:1032/1‑149 (MQ=255)
                                                                                                            cACCTGCACGGTTTGCGGGCTGACAATCACAGCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGc                                          <  3:390572/149‑1 (MQ=255)
                                                                                                                 gCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGc                                                                                                                                             >  2:521408/1‑45 (MQ=255)
                                                                                                                 gCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGc                                                                                                                                             <  1:521408/45‑1 (MQ=255)
                                                                                                                  cACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAt                                    >  1:836945/1‑149 (MQ=255)
                                                                                                                        ttGCGGGCTGACAGTCACATCGTCACGCTGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGTTAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAATATCAGGCTAGCATTTGCTg                              <  2:247827/149‑1 (MQ=255)
                                                                                                                         tGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTgg                             <  2:369044/149‑1 (MQ=255)
                                                                                                                           cGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcc                           <  4:1354728/149‑1 (MQ=255)
                                                                                                                              gCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCa                        >  2:446169/1‑149 (MQ=255)
                                                                                                                                tGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                      <  4:1493026/149‑1 (MQ=255)
                                                                                                                                 gACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACAGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAt                     >  4:1088783/1‑149 (MQ=255)
                                                                                                                                    aaTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                  >  3:1504961/1‑149 (MQ=255)
                                                                                                                                    aaTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                  >  1:877025/1‑149 (MQ=255)
                                                                                                                                       cacaTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGTTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATaa               <  2:1003234/149‑1 (MQ=255)
                                                                                                                                       cacaTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATaa               <  2:705756/149‑1 (MQ=255)
                                                                                                                                       cacaTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATaa               >  4:660209/1‑149 (MQ=255)
                                                                                                                                       cacaCCGTCACGCGGCGTGACGCTGCGCTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATaa               <  1:1224757/149‑1 (MQ=255)
                                                                                                                                          agcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCTCCAGATTGATAACaa            >  2:1177316/3‑149 (MQ=255)
                                                                                                                                          agcGTCACGCCACGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTccc                                                           >  3:1014320/3‑102 (MQ=255)
                                                                                                                                          aTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGcc                                                                                                             >  2:46340/1‑52 (MQ=255)
                                                                                                                                           tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa           <  3:1080635/149‑1 (MQ=255)
                                                                                                                                           gcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa           >  1:1162302/2‑149 (MQ=255)
                                                                                                                                             gTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGt         >  2:1521659/1‑149 (MQ=255)
                                                                                                                                             gTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATAGATAACAAAGt         <  3:1256734/149‑1 (MQ=255)
                                                                                                                                                 cgcgGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgc     <  2:521770/149‑1 (MQ=255)
                                                                                                                                                  gcggcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgcg    >  4:205910/1‑149 (MQ=255)
                                                                                                                                                  gcggcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgcg    >  2:88591/1‑149 (MQ=255)
                                                                                                                                                    ggcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGtt  >  3:922933/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTT  >  NC_000913/3101393‑3101689

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: