Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,018,983 A→T 6.2% intergenic (+46/+30) matP → / ← ompA Ter macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,018,9830AT 6.2% 313.1 / 4.6 114intergenic (+46/+30)matP/ompATer macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base A (53/53);  new base T (3/5);  total (56/58)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 120 total aligned reads displayed.

AAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGT  >  NC_000913/1018834‑1019126
                                                                                                                                                      |                                                                                                                                               
aaaCGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCgca                                                                                                                                                 >  1:1508106/1‑150 (MQ=255)
 aaCGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGAAGTATTGGCAGGCATTGAGGGTAAAGAAATAAATCAC‑AAAGTAAAGTTAAATTAAATAAAGCAAAAAAAAACCCCgcac                                                                                                                                                >  2:1565954/1‑149 (MQ=255)
   cGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGC‑AAAAAAACCCCGcagcag                                                                                                                                             >  1:1277134/1‑150 (MQ=255)
    gATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGcagcag                                                                                                                                             >  2:1557262/1‑150 (MQ=255)
     aTTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCagcagc                                                                                                                                            >  1:792126/1‑150 (MQ=255)
         ttCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCgggg                                                                                                                                        <  1:1494293/150‑1 (MQ=255)
         ttCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCgggg                                                                                                                                        <  1:258796/150‑1 (MQ=255)
         ttCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCgggg                                                                                                                                        <  1:42067/150‑1 (MQ=255)
           cAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGGtt                                                                                                                                      <  2:525300/150‑1 (MQ=255)
            aGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGTGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGGttt                                                                                                                                     <  1:307972/150‑1 (MQ=255)
              cTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGttttt                                                                                                                                   >  1:65197/1‑150 (MQ=255)
               tGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTc                                                                                                                                  <  1:1205183/150‑1 (MQ=255)
                aaatGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATTTCTTCTCTTAATCAGGATCTGCAGGCATTGCTGGGTAAGTAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCGCGCAGCAGCGGGGTTTTTCt                                                                                                                                 <  2:877558/147‑1 (MQ=255)
                  ttGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAc                                                                                                                               <  2:806178/150‑1 (MQ=255)
                   tGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGGttttttttgc                                                                                                                              <  1:1088635/150‑6 (MQ=255)
                   tGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAcc                                                                                                                              <  2:1649562/150‑1 (MQ=255)
                   tGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCACGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAcc                                                                                                                              <  1:1041561/150‑1 (MQ=255)
                        aTGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACg                                                                                                                         <  1:1557262/150‑1 (MQ=255)
                           cGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACgaga                                                                                                                      >  2:461153/1‑150 (MQ=255)
                             gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAAc                                                                                                                    >  2:1216110/1‑150 (MQ=255)
                               aaaCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACtt                                                                                                                  >  1:1211154/1‑150 (MQ=255)
                               aaaCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGAATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACtt                                                                                                                  >  1:808394/1‑150 (MQ=255)
                                 aCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                >  1:653442/1‑150 (MQ=255)
                                  cAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAg                                                                                                               <  2:1718329/150‑1 (MQ=255)
                                  cAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAg                                                                                                               >  1:1189644/1‑150 (MQ=255)
                                  cAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAg                                                                                                               >  1:1639580/1‑150 (MQ=255)
                                   aaaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                              <  1:1465977/150‑1 (MQ=255)
                                   aaaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                              >  2:461857/1‑150 (MQ=255)
                                    aaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGcc                                                                                                             >  2:469118/1‑150 (MQ=255)
                                         aaaTACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCgg                                                                                                        >  1:1207364/1‑150 (MQ=255)
                                         aaaTACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCgg                                                                                                        <  2:521311/150‑1 (MQ=255)
                                          aaTACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACAAGACGAGAACTTAAGCATGCgga                                                                                                       >  1:168593/1‑149 (MQ=255)
                                              cgcgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAg                                                                                                   <  2:748416/150‑1 (MQ=255)
                                               gcgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGt                                                                                                  <  1:452616/150‑1 (MQ=255)
                                                cgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGtt                                                                                                 >  2:160680/1‑150 (MQ=255)
                                                cgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGtt                                                                                                 >  2:414403/1‑150 (MQ=255)
                                                 gAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTa                                                                                                <  2:1333378/150‑1 (MQ=255)
                                                   aTAAAATGGCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACa                                                                                              <  1:825408/150‑1 (MQ=255)
                                                     aaaaTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAAc                                                                                            <  1:622107/150‑1 (MQ=255)
                                                     aaaaTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGATCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAACAGACGAGAACTTAAGCATGCGGCTGAGTTACAAc                                                                                            >  1:1690845/1‑150 (MQ=255)
                                                      aaaTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACg                                                                                           >  1:1435967/1‑150 (MQ=255)
                                                        aTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACGACGtg                                                                                         >  2:1631660/1‑149 (MQ=255)
                                                           tcttctCTGAAGCAGGTTCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CATTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                      <  2:871060/150‑1 (MQ=255)
                                                             ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGTAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGa                                                                                    <  1:438763/150‑1 (MQ=255)
                                                             ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGa                                                                                    >  1:1478441/1‑150 (MQ=255)
                                                               ctctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATa                                                                                  >  1:1503494/1‑150 (MQ=255)
                                                                 ctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGa                                                                                                    <  1:1236138/130‑1 (MQ=255)
                                                                 ctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGa                                                                                                    >  2:1236138/1‑130 (MQ=255)
                                                                    aaGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCttt                                                                             <  1:801351/150‑1 (MQ=255)
                                                                     aGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTa                                                                            >  2:307080/1‑150 (MQ=255)
                                                                      gCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTAAAGCCTGCGGCTGAGTTACAACGTCTTTGATACATTTaa                                                                           >  2:146866/1‑150 (MQ=255)
                                                                         ggATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACtt                                                                        >  2:808591/1‑150 (MQ=255)
                                                                          gATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTc                                                                       >  1:194846/1‑150 (MQ=255)
                                                                              tGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtc                                                                   <  1:1015058/150‑1 (MQ=255)
                                                                              tGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtc                                                                   <  1:1588091/150‑1 (MQ=255)
                                                                              tGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtc                                                                   >  1:756211/1‑150 (MQ=255)
                                                                               gCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGAAGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtct                                                                  >  1:885099/1‑150 (MQ=255)
                                                                                cAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCACCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtctc                                                                 >  2:136607/1‑150 (MQ=255)
                                                                                 aGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATctct                                                                <  2:653442/150‑1 (MQ=255)
                                                                                    cATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAcg                                                             <  1:893014/150‑1 (MQ=255)
                                                                                     aTTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAcgc                                                            <  2:252676/150‑1 (MQ=255)
                                                                                         cTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcgacg                                                        <  2:792126/150‑1 (MQ=255)
                                                                                            ggTAAGGAATAACTGTCGAAAGT‑CATTTCAATTTACTAAAGGAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATc                                                     <  2:817294/150‑1 (MQ=255)
                                                                                            ggTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGGAGCGGGGTTTTTCTACCAGAGGAGAACTTAAGCCTGCGGCGGAGTTACAACGTCTTTGATTCCTTTAACTTCGATCTCTACGCGACGATc                                                     >  1:1132817/1‑150 (MQ=255)
                                                                                            ggTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCttt                                                                             >  1:1320561/1‑126 (MQ=255)
                                                                                            ggTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCttt                                                                             <  2:1320561/126‑1 (MQ=255)
                                                                                             gTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATcc                                                    >  2:652572/1‑150 (MQ=255)
                                                                                             gTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATcc                                                    >  1:405069/1‑150 (MQ=255)
                                                                                              tAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCg                                                   <  2:533468/150‑1 (MQ=255)
                                                                                              tAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCg                                                   <  1:1274702/150‑1 (MQ=255)
                                                                                                aGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGa                                                 >  1:391633/1‑150 (MQ=255)
                                                                                                 ggAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAg                                                <  1:1173555/150‑1 (MQ=255)
                                                                                                       aCTGACGAAAGT‑CAGTTCAATTTACTAAAGTCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGAGACGATCCGGAGCCAGGc                                          <  2:203702/150‑1 (MQ=255)
                                                                                                       aCTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGc                                          <  2:1508106/150‑1 (MQ=255)
                                                                                                        cTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCa                                         >  2:543818/1‑150 (MQ=255)
                                                                                                         tGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACAAGACGAGAACTTAAGCCTGAGGCTGAGTTACAACGTCTTAGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAg                                        >  1:832282/1‑150 (MQ=255)
                                                                                                          gACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGt                                       <  2:1207364/150‑1 (MQ=255)
                                                                                                                     aGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgc                             <  2:762395/150‑1 (MQ=255)
                                                                                                                      gTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                            <  1:160680/150‑1 (MQ=255)
                                                                                                                          aaTTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCAc                        <  1:461857/150‑1 (MQ=255)
                                                                                                                           aTTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACg                       >  2:1050801/1‑150 (MQ=255)
                                                                                                                           aTTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACg                       >  1:1319762/1‑150 (MQ=255)
                                                                                                                             ttACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCt                     >  1:602587/1‑150 (MQ=255)
                                                                                                                               aCTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGt                   <  1:307080/150‑1 (MQ=255)
                                                                                                                                  aaaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTc                <  2:1478441/150‑1 (MQ=255)
                                                                                                                                   aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCa               <  2:1319762/150‑1 (MQ=255)
                                                                                                                                      gCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGt            <  1:461153/150‑1 (MQ=255)
                                                                                                                                       cAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACACTGTTTCACGtt           >  2:1675905/1‑150 (MQ=255)
                                                                                                                                       cAAAAAAAACACCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGGCGATAAGTTCAGCAAGATGTTTCACGtt           >  1:957098/1‑150 (MQ=255)
                                                                                                                                         aaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt         >  2:571258/1‑150 (MQ=255)
                                                                                                                                         aaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt         <  2:1277134/150‑1 (MQ=255)
                                                                                                                                         aaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt         <  2:1189644/150‑1 (MQ=255)
                                                                                                                                         aaaaaaaCCACGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCGGCGGCTGAGTTACTACGTCTTTGATACCAATAAATACGATCTCTACGCGCCGATCCGGAGCAAGGCAGTATATCAGTGCAGAACGCTGTTTCACTTTGt         >  2:121332/1‑150 (MQ=255)
                                                                                                                                          aaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTc        >  1:132484/1‑150 (MQ=255)
                                                                                                                                           aaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTca       <  2:1639580/150‑1 (MQ=255)
                                                                                                                                            aaaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcac      >  2:161246/1‑150 (MQ=255)
                                                                                                                                              aaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAg    >  1:924506/1‑150 (MQ=255)
                                                                                                                                               aCCCCGCTGCAGGGGGGTTTTTCTACCAGACGTGAACTTAAGCCTGCGGCTGAGTTACTACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAgg   <  1:146866/150‑1 (MQ=255)
                                                                                                                                                ccccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATTCCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAggg  >  1:42928/1‑149 (MQ=255)
                                                                                                                                                ccccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGgt  <  2:194846/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
AAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGT‑CAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGT  >  NC_000913/1018834‑1019126

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