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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,5410GC11.4% 158.1 / 11.5 58P28R (CCG→CGG) yggNputative EcfF
Reads supporting (aligned to +/- strand):  ref base G (32/19);  new base C (5/2);  total (37/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01

CCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGT  >  NC_000913/3101401‑3101682
                                                                                                                                            |                                                                                                                                             
cgtgcGCCGGCTCGCGCGGGGCGGCAGGCGGGGAATATCGCTTACCGTTATACATCCCGTGGCCGTCGGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACgcggc                                                                                                                                             <  1:859477/139‑1 (MQ=255)
ccTTCGCCTGCTCGCCCTGGCCGGCTTCGCTGGGATATTGCTTCCCGTTATACCTCACGTGGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACgcggc                                                                                                                                             <  1:856071/143‑1 (MQ=255)
      cccgcgcgcgCTGGGCGGCAGTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACg                                                                                                                                       <  2:113940/137‑1 (MQ=255)
      ccTTCTCCCGCTGGGCGGCAGTCAGGGAAGATTGCTTACCGTTCTACGTCACGTTGCCGTCTGGCGTGACCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACg                                                                                                                                       <  1:438359/143‑1 (MQ=255)
       cTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGc                                                                                                                                      >  1:1018839/1‑143 (MQ=255)
        tGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGCTTTCGATCACTTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACgat                                                                                                                                     <  2:780413/143‑3 (MQ=255)
          cTCGCGCTGGGCGGCATTCGGGGAATATTGCTTACCGTTATCCATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACGATgt                                                                                                                                   <  1:726688/143‑2 (MQ=255)
          cTCGCGCTGGGCGGCATTCAGGGCATATTGCTTAACGTCATACATCACGTTGCCGACTGGCGTGAGCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGATCTGCGGGCCGACAGCAACACCGTCACACCGCGAGAccgcgc                                                                                                                                   >  2:139355/1‑138 (MQ=255)
          cTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGc                                                                                                                                   >  1:1005204/1‑143 (MQ=255)
           tCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGCTTGCCGTTTTCGCCCTTCACCTGCACGGTTCGCGGGCTGACAATCACATCCTCACGCGGCGCGACGcgccc                                                                                                                                  >  2:242295/1‑139 (MQ=255)
               gcTGGGCGGCACTCAGGGAATATTGCTTACCGTTATCCATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCCGCACTgg                                                                                                                              <  2:592907/143‑1 (MQ=255)
                   ggcggcATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCCCCTCCCCGGTTTGCGAGCTGACAATCACATCGCCACGCGGCGTGACGCTGCCCCGGTAGt                                                                                                                          >  2:216407/1‑143 (MQ=255)
                     cggcATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCg                                                                                                                        >  1:433000/1‑143 (MQ=255)
                         cttCAGGGACTATTGCTTCCCGTTATCCCTCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCTCGCTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGt                                                                                                                    <  1:101989/142‑1 (MQ=255)
                              gggAATATTGCTTACCGTTATACTTCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCg                                                                                                               <  1:376182/143‑1 (MQ=255)
                              gggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCCCATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGGCTGAGCg                                                                                                               >  1:335047/1‑143 (MQ=255)
                              gggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTCGTCGCCGGCGGGACCg                                                                                                               >  2:23614/1‑143 (MQ=255)
                               ggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCgg                                                                                                              >  2:901362/1‑143 (MQ=255)
                                    atctcTTCCCGTTATCCCTCCCGTTGCCGTCTGGCGTGATCCCCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACACTCATATCCTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCAtt                                                                                                         <  1:758202/139‑1 (MQ=255)
                                     ttGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTCCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTg                                                                                                        >  1:99711/1‑143 (MQ=255)
                                             cGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGCGCGGTCATTGCCGTCACt                                                                                                >  2:167935/1‑143 (MQ=255)
                                               ttATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCCGTCATTGCAGTCCCTGa                                                                                              >  1:1089934/1‑143 (MQ=255)
                                                 ataCATCACGTTGCCGTCCGGCGTGATCCCCTGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGaaa                                                                                            <  2:867057/143‑1 (MQ=255)
                                                 ataCATCACGTCGCCGTCTGGCGTGATCACCAGATTGCCGTCTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGaaa                                                                                            <  1:676473/143‑1 (MQ=255)
                                                     aTCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtgc                                                                                        >  1:1040658/1‑143 (MQ=255)
                                                       cACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCCGTCATTGCCGACAATGAAAGtgctg                                                                                      >  2:168796/1‑143 (MQ=255)
                                                       cACGTTGCCGTCTGGCGTGATCACCAGATTCCCGCTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCACTGCCGTCACTGAAAGtgctg                                                                                      >  2:673512/1‑143 (MQ=255)
                                                               cGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGCGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTCGTCGGCGTGAGCGGCCATCGCCGTCACTGCAAGGACTGCCCCcagc                                                                              >  2:331068/1‑143 (MQ=255)
                                                                gTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCTCATCGGCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCcagca                                                                             <  2:1095388/143‑1 (MQ=255)
                                                                gTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTTGTCGGCGTGAGCGGTCATTGCCGTCACTGGAATTGCTGCCGCcagca                                                                             >  1:782127/1‑143 (MQ=255)
                                                                       gTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCAgttt                                                                      >  2:51547/1‑143 (MQ=255)
                                                                         gATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTgc                                                                    <  1:319557/143‑1 (MQ=255)
                                                                             aCCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGcat                                                                >  2:1074646/1‑143 (MQ=255)
                                                                              ccAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGAGcctc                                                               >  2:1056734/1‑143 (MQ=255)
                                                                               caagttGCCGCTCTCGCCCTTCACCGGCACGGGTGGCGGGCTGACAATCACATCGCCACGCGGCGTGACGTTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGcatca                                                              <  1:56852/139‑1 (MQ=255)
                                                                                       cGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTCTGCGCATCATAGTCTTc                                                      >  1:782949/1‑143 (MQ=255)
                                                                                        gTTTTCGCCCTCCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTcc                                                     <  1:464602/143‑1 (MQ=255)
                                                                                        gTTTTCGCCCCTCGCCTGCACGGTTTGCGGGCTGACAGTCACGTCGTCACGCGGCGTGACGCGGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTcc                                                     <  1:357297/143‑1 (MQ=255)
                                                                                         tggtCGCCCTTCACCGGCACGGTTTGCGGGCTGACAATCACATCGCCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTccc                                                    <  1:627787/140‑1 (MQ=255)
                                                                                               cccTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTAACGCGACGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCGGTCACTGAAAGTGCTGCCGCCAGCAGCAATGTGCGGATCATGGACTGCCCTCAAGa                                              >  2:86274/1‑143 (MQ=255)
                                                                                                  ttCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCAGTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATGTTGCGCATCATAGTCTGCCCGCAAGaaaa                                           >  2:853252/1‑143 (MQ=255)
                                                                                                   tCACCGGCACGGTTTGCGGGCTGACAATCGCACCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaa                                          <  1:473350/143‑1 (MQ=255)
                                                                                                          cACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGc                                   >  1:1012732/1‑143 (MQ=255)
                                                                                                            cGGTTGTCGGGCGGACACTGCCATCGTCACGCGGCGGGACGCTGCGCTCGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCa                                 <  1:126280/143‑1 (MQ=255)
                                                                                                               tttGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCa                              >  1:836537/1‑143 (MQ=255)
                                                                                                                   cGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCGGGGTTAGTCGGCGTGAGCGGTCATTGCCGTCCCTGAAAGTGGCGGCGCCGGCGGGGGATTGCGCATCAGAGTCTTCCCGCCACAAGAAAGCAGCCCACCAtttc                          >  2:20746/1‑142 (MQ=255)
                                                                                                                        gggcAATCACATCGGCACGGGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGc                     <  2:811049/140‑1 (MQ=255)
                                                                                                                        ggACAAGCACATCGTCACGCGGCGTGACGCTGCGCTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGc                     <  2:868626/142‑1 (MQ=255)
                                                                                                                            aaTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGCAAGTGCTGCCGCCAGCAGCCTTTCGCGCATCATCGTCTCCCCTCAAGAAAAAATCAGGCCAGCATATGCTGGcccaa                 >  1:39183/1‑141 (MQ=255)
                                                                                                                              tCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCGTAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATATGCTGGCCCCAGa               >  1:613174/1‑143 (MQ=255)
                                                                                                                               cacaTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAt              <  2:1040658/143‑1 (MQ=255)
                                                                                                                               cacaTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGACCTCCCTCCTGAATATATCAGGCCAGCATTTGCTGGCCCCAGAt              >  1:750014/1‑143 (MQ=255)
                                                                                                                                 cagcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGCCCCCAGCTTg            >  2:349857/4‑143 (MQ=255)
                                                                                                                                 caTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTg            >  1:539055/1‑143 (MQ=255)
                                                                                                                                 caTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTg            >  1:17578/1‑143 (MQ=255)
                                                                                                                                  agcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCTTCAAGAAAAAATCAGGCCAGCATTTGCTGGACCCAGATTGa           >  1:527780/3‑143 (MQ=255)
                                                                                                                                  agcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa           >  2:512466/3‑143 (MQ=255)
                                                                                                                                   gcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCCGCATCAGGTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGGTTGAt          >  1:929499/2‑143 (MQ=255)
                                                                                                                                      tCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATGAGGCGAGCACTTGCTGGCCCCAGATTGATAAc       >  2:1066213/1‑143 (MQ=255)
                                                                                                                                      tCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCAGACACCAGATTGATAAc       >  1:146568/1‑143 (MQ=255)
                                                                                                                                         cgcgGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa    <  1:306603/143‑1 (MQ=255)
                                                                                                                                          gcggcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAACAAATAATTCCTGCATTTGCTGGCCCCCGATTGATAACAAAg   >  2:486484/1‑143 (MQ=255)
                                                                                                                                          gcggcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGGTTGATAACAAAg   >  1:754158/1‑143 (MQ=255)
                                                                                                                                          gcggcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGACCCAAAATGGATAACacaa   >  2:253589/1‑140 (MQ=255)
                                                                                                                                          gcggcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGCCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg   >  1:813401/1‑143 (MQ=255)
                                                                                                                                           cggcgTGACGCTGCCCTGGTAGCCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGt  <  1:729158/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                             
CCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGT  >  NC_000913/3101401‑3101682

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: