Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,967,449 C→T 19.4% 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,4490CT19.4% 301.8 / 61.5 140L345F (CTC→TTC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (60/52);  new base T (16/11);  total (77/63)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

Only 100 of 147 total aligned reads displayed.

TCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGA  >  NC_000913/3967301‑3967597
                                                                                                                                                    |                                                                                                                                                    
tCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACct                                                                                                                                                     >  2:1142485/1‑150 (MQ=255)
 cGTCCGTAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACttc                                                                                                                                                    <  2:657837/150‑1 (MQ=255)
   tCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACctctc                                                                                                                                                  <  1:1179157/150‑1 (MQ=255)
    ccGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTTCAGGCAACATGGAACTGCACctctct                                                                                                                                                 >  1:492371/1‑150 (MQ=255)
    ccGAAACGCGTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGCACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACctctct                                                                                                                                                 <  1:383211/150‑1 (MQ=255)
        aaCGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTa                                                                                                                                             <  2:984404/150‑1 (MQ=255)
        aaCGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTa                                                                                                                                             <  1:1662293/150‑1 (MQ=255)
            cttcttTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAACGGTTCAGGCAACATGGACCGGGACTTCTCTCGTAAGAt                                                                                                                                         >  1:266642/1‑150 (MQ=255)
               cttTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGc                                                                                                                                      >  2:1072027/1‑150 (MQ=255)
                 ttGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTg                                                                                                                                    >  1:813702/1‑150 (MQ=255)
                 ttGGTGCGGCGCGTAACGTGGAAGAGGGCGACAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTg                                                                                                                                    <  1:250718/150‑1 (MQ=255)
                   ggTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGaa                                                                                                                                  >  2:2012733/1‑150 (MQ=255)
                    gTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGACCTGCACCTCTCTCGTAAGATCGCTGaaa                                                                                                                                 >  1:677328/1‑150 (MQ=255)
                      gcggcgCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGaaaaa                                                                                                                               >  1:510596/1‑150 (MQ=255)
                      gcggcgCGTAACGTGGAAGAGGGCGGCAGCCGGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACTAAGAGTTTAAAGGTACAGGCAACATGGAACTGGACCTCTCTCGTAAGATCGCTGaaaaa                                                                                                                               <  1:1823759/150‑1 (MQ=255)
                        ggcgCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAAcg                                                                                                                             <  1:417607/150‑1 (MQ=255)
                        ggcgCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAAcg                                                                                                                             <  2:322895/150‑1 (MQ=255)
                        ggcgCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTGCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAAcg                                                                                                                             >  2:329090/1‑150 (MQ=255)
                        ggcgCGTAACGTGGAAGAGGGCGGCAGCCGGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAAcg                                                                                                                             <  2:1906211/150‑1 (MQ=255)
                         gcgcGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAAcgc                                                                                                                            <  1:1329491/150‑1 (MQ=255)
                         gcgcGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAAcgc                                                                                                                            >  2:1925631/1‑150 (MQ=255)
                         gcgcGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAAcgc                                                                                                                            <  1:457038/150‑1 (MQ=255)
                          cgcgTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAAcgcg                                                                                                                           <  1:318686/150‑1 (MQ=255)
                            cgTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTc                                                                                                                         <  2:591708/150‑1 (MQ=255)
                                  gTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCg                                                                                                                   >  1:1662899/1‑150 (MQ=255)
                                    ggAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGc                                                                                                                 >  2:1531604/1‑150 (MQ=255)
                                    ggAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGc                                                                                                                 <  1:1713668/150‑1 (MQ=255)
                                          gggCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGCTACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGa                                                                                                           >  2:1947641/1‑150 (MQ=255)
                                          gggCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGa                                                                                                           >  2:574602/1‑150 (MQ=255)
                                          gggCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGa                                                                                                           >  1:1327120/1‑150 (MQ=255)
                                           ggcggcAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGAc                                                                                                          >  1:1034626/1‑150 (MQ=255)
                                           ggcggcAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGCTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGAc                                                                                                          <  2:1122693/150‑1 (MQ=255)
                                                 aGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGATCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAAc                                                                                                    <  1:1832123/150‑1 (MQ=255)
                                                      gACCATTATCGCGACGGAGCTTATCGATACCGGTTCTAATATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTc                                                                                               <  2:255151/150‑1 (MQ=255)
                                                        ccATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTg                                                                                             <  2:1542829/150‑1 (MQ=255)
                                                        ccATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTg                                                                                             <  1:737700/150‑1 (MQ=255)
                                                         cATTATCGCGTCGGGGATTATCGATACCGGTTCTAAAATTGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAATAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTgg                                                                                            <  2:1303572/150‑1 (MQ=255)
                                                         cATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTgg                                                                                            >  2:121650/1‑150 (MQ=255)
                                                          aTTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGt                                                                                           >  1:707136/1‑150 (MQ=255)
                                                             aTCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTAcc                                                                                        <  2:453762/150‑1 (MQ=255)
                                                                       cgcTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGCCTACAACCGTTCTGGTACCCGTAaagaag                                                                              >  2:1747722/1‑150 (MQ=255)
                                                                         cTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAagag                                                                            <  1:553811/150‑1 (MQ=255)
                                                                          ttttCGATACCGGTTCTAAAATGGACGATGTTATCTACGAATAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgc                                                                           <  2:67734/147‑1 (MQ=255)
                                                                           tATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgct                                                                          >  2:1586409/1‑150 (MQ=255)
                                                                             tagaaaCCGGAGCTAAAATTGAAGAAGTTATATACTAAGATTTTAAAGGTACTGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgctgc                                                                        <  2:176476/145‑1 (MQ=255)
                                                                             tCGATACCGGTTCTAAAATGGTCGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGTTACCCGTAAAGAAGAgctgc                                                                        <  2:198616/150‑1 (MQ=255)
                                                                             tCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgctgc                                                                        >  1:1722677/1‑150 (MQ=255)
                                                                              cGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgctgct                                                                       >  2:666412/1‑150 (MQ=255)
                                                                              cGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgctgcg                                                                       >  2:899977/1‑149 (MQ=255)
                                                                               gATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTc                                                                      >  1:1434049/1‑150 (MQ=255)
                                                                                aTACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCa                                                                     <  1:767374/150‑1 (MQ=255)
                                                                                aTACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCa                                                                     <  2:61271/150‑1 (MQ=255)
                                                                                 aaCCGGTTCTAAAATGGAAGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGTTACTGCACCTCTCTCGTAAGAGCGCTGAAAAACGCGTGTTCCCGGCTATCGACTACAACCTTTCTGGTACCCGTAAAGAAGAGCTGCTCAc                                                                    <  1:1295222/149‑1 (MQ=255)
                                                                                    cGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGAc                                                                 >  1:146274/1‑150 (MQ=255)
                                                                                    cGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGGTCTGGTACCCGTAAAGAAGAGCTGCTCACGAc                                                                 >  2:419449/1‑150 (MQ=255)
                                                                                      gTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTc                                                               <  2:1034626/150‑1 (MQ=255)
                                                                                        tatAAAATGGACTAAGTTATCTACGAATAGTTTAAAGTTACAGGCAACATGGAAATGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTATTCCCGGCTATCGAATACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAg                                                             <  1:993793/148‑1 (MQ=255)
                                                                                        tCTAAATTGGACGAAGTAATCTACGAAGAGTTTAAAGGTACAGGCACCATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAg                                                             <  2:1616140/150‑1 (MQ=255)
                                                                                         cTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCTCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAgg                                                            <  2:83192/150‑1 (MQ=255)
                                                                                          tAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACATCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGga                                                           >  2:1157417/1‑150 (MQ=255)
                                                                                           aaaaTGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaa                                                          <  1:1505729/150‑1 (MQ=255)
                                                                                            aaaTGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaag                                                         >  2:1610143/1‑150 (MQ=255)
                                                                                             aaTGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACATCCGTTCTGGTACCCGTAAAGCAGAGCTGCTCACGACTCAGgaaga                                                        >  2:845899/1‑150 (MQ=255)
                                                                                                 gACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTg                                                    >  2:1170917/1‑150 (MQ=255)
                                                                                                 gACGAACTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTg                                                    >  2:1499419/1‑150 (MQ=255)
                                                                                                  ccGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGTAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAATGAAGAGCTGCTCACGACTCAGGAAGAACTGc                                                   <  1:1546034/149‑1 (MQ=255)
                                                                                                  aCGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGc                                                   <  2:1331384/150‑1 (MQ=255)
                                                                                                   cGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCCTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCa                                                  >  2:280140/1‑150 (MQ=255)
                                                                                                     aaGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGa                                                >  2:1873978/1‑150 (MQ=255)
                                                                                                      aGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGAGTCTTCCCGGCTATCGAATACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGaa                                               >  2:516186/1‑150 (MQ=255)
                                                                                                      aGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAATAGCTGCTCACGACTCAGGAAGAACTGCAGaa                                               <  2:1327120/150‑1 (MQ=255)
                                                                                                         tATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAt                                            >  2:867869/1‑150 (MQ=255)
                                                                                                           tCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCCCGACGCAGGAAGCACTGCAGACAAtgg                                          >  2:1281055/1‑149 (MQ=255)
                                                                                                            cTACGAAGAGTTTAAAGGTACACGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAtgtg                                         <  2:806818/150‑1 (MQ=255)
                                                                                                               cGAAGAGTTTAAAGTTACAGGAAACATTGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGATGAACTGCAGAAAATGTGGAt                                      <  2:158732/150‑1 (MQ=255)
                                                                                                                 aaGAGTTTAAAGGTACAGCCAAAATGTAACTGCACTTCTCTCGTAAGATCGCGGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATcc                                    <  2:83706/150‑1 (MQ=255)
                                                                                                                    agTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCTGAAAATGTGGATCCTgc                                 >  1:1076764/1‑150 (MQ=255)
                                                                                                                    agTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTgc                                 <  1:574602/150‑1 (MQ=255)
                                                                                                                     gTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTgcg                                >  1:530866/1‑150 (MQ=255)
                                                                                                                       ttAATGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCCCa                              <  2:684426/150‑1 (MQ=255)
                                                                                                                       ttAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCa                              >  1:200930/1‑150 (MQ=255)
                                                                                                                        tAAAGGTACAGGCAACATGGAACTGCACCTCTCTAGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGATCTGCTCACGACTCAGGAAGAACTGCATAAAATGTGGATCCTGCGCaa                             >  2:1569803/1‑150 (MQ=255)
                                                                                                                         aaaGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGACCTGCAGAAAATGTGGATCCTGCGCaaa                            >  1:163960/1‑150 (MQ=255)
                                                                                                                         aaaGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCaaa                            >  2:798600/1‑150 (MQ=255)
                                                                                                                           aGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAAt                          <  2:1989854/150‑1 (MQ=255)
                                                                                                                           aGGTACAGGCAACATGGAACTGCACCGCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGCGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAAt                          <  1:1989680/150‑1 (MQ=255)
                                                                                                                            ggTACCGGCAACATGTAACTGCCCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATc                         <  2:424460/150‑1 (MQ=255)
                                                                                                                            ggTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATc                         >  1:1838255/1‑150 (MQ=255)
                                                                                                                              tACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCAt                       <  1:67691/150‑1 (MQ=255)
                                                                                                                              tACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCAt                       >  2:1531384/1‑150 (MQ=255)
                                                                                                                                cAGGCAACATGGAACTGCACCTCTCTGGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACGGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTc                     <  2:1281965/150‑1 (MQ=255)
                                                                                                                                cAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGCGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTc                     >  2:1255177/1‑150 (MQ=255)
                                                                                                                                  ggCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCAc                   >  2:485923/1‑150 (MQ=255)
                                                                                                                                        aTGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATg             >  2:1878016/1‑150 (MQ=255)
                                                                                                                                           taaCTGCACCTCTCTCGTAAGATCGCTCAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGc          <  1:1518348/149‑1 (MQ=255)
                                                                                                                                           gAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGc          >  2:576771/1‑150 (MQ=255)
                                                                                                                                           gAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGc          >  2:1491638/1‑150 (MQ=255)
                                                                                                                                                gCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAAt     >  1:868444/1‑150 (MQ=255)
                                                                                                                                                 cACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATc    <  2:2061325/150‑1 (MQ=255)
                                                                                                                                                   ccTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCACGGCTATCGACTACAACCTTTCTGATACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACAAGATGGGCGAAATCGa  >  1:114856/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGA  >  NC_000913/3967301‑3967597

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑