Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,062,811 C→A 6.7% L189M (CTG→ATG) ‡ yihN → MFS transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,062,8110CA 6.7% 639.1 / 25.9 179L189M (CTG→ATG) ‡yihNMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base C (87/80);  new base A (5/7);  total (92/87)
Fisher's exact test for biased strand distribution p-value = 5.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 180 total aligned reads displayed.

AGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCT  >  NC_000913/4062666‑4062959
                                                                                                                                                 |                                                                                                                                                    
aGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGa                                                                                                                                                   >  4:392600/1‑149 (MQ=255)
     gTGGCATTGTCGATACCGTGGTGGCGCTTCCTGCGTAGGCAGTATTTACCTGGTTTGGCCGTGGCTTATGAGGTATTAAACCAGGCATCTGGTTCTATTCCCTTATTGTGATTACCGTCGGCATGATTATTTTGTTTGTCCTGAATGAc                                                                                                                                              <  4:1174319/149‑1 (MQ=255)
           ttGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAaagaa                                                                                                                                        >  1:1052685/1‑149 (MQ=255)
           ttGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAaagaa                                                                                                                                        <  4:1687886/149‑1 (MQ=255)
           ttGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAaagaa                                                                                                                                        <  3:392600/149‑1 (MQ=255)
             gTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTg                                                                                                                                                    <  3:563312/135‑1 (MQ=255)
              tCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAagag                                                                                                                                     >  3:65518/1‑149 (MQ=255)
                gATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGaa                                                                                                                                                  >  1:1253357/1‑134 (MQ=255)
                        ggtggCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAAt                                                                                                                                                 <  1:1612750/127‑1 (MQ=255)
                         gtggCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAccgtccgt                                                                                                                          <  1:376794/149‑1 (MQ=255)
                          tggCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGtt                                                                                                                         >  2:1646352/1‑149 (MQ=255)
                          tggCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGtt                                                                                                                         >  4:629385/1‑149 (MQ=255)
                            gCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACaaa                                                                                                                                           <  1:1721143/129‑1 (MQ=255)
                            gCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACaaa                                                                                                                                           >  2:1721143/1‑129 (MQ=255)
                            gCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACaaa                                                                                                                                           >  2:1773643/1‑129 (MQ=255)
                             cATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAg                                                                                                                      <  3:629385/149‑1 (MQ=255)
                              aTTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGa                                                                                                                                               >  3:992041/1‑123 (MQ=255)
                               ttttCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGt                                                                                                                    <  4:65518/149‑1 (MQ=255)
                               ttttCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGt                                                                                                                    >  1:1051489/1‑149 (MQ=255)
                                  tCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGa                                                                                                                                               >  2:1187522/1‑119 (MQ=255)
                                           gCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAagaaga                                                                                                        >  3:94571/1‑149 (MQ=255)
                                           gCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAagaaga                                                                                                        >  3:411655/1‑149 (MQ=255)
                                            cAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGAc                                                                                                       <  2:1051489/149‑1 (MQ=255)
                                                      cTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAagag                                                                                                                                     <  1:1410117/109‑1 (MQ=255)
                                                      cTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGa                                                                                             >  1:270812/1‑149 (MQ=255)
                                                        ggTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTTATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGaaa                                                                                           >  1:1610827/1‑149 (MQ=255)
                                                           ttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAac                                                                                        <  4:411655/149‑1 (MQ=255)
                                                               cAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCa                                                                                    >  1:1703549/1‑149 (MQ=255)
                                                               cAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCa                                                                                    <  2:1610827/149‑1 (MQ=255)
                                                                    gCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGa                                                                                                                                               >  4:564408/1‑85 (MQ=255)
                                                                      ttattaGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGa                                                                                                                       >  3:1681303/1‑107 (MQ=255)
                                                                      ttattaGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGa                                                                                                                       <  4:1681303/107‑1 (MQ=255)
                                                                      ttattaGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGAc                                                                             >  3:852270/1‑149 (MQ=255)
                                                                        attaGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAaca                                                                                       >  3:1847322/1‑137 (MQ=255)
                                                                         ttaGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGa                                                                                                    >  2:1721644/1‑123 (MQ=255)
                                                                           aGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGc                                                                                                                                   >  4:537097/1‑90 (MQ=255)
                                                                            ggTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGGTTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGt                                                                       >  2:597507/1‑149 (MQ=255)
                                                                            ggTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGa                                                                                                    <  3:1288234/120‑1 (MQ=255)
                                                                             gTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAagaaga                                                                                                                                      <  1:1499693/85‑1 (MQ=255)
                                                                             gTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGaaa                                                                                           >  2:1435201/1‑128 (MQ=255)
                                                                              ttttAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAagag                                                                                                                                     >  4:1446088/1‑85 (MQ=255)
                                                                              ttttAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTg                                                                      <  4:923617/148‑1 (MQ=255)
                                                                              ttttAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTg                                                                      >  3:923617/1‑148 (MQ=255)
                                                                              ttttAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGc                                                                     >  3:744501/1‑149 (MQ=255)
                                                                              ttttAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGc                                                                     >  3:1891433/1‑149 (MQ=255)
                                                                               tttAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCt                                                                    <  4:847499/149‑1 (MQ=255)
                                                                                  aaaGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGaa                                                                                                                                                  <  2:966298/68‑1 (MQ=255)
                                                                                  aaaGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGaa                                                                                                                                                  >  1:966298/1‑68 (MQ=255)
                                                                                   aaGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGaaa                                                                <  4:1891433/149‑1 (MQ=255)
                                                                                    aGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAacac                                                                                      <  2:816042/126‑1 (MQ=255)
                                                                                    aGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAacac                                                                                      >  1:816042/1‑126 (MQ=255)
                                                                                     gCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCa                                                                                    <  2:1301063/127‑1 (MQ=255)
                                                                                     gCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGa                                                              <  4:981704/149‑1 (MQ=255)
                                                                                          cATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACaaaa                                                         <  4:94571/149‑1 (MQ=255)
                                                                                           aTCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAaagaa                                                                                                          <  4:642348/99‑1 (MQ=255)
                                                                                            tCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGaaa                                                                >  3:653427/1‑140 (MQ=255)
                                                                                             cTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGaaaaa                                                                                         >  4:115761/1‑114 (MQ=255)
                                                                                             cTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGaaaaa                                                                                         <  3:115761/114‑1 (MQ=255)
                                                                                               ggTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATc                                                    <  2:270812/149‑1 (MQ=255)
                                                                                                gTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAagaaga                                                                                                                                      >  1:192953/1‑66 (MQ=255)
                                                                                                       tCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAacac                                                                                      <  4:517283/107‑1 (MQ=255)
                                                                                                       tCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAacac                                                                                      >  3:517283/1‑107 (MQ=255)
                                                                                                        cccTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCa                                                                                    <  2:984788/108‑1 (MQ=255)
                                                                                                        cccTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCa                                                                                    >  1:984788/1‑108 (MQ=255)
                                                                                                         gcTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCg                                          >  2:1612370/2‑149 (MQ=255)
                                                                                                         ccTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGt                                                                       >  3:1043624/1‑120 (MQ=255)
                                                                                                         ccTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGt                                                                       <  4:1043624/120‑1 (MQ=255)
                                                                                                         ccTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGt                                                                       <  2:1702637/120‑1 (MQ=255)
                                                                                                         ccTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGt                                                                       <  4:1388848/120‑1 (MQ=255)
                                                                                                           ttATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACg                                                                                                      <  3:1403536/87‑1 (MQ=255)
                                                                                                             aTTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAccgtccgt                                                                                                                          >  4:762506/1‑65 (MQ=255)
                                                                                                              ttGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTAt                                            >  2:1383917/1‑142 (MQ=255)
                                                                                                                 tgATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAAc                                  <  4:720395/149‑1 (MQ=255)
                                                                                                                    ttGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGa                                                                              <  2:301759/102‑1 (MQ=255)
                                                                                                                    ttGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGa                                                                              >  1:301759/1‑102 (MQ=255)
                                                                                                                         gTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAgg                                                                                                                                    >  1:1510966/1‑43 (MQ=255)
                                                                                                                         gTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAgg                                                                                                                                    <  2:1510966/43‑1 (MQ=255)
                                                                                                                         gTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTcttctt                          <  2:985877/149‑1 (MQ=255)
                                                                                                                           aGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTcttctt                          >  3:1452453/1‑147 (MQ=255)
                                                                                                                           aGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCg                        >  3:134461/1‑149 (MQ=255)
                                                                                                                               attattattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACaaaa                                                         >  3:791638/1‑112 (MQ=255)
                                                                                                                               attattattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCtt                                       <  1:602199/130‑1 (MQ=255)
                                                                                                                               attattattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCtt                                       >  2:602199/1‑130 (MQ=255)
                                                                                                                               attattattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTa                    >  4:1128045/1‑149 (MQ=255)
                                                                                                                                 tattattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTAcg                  >  4:1324804/1‑149 (MQ=255)
                                                                                                                                  attattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCtt                                       <  3:528860/127‑1 (MQ=255)
                                                                                                                                   ttattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTAcgcg                >  1:956467/1‑149 (MQ=255)
                                                                                                                                   ttattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTAcgcg                >  1:1649447/1‑149 (MQ=255)
                                                                                                                                      ttttCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGtt             >  1:173628/1‑149 (MQ=255)
                                                                                                                                      tcttCTTTGTCATTCAGGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGaaaaa                                                                                         >  3:1800551/3‑73 (MQ=37)
                                                                                                                                      tcttCTTTGTCATTCAGGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGaaaaa                                                                                         <  4:1800551/71‑1 (MQ=37)
                                                                                                                                        ttCTTTGTCATTCAGGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCtt                                       <  4:837258/121‑1 (MQ=255)
                                                                                                                                         tCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGttt                     <  2:724763/138‑1 (MQ=255)
                                                                                                                                         tCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGttt                     >  1:724763/1‑138 (MQ=255)
                                                                                                                                           tttGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAg                                                               >  1:1303896/1‑94 (MQ=255)
                                                                                                                                           tttGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAg                                                               <  2:1303896/94‑1 (MQ=255)
                                                                                                                                           tttGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACtg        >  4:452990/1‑149 (MQ=255)
                                                                                                                                               tCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTacc                                  >  2:233006/1‑117 (MQ=255)
                                                                                                                                               tCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGc    >  3:590038/1‑149 (MQ=255)
                                                                                                                                                 cTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCt  <  3:1128045/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCT  >  NC_000913/4062666‑4062959

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