Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,766,414 G→A 5.3% R166Q (CGA→CAA)  rnlA → CP4‑57 prophage; RNase LS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,766,4140GA 5.3% 669.0 / 16.1 188R166Q (CGA→CAA) rnlACP4‑57 prophage; RNase LS
Reads supporting (aligned to +/- strand):  ref base G (100/78);  new base A (5/5);  total (105/83)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 190 total aligned reads displayed.

TGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAGACTG  >  NC_000913/2766265‑2766557
                                                                                                                                                     |                                                                                                                                               
tGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATATCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCg                                                                                                                                                 >  2:1859038/1‑150 (MQ=255)
      aTCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTTGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCt                                                                                                                                           <  1:2419509/150‑1 (MQ=255)
      aTCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCt                                                                                                                                           >  2:1092713/1‑150 (MQ=255)
         agcagATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTc                                                                                                                                        >  2:841005/1‑150 (MQ=255)
         agcagATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTc                                                                                                                                        >  2:793638/1‑150 (MQ=255)
            agATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATg                                                                                                                                     >  1:1074022/1‑150 (MQ=255)
             gATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGt                                                                                                                                    >  2:1594641/1‑150 (MQ=255)
                 aaTCGCATATAGCGGTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACa                                                                                                                                <  2:1732715/150‑1 (MQ=255)
                 aaTCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACa                                                                                                                                <  1:2496478/150‑1 (MQ=255)
                 aaTCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACa                                                                                                                                <  2:1513794/150‑1 (MQ=255)
                 aaTCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACa                                                                                                                                <  2:2147246/150‑1 (MQ=255)
                 aaTCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACa                                                                                                                                <  2:2147321/150‑1 (MQ=255)
                 aaTCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACa                                                                                                                                <  1:1694157/150‑1 (MQ=255)
                   tCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACaga                                                                                                                              <  2:1089303/150‑1 (MQ=255)
                   tCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACaga                                                                                                                              >  2:2027504/1‑150 (MQ=255)
                   tCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACaga                                                                                                                              >  2:2316058/1‑150 (MQ=255)
                   tCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACaga                                                                                                                              >  2:1832501/1‑150 (MQ=255)
                   tCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACaga                                                                                                                              >  2:2392450/1‑150 (MQ=255)
                   tCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACaga                                                                                                                              <  1:923206/150‑1 (MQ=255)
                   tCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACaga                                                                                                                              >  2:263221/1‑150 (MQ=255)
                   tCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACaga                                                                                                                              <  1:521698/150‑1 (MQ=255)
                      cATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGtt                                                                                                                           <  1:2721419/150‑1 (MQ=255)
                       atatAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGttt                                                                                                                          >  2:1233327/1‑150 (MQ=255)
                       atatAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGttt                                                                                                                          <  1:1914187/150‑1 (MQ=255)
                       atatAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGttt                                                                                                                          >  2:2437845/1‑150 (MQ=255)
                       atatAGAGTTCAGAGAACACACTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGttt                                                                                                                          <  1:1862010/150‑1 (MQ=255)
                         ataGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTc                                                                                                                        >  1:2477609/1‑150 (MQ=255)
                          taGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCa                                                                                                                       >  1:2124182/1‑150 (MQ=255)
                          taGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCa                                                                                                                       >  1:594162/1‑150 (MQ=255)
                               ttCAGAGAACACTCGCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGGTACAGAGTGTTCACGtgt                                                                                                                  >  2:466006/1‑148 (MQ=255)
                                  agagAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAAc                                                                                                               >  2:2714445/1‑150 (MQ=255)
                                    agaACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCAACCGCTGTCATGTTACAGAGTTTTCACGTTTAACtt                                                                                                             >  1:1018410/1‑150 (MQ=255)
                                    agaACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCAACCGCTGTCATGTTACAGAGTTTTCACGTTTAACtt                                                                                                             >  1:2076478/1‑150 (MQ=255)
                                             tcGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTTATCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCt                                                                                                    <  1:1244973/150‑1 (MQ=255)
                                             tcGTAACACTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCt                                                                                                    <  1:937775/150‑1 (MQ=255)
                                                   cGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCATGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGa                                                                                              <  1:2211238/150‑1 (MQ=255)
                                                   cGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGa                                                                                              <  1:40467/150‑1 (MQ=255)
                                                    gCTCATACCGTAGTGTTGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACATGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCATCCCTGCTTGAt                                                                                             >  2:1191180/1‑150 (MQ=255)
                                                    gCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGAt                                                                                             >  1:2573463/1‑150 (MQ=255)
                                                    gCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGAt                                                                                             >  1:2243613/1‑150 (MQ=255)
                                                    gCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGAt                                                                                             >  1:1123981/1‑150 (MQ=255)
                                                      tCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGTGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGAttt                                                                                           >  2:2116898/1‑150 (MQ=255)
                                                      tCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGAttt                                                                                           >  2:1112282/1‑150 (MQ=255)
                                                      tCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGAttt                                                                                           >  2:1572323/1‑150 (MQ=255)
                                                       cATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTa                                                                                          <  1:2313589/150‑1 (MQ=255)
                                                         tACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACa                                                                                        >  1:2691159/1‑150 (MQ=255)
                                                            cGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAggg                                                                                     >  1:810466/1‑150 (MQ=255)
                                                            cGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAggg                                                                                     >  2:2539057/1‑150 (MQ=255)
                                                              tAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGtt                                                                                   <  1:1448911/150‑1 (MQ=255)
                                                                     gAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGaa                                                                            >  1:1016376/1‑150 (MQ=255)
                                                                        aaTTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGt                                                                         >  1:1515365/1‑150 (MQ=255)
                                                                          ttattTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGc                                                                       >  2:1847759/1‑150 (MQ=255)
                                                                          ttattTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGc                                                                       >  1:1284538/1‑150 (MQ=255)
                                                                            attTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTa                                                                     <  1:866494/150‑1 (MQ=255)
                                                                            attTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTa                                                                     <  2:1189168/150‑1 (MQ=255)
                                                                             tttCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTat                                                                    <  2:956954/150‑2 (MQ=255)
                                                                             tttCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTaa                                                                    <  2:2668417/150‑1 (MQ=255)
                                                                             tttCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTaa                                                                    <  2:1803076/150‑1 (MQ=255)
                                                                             tttCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTaa                                                                    <  2:1971179/150‑1 (MQ=255)
                                                                             tttCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTaa                                                                    <  1:196453/150‑1 (MQ=255)
                                                                                   ccaGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGcc                                                              >  2:2244055/1‑150 (MQ=255)
                                                                                   ccaGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGcc                                                              >  2:614399/1‑150 (MQ=255)
                                                                                       cTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCagga                                                          >  2:521419/1‑150 (MQ=255)
                                                                                             ggACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATgg                                                    >  1:347289/1‑150 (MQ=255)
                                                                                              gACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGt                                                   >  1:60411/1‑150 (MQ=255)
                                                                                                cGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTaa                                                 <  2:1515365/150‑1 (MQ=255)
                                                                                                      gACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTaa                                           <  1:1775016/150‑1 (MQ=255)
                                                                                                            gAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGt                                     >  2:700089/1‑150 (MQ=255)
                                                                                                            gAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGt                                     >  2:1726415/1‑150 (MQ=255)
                                                                                                                   cATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACag                              <  1:1703423/150‑1 (MQ=255)
                                                                                                                     tatCACAACAGGTAAGCTCCAGATTCAGGGCCGACCTCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTTATATTGTTCAACagga                            <  2:1652479/150‑1 (MQ=255)
                                                                                                                     tatCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACagga                            >  2:192410/1‑150 (MQ=255)
                                                                                                                     tatCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACagga                            >  2:32291/1‑150 (MQ=255)
                                                                                                                     tatCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACagga                            >  1:2196650/1‑150 (MQ=255)
                                                                                                                         acaacaGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGtt                        >  1:2030020/1‑150 (MQ=255)
                                                                                                                         acaacaGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGtt                        >  1:2588569/1‑150 (MQ=255)
                                                                                                                                gTAAGCTCCAGATTCAGGGCCAACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCa                 <  2:1018410/150‑1 (MQ=255)
                                                                                                                                gTAAGCTCCAGATTCAGGGCCAACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCa                 <  2:2076478/150‑1 (MQ=255)
                                                                                                                                gTAAGCTCCAGATTCAGGGCCAACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCa                 <  2:1369274/150‑1 (MQ=255)
                                                                                                                                 tAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCac                >  2:1172826/1‑150 (MQ=255)
                                                                                                                                        cAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTg         <  2:1414256/150‑1 (MQ=255)
                                                                                                                                        cAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTg         <  2:2263880/150‑1 (MQ=255)
                                                                                                                                          gATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCa       <  1:317697/150‑1 (MQ=255)
                                                                                                                                          gATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCa       <  1:1572323/150‑1 (MQ=255)
                                                                                                                                          gATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCa       <  1:1112269/150‑1 (MQ=255)
                                                                                                                                          gATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCa       <  2:1288666/150‑1 (MQ=255)
                                                                                                                                          gATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCa       <  2:1334315/150‑1 (MQ=255)
                                                                                                                                          gATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCa       <  1:2116898/150‑1 (MQ=255)
                                                                                                                                           aTTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAg      >  1:246853/1‑150 (MQ=255)
                                                                                                                                           aTTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAg      >  1:727194/1‑150 (MQ=255)
                                                                                                                                           aTTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAg      >  1:727215/1‑150 (MQ=255)
                                                                                                                                           aTTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAg      >  1:254000/1‑150 (MQ=255)
                                                                                                                                            ttCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCATGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAGa     <  1:498286/150‑1 (MQ=255)
                                                                                                                                            ttCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAGa     <  2:2738094/150‑1 (MQ=255)
                                                                                                                                            ttCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAGa     <  1:897090/150‑1 (MQ=255)
                                                                                                                                              cAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAGACt   >  1:856405/1‑150 (MQ=255)
                                                                                                                                              cAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAGACt   <  1:840541/150‑1 (MQ=255)
                                                                                                                                              cAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAGACt   >  2:1096403/1‑150 (MQ=255)
                                                                                                                                              cAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAGACt   >  2:1738369/1‑150 (MQ=255)
                                                                                                                                               aGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAGACTg  <  2:1284538/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATTTCCACCAGCTATCAGGACGAATTGACTGTGAGCCTGCATATCACAACAGGTAAGCTCCAGATTCAGGGCCGACCGCTGTCATGTTACAGAGTTTTCACGTTTAACTTGGCAGCCCTGCTTGATTTACAGGGTTTGGAGAAAGTGCTAATCCGCCAGGAGGATGGTAAAGCTAATATTGTTCAACAGGAGGTTGCCCGCACTTACTTGCAGACTG  >  NC_000913/2766265‑2766557

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