Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,967,449 C→T 18.9% L345F (CTC→TTC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,967,4490CT18.9% 889.6 / 190.1 375L345F (CTC→TTC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (140/163);  new base T (34/37);  total (174/201)
Fisher's exact test for biased strand distribution p-value = 8.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 396 total aligned reads displayed.

ATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGAT  >  NC_000913/3967300‑3967598
                                                                                                                                                     |                                                                                                                                                     
aTCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCAcc                                                                                                                                                       <  1:1252043/150‑1 (MQ=255)
 tCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACct                                                                                                                                                      <  1:3436670/150‑1 (MQ=255)
      cGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCtctctc                                                                                                                                                 <  1:1194327/150‑1 (MQ=255)
        aaaCGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGt                                                                                                                                               <  2:737107/150‑1 (MQ=255)
         aaCGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGTTACCGGTTCTATAATGGACGTAGTTATCTACGAAGAGTTTAATGGTACAGGCAACATGGAACTGCACTTCTCTCGTa                                                                                                                                              >  2:395658/1‑150 (MQ=255)
             cttcttTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGAt                                                                                                                                          <  1:3307185/150‑1 (MQ=255)
              ttcttTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATc                                                                                                                                         >  2:5117617/1‑150 (MQ=255)
              ttcttTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATc                                                                                                                                         >  1:1661501/1‑150 (MQ=255)
               tcttTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCg                                                                                                                                        >  1:5132674/1‑150 (MQ=255)
                 tttGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCt                                                                                                                                      <  1:5943376/150‑1 (MQ=255)
                 tttGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCt                                                                                                                                      <  2:2079025/150‑1 (MQ=255)
                 tttGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCt                                                                                                                                      <  1:1783641/150‑1 (MQ=255)
                 tttGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCt                                                                                                                                      <  1:70178/150‑1 (MQ=255)
                 tttGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCt                                                                                                                                      <  2:2211511/150‑1 (MQ=255)
                 ttgggTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCt                                                                                                                                      <  1:5806722/147‑1 (MQ=255)
                    ggTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGaa                                                                                                                                   >  1:2730932/1‑150 (MQ=255)
                          gcgcTTATCGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAAcgc                                                                                                                             <  2:3618872/150‑1 (MQ=255)
                          gcgcGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAAcgc                                                                                                                             <  1:1464939/150‑1 (MQ=255)
                             cgTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTc                                                                                                                          >  2:3636763/1‑150 (MQ=255)
                             cgTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTc                                                                                                                          >  2:5750313/1‑150 (MQ=255)
                              gTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCt                                                                                                                         <  1:2623142/150‑1 (MQ=255)
                                  cGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTccc                                                                                                                     <  1:341175/150‑1 (MQ=255)
                                    tGGAAGAGGTCGGCAGCCTGACCATTATCGAGACGTAGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTAAGAAGAGTTTAAAGGTACAGGAAACATGGATCTGCACATCTCTCTTAAGATCGCTGAAAAACGCGTCTTCCCgg                                                                                                                   <  2:1820329/150‑1 (MQ=255)
                                    tGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCgg                                                                                                                   >  1:4093877/1‑150 (MQ=255)
                                        agagGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTAt                                                                                                               <  2:2288870/150‑1 (MQ=255)
                                        agagGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTAt                                                                                                               >  1:3969852/1‑150 (MQ=255)
                                           gggCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGa                                                                                                            <  1:6006281/150‑1 (MQ=255)
                                              cggcAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTa                                                                                                         >  1:833270/1‑150 (MQ=255)
                                              cggcAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTa                                                                                                         <  2:2317624/150‑1 (MQ=255)
                                              cggcAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTa                                                                                                         <  2:4873046/150‑1 (MQ=255)
                                              cggcAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTa                                                                                                         >  1:4051464/1‑150 (MQ=255)
                                                gcAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACa                                                                                                       <  1:3945912/150‑1 (MQ=255)
                                                   gCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAAcc                                                                                                    <  1:4364452/150‑1 (MQ=255)
                                                    ccTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCg                                                                                                   >  1:5517047/1‑150 (MQ=255)
                                                     cTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGt                                                                                                  >  2:628442/1‑150 (MQ=255)
                                                      tGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCACCATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGtt                                                                                                 >  2:3068603/1‑150 (MQ=255)
                                                        aCCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCt                                                                                               <  2:408618/150‑1 (MQ=255)
                                                        aCCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCt                                                                                               >  1:2976246/1‑150 (MQ=255)
                                                        aCCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCt                                                                                               <  2:3329193/150‑1 (MQ=255)
                                                         ccATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTg                                                                                              <  1:466623/150‑1 (MQ=255)
                                                          cATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTgg                                                                                             >  2:5154153/1‑150 (MQ=255)
                                                          cATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTgg                                                                                             <  2:2244427/150‑1 (MQ=255)
                                                               tCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTAccc                                                                                        >  2:3456721/1‑150 (MQ=255)
                                                               tCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTAccc                                                                                        <  1:3196612/150‑1 (MQ=255)
                                                               tCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTAccc                                                                                        <  2:1622852/150‑1 (MQ=255)
                                                                  cgACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGATGAAAAACGCGTCTTCCCGGCTATCGACTACAAACGTTCTTGTACCCGTa                                                                                     >  1:524115/1‑150 (MQ=255)
                                                                    aCGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTaaa                                                                                   <  2:1017188/150‑1 (MQ=255)
                                                                         gcTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAagaaga                                                                              <  1:3613732/150‑1 (MQ=255)
                                                                         gcTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAagaaga                                                                              <  1:4555821/150‑1 (MQ=255)
                                                                         gcTGATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAagaaga                                                                              <  1:2525807/150‑1 (MQ=255)
                                                                          cTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAagag                                                                             <  2:1075385/150‑1 (MQ=255)
                                                                           ttATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgc                                                                            >  2:4423440/1‑150 (MQ=255)
                                                                             aTCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgctg                                                                          >  2:3284562/1‑150 (MQ=255)
                                                                                gATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTc                                                                       >  2:5614727/1‑150 (MQ=255)
                                                                                    ccGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGa                                                                   <  2:2709869/150‑1 (MQ=255)
                                                                                    ccGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGa                                                                   <  2:3853275/150‑1 (MQ=255)
                                                                                     cGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGAc                                                                  >  1:1516585/1‑150 (MQ=255)
                                                                                      ggTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACt                                                                 >  1:137471/1‑150 (MQ=255)
                                                                                        ttCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCa                                                               >  2:5964961/1‑150 (MQ=255)
                                                                                          atAAAATGGACGAATTTATCTACTAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAATATCGCTTAAAAACGCGTCTTCCCGGCTATCGACTACACCCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAgg                                                             <  1:4854137/149‑1 (MQ=255)
                                                                                            aaaaTGGACGAAGTTATCTACTAATAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaa                                                           <  1:3124452/150‑1 (MQ=255)
                                                                                            aaaaTGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaa                                                           <  2:1239512/150‑1 (MQ=255)
                                                                                               aTGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaagaa                                                        >  2:4172837/1‑150 (MQ=255)
                                                                                                tGGACGAAGTTATCTACTAAGTGTTTAAAGGTACAGGCAACATGGATCTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAAc                                                       <  1:4879194/150‑1 (MQ=255)
                                                                                                  gACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTg                                                     >  2:1985657/1‑150 (MQ=255)
                                                                                                  gACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTg                                                     >  1:2481635/1‑150 (MQ=255)
                                                                                                          tATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGTACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAt                                             >  2:2574821/1‑150 (MQ=255)
                                                                                                           aTCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAtg                                            >  1:2449745/1‑150 (MQ=255)
                                                                                                            tCTACGAAGAGTTTAAAGGTACTGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGATATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAtgt                                           <  1:4911072/150‑1 (MQ=255)
                                                                                                            tCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTATTATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGGTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAtgt                                           <  2:4601985/150‑1 (MQ=255)
                                                                                                             cTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAtgtg                                          <  1:2071501/150‑1 (MQ=255)
                                                                                                              tACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTgg                                         >  1:225158/1‑150 (MQ=255)
                                                                                                               aCGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGa                                        >  1:1074873/1‑150 (MQ=255)
                                                                                                                  aaGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATcc                                     <  1:4682350/150‑1 (MQ=255)
                                                                                                                     agTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTgc                                  <  2:5537592/150‑1 (MQ=255)
                                                                                                                       tttAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTgcgc                                <  2:5913711/150‑1 (MQ=255)
                                                                                                                       tttAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTgcgc                                >  2:2157951/1‑150 (MQ=255)
                                                                                                                       tttAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTgcgc                                >  2:4747386/1‑150 (MQ=255)
                                                                                                                         tAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCaa                              >  2:3217210/1‑150 (MQ=255)
                                                                                                                          aaaGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCaaa                             >  2:1381131/1‑150 (MQ=255)
                                                                                                                             ggTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATc                          <  1:3837912/150‑1 (MQ=255)
                                                                                                                              gTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCa                         >  1:3166636/1‑150 (MQ=255)
                                                                                                                              gTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCa                         >  2:2868615/1‑150 (MQ=255)
                                                                                                                              gTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCa                         >  2:2446428/1‑150 (MQ=255)
                                                                                                                              gTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCa                         >  2:2012577/1‑150 (MQ=255)
                                                                                                                              gTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCa                         >  1:909916/1‑150 (MQ=255)
                                                                                                                                   ggCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCAc                    <  2:5288781/150‑1 (MQ=255)
                                                                                                                                    gCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCAcc                   <  1:1009555/150‑1 (MQ=255)
                                                                                                                                     cAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCAccc                  >  1:3471337/1‑150 (MQ=255)
                                                                                                                                     cAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAGGTGGATCCTGCGCAAAATCATTCAccc                  >  1:3723231/1‑150 (MQ=255)
                                                                                                                                       aCATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGa                <  1:664328/150‑1 (MQ=255)
                                                                                                                                        cATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGAt               <  1:443830/150‑1 (MQ=255)
                                                                                                                                          tGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACt                                                                                                          <  2:3758323/57‑1 (MQ=255)
                                                                                                                                            gAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGc           <  1:326832/150‑1 (MQ=255)
                                                                                                                                            gAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGc           >  1:2312975/1‑150 (MQ=255)
                                                                                                                                             aaCTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCg          >  1:5711731/1‑150 (MQ=255)
                                                                                                                                             aaCTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCg          >  2:835892/1‑150 (MQ=255)
                                                                                                                                               cTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGaa        <  1:1597243/150‑1 (MQ=255)
                                                                                                                                                  cACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATc     <  1:5329413/150‑1 (MQ=255)
                                                                                                                                                     ctctctCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGAt  <  2:4144090/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGAT  >  NC_000913/3967300‑3967598

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