Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,993,896 T→A 5.4% intergenic (+40/‑43) ygeI → / → pbl uncharacterized protein/pseudogene, peptidoglycan‑binding enzyme family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,993,8960TA 5.4% 340.5 / 6.9 130intergenic (+40/‑43)ygeI/pbluncharacterized protein/pseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base T (60/63);  new base A (2/5);  total (62/68)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

Only 100 of 132 total aligned reads displayed.

ATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAACAATGT  >  NC_000913/2993747‑2994044
                                                                                                                                                     |                                                                                                                                                    
aTACAGATCAAATACAAGAATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                                                                                      <  1:786144/150‑1 (MQ=255)
  aCAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                                                                                                                                                    <  1:1230999/150‑1 (MQ=255)
    aGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa                                                                                                                                                  >  1:180658/1‑150 (MQ=255)
    aGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTa                                                                                                                                                  <  2:135123/150‑1 (MQ=255)
       tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtct                                                                                                                                               <  1:1275231/150‑1 (MQ=255)
       tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtct                                                                                                                                               <  2:207176/150‑1 (MQ=255)
        cAAATACAAGCATTTATAGTTTCCACATGGATGGCGCCTTTTCAAAATGATATGTATTAAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATAATTTTATAATAGTAtcta                                                                                                                                              >  2:1062186/1‑149 (MQ=255)
         aaaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCg                                                                                                                                             >  2:303642/1‑150 (MQ=255)
          aaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                                                                            <  2:1489207/150‑1 (MQ=255)
          aaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                                                                            <  2:621850/150‑1 (MQ=255)
          aaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                                                                            <  2:1458776/150‑1 (MQ=255)
           aTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGtt                                                                                                                                           <  1:1409491/150‑1 (MQ=255)
            tACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTc                                                                                                                                          <  1:907178/150‑1 (MQ=255)
             aCAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttat                                                                                                                                         <  1:265517/150‑3 (MQ=255)
              aaaGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatg                                                                                                                                        <  2:1551162/149‑4 (MQ=255)
                aGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                                                                                      >  2:241374/1‑150 (MQ=255)
                aGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                                                                                      <  2:1198144/150‑1 (MQ=255)
                 tcATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                                                                                     <  1:792947/149‑1 (MQ=255)
                 tcATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                                                                                     <  2:1392193/149‑1 (MQ=255)
                      tataGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattca                                                                                                                                >  2:128540/1‑150 (MQ=255)
                         aGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATc                                                                                                                             <  1:248821/150‑1 (MQ=255)
                         aGTTTCAACATGGATGGAGCCTTTTCAAAATGATATGTATTCAGAAGATAATACTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATAAATTTATAATAGTATCTCGTTCCTTTCATTAAAtg                                                                                                                             >  1:93625/1‑149 (MQ=255)
                          gTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                                                                            >  1:1495949/1‑150 (MQ=255)
                          gTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                                                                            >  2:818717/1‑150 (MQ=255)
                              cAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcac                                                                                                                        >  1:54172/1‑150 (MQ=255)
                               aaCATGGATGTAGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCATACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcaca                                                                                                                       <  2:46567/150‑1 (MQ=255)
                               aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcaca                                                                                                                       <  2:222770/150‑1 (MQ=255)
                               aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTAAAATAGTATATCGTTCCTTTCATTAAATCCTccca                                                                                                                       >  2:1487026/1‑150 (MQ=255)
                                    ggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGa                                                                                                                  <  2:297328/150‑1 (MQ=255)
                                     gatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaa                                                                                                                 >  2:133333/1‑150 (MQ=255)
                                     gatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaa                                                                                                                 >  1:1377007/1‑150 (MQ=255)
                                     gatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGATCCTTTCATTCAATCCTCACACATGaa                                                                                                                 <  2:1638354/150‑1 (MQ=255)
                                       tggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaaa                                                                                                               <  1:783067/150‑1 (MQ=255)
                                            ccTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                                                                                     >  2:931960/1‑123 (MQ=255)
                                            ccTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAatat                                                                                                          >  1:1551471/1‑150 (MQ=255)
                                             cTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATg                                                                                                         <  2:51233/150‑1 (MQ=255)
                                                  cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTc                                                                                                    >  2:620485/1‑150 (MQ=255)
                                                  cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTc                                                                                                    >  1:477462/1‑150 (MQ=255)
                                                  cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTc                                                                                                    >  2:1641019/1‑150 (MQ=255)
                                                      aTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCa                                                                                                <  1:1001795/150‑1 (MQ=255)
                                                         atatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGAACAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaaaa                                                                                             >  1:212936/1‑150 (MQ=255)
                                                         atatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGAAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaaaa                                                                                             <  2:1120924/150‑1 (MQ=255)
                                                            tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACTCATGAAAAAATATGCGCTCACCAAAAGtt                                                                                          >  1:1129968/1‑150 (MQ=255)
                                                            tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTAAATCCTCACACATGAAAAAATATGCGCTCAACAAAAGtt                                                                                          >  2:542074/1‑150 (MQ=255)
                                                               aTTCAGAAGATAATCCTATCTCACCTTACAATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGATCACCAAAAGTTaaa                                                                                       >  2:1323018/1‑150 (MQ=255)
                                                                ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTTGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAAt                                                                                      <  2:844358/150‑1 (MQ=255)
                                                                ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAAt                                                                                      <  1:363598/150‑1 (MQ=255)
                                                                ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAAt                                                                                      >  1:1036251/1‑150 (MQ=255)
                                                                 tCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCGTCACACATGAAAAAATATGAGCTAACCAAAATTTAAAtt                                                                                     >  2:1526249/1‑150 (MQ=255)
                                                                 tCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATTCGCTCACCAAAAGTTAAAtt                                                                                     >  2:316217/1‑150 (MQ=255)
                                                                 tCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAAtt                                                                                     <  1:956677/150‑1 (MQ=255)
                                                                  cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                                                                                    >  2:58179/1‑150 (MQ=255)
                                                                      agaTAATCCTATCTCACCTTACTATAATATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTaa                                                                                <  2:910686/150‑1 (MQ=255)
                                                                         tAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTGACCa                                                                             <  1:562599/150‑1 (MQ=255)
                                                                          aaTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCAt                                                                            >  2:1640433/1‑150 (MQ=255)
                                                                           aTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTAACAAAAAGTTAAATACTTAACCAta                                                                           >  2:583657/1‑149 (MQ=255)
                                                                                    cACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGttt                                                                  >  2:151181/1‑150 (MQ=255)
                                                                                      ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCCTAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCCCCAAAAGTTAAATTCTTAACCATCTTCACGtattt                                                                >  2:340127/1‑149 (MQ=255)
                                                                                       cTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAACATCTTCACGTTTTGt                                                               >  1:824415/1‑150 (MQ=255)
                                                                                         tACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtat                                                             <  2:1036251/150‑1 (MQ=255)
                                                                                           cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACGATCTTCACGTTTTGTatat                                                           <  2:916898/150‑1 (MQ=255)
                                                                                           cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTatat                                                           >  1:1215859/1‑150 (MQ=255)
                                                                                           cTATAAAATTGAGTGGTAGAAGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCTTTCCTTTAATTCAATCCTCACACATGAAAAAATAAGCGCTCACCAAAAGTTAAATTCTTAAACATCTTCACGTTTTCtataa                                                           >  1:1186019/1‑149 (MQ=255)
                                                                                            tataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatt                                                          <  2:1397722/150‑1 (MQ=255)
                                                                                            tataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATCCGCTCACAAAAAGTTAAATTCTTAACCATCTTCACTTTTTGTAtatt                                                          >  2:248921/1‑150 (MQ=255)
                                                                                             ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtattt                                                         >  2:1719785/1‑150 (MQ=255)
                                                                                                aaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatttatt                                                      <  1:1561384/150‑1 (MQ=255)
                                                                                                 aaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATttatta                                                     <  1:1641019/150‑1 (MQ=255)
                                                                                                 aaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTGCACGTTTTGTATATttatta                                                     <  1:58179/150‑1 (MQ=255)
                                                                                                  aTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTATATATttattaa                                                    >  1:395237/1‑149 (MQ=255)
                                                                                                   ttGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACa                                                   <  2:477462/150‑1 (MQ=255)
                                                                                                    tGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACaa                                                  >  1:889872/1‑150 (MQ=255)
                                                                                                     gAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACaaa                                                 >  2:626669/1‑150 (MQ=255)
                                                                                                       ttGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGGTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACaaaaa                                               <  1:613493/149‑1 (MQ=255)
                                                                                                         ggTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATg                                             <  2:1009904/150‑1 (MQ=255)
                                                                                                         ggTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATg                                             <  2:391452/150‑1 (MQ=255)
                                                                                                          gTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGt                                            >  1:1630059/1‑150 (MQ=255)
                                                                                                           tAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTc                                           <  1:102831/150‑1 (MQ=255)
                                                                                                            aGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCt                                          >  1:1568227/1‑150 (MQ=255)
                                                                                                               ctaTTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTc                                       <  2:1186019/148‑1 (MQ=255)
                                                                                                               cGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTc                                       >  2:1720508/1‑150 (MQ=255)
                                                                                                                 aTTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAAAATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGc                                     >  2:1431584/1‑150 (MQ=255)
                                                                                                                  ttGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGcc                                    <  1:133333/150‑1 (MQ=255)
                                                                                                                    gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAg                                  >  2:6919/1‑150 (MQ=255)
                                                                                                                    gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAg                                  <  2:532714/150‑1 (MQ=255)
                                                                                                                     aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACa                                                   <  1:661356/132‑1 (MQ=255)
                                                                                                                     aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACa                                                   >  2:661356/1‑132 (MQ=255)
                                                                                                                        aaataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAAc                              >  1:1652786/1‑150 (MQ=255)
                                                                                                                         aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAAct                             >  2:448299/1‑150 (MQ=255)
                                                                                                                           taaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAActct                           <  1:178470/150‑1 (MQ=255)
                                                                                                                              atCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGc                        <  1:1171457/150‑1 (MQ=255)
                                                                                                                                aggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGctc                      >  1:594833/4‑150 (MQ=255)
                                                                                                                                 ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTTTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGctct                     >  1:142065/3‑150 (MQ=255)
                                                                                                                                     cGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGATTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAAtt                 >  2:82937/1‑149 (MQ=255)
                                                                                                                                         aTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGc             <  1:1681346/150‑1 (MQ=255)
                                                                                                                                          tACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCa            >  2:742130/1‑150 (MQ=255)
                                                                                                                                               ttATAATAGTATCTCGTTCCTTTCATTCTATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGaat       <  1:248921/150‑2 (MQ=255)
                                                                                                                                                tataATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAACa      >  2:1073489/1‑150 (MQ=255)
                                                                                                                                                tataATAGTATCTCGTTCCTTTCATTCAATCATCACACATTAAAAAATATGATCTAACAAAAATTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTATTTCTCCTGCAACTCATGATCTAATGAAGCATGaaaa      >  1:136970/1‑148 (MQ=255)
                                                                                                                                                    atatTATCTCTTTCCTTTCATTAAATCCTCACACATGAAAAAATATGCGCTCAAAAAAAGTTAAATTCTTTACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGACAGCAACTCTTGCTCTAATGAAGCAGGAACAATGt  <  1:1062186/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAACAATGT  >  NC_000913/2993747‑2994044

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: