Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,993,896 T→A 9.4% intergenic (+40/‑43) ygeI → / → pbl protein YgeI/transglycosylase domain‑containing protein Pbl

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,993,8960TA 9.4% 487.6 / 32.9 178intergenic (+40/‑43)ygeI/pblprotein YgeI/transglycosylase domain‑containing protein Pbl
Reads supporting (aligned to +/- strand):  ref base T (88/73);  new base A (10/7);  total (98/80)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

Only 100 of 184 total aligned reads displayed.

ATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAACAAT  >  NC_000913/2993747‑2994042
                                                                                                                                                     |                                                                                                                                                  
 tACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataata                                                                                                                                                   >  1:616200/1‑150 (MQ=255)
  aCAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                                                                                                                                                  <  1:182776/150‑1 (MQ=255)
    aGATCAAATACAAGCATTTAAAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAATATAAACCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa                                                                                                                                                <  2:1754621/150‑1 (MQ=255)
       tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtct                                                                                                                                             >  2:235394/1‑150 (MQ=255)
          aaTACAATCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAATATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                                                                          <  2:581091/150‑1 (MQ=255)
          aaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCATAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                                                                          <  1:1717047/150‑1 (MQ=255)
          aaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                                                                          >  2:212714/1‑150 (MQ=255)
          aaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                                                                          <  1:540842/150‑1 (MQ=255)
            tACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTc                                                                                                                                        <  1:714465/150‑1 (MQ=255)
            tACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGtta                                                                                                                                        <  1:481664/150‑2 (MQ=255)
              cAACCATTTATAGTCTCAACAGGGATGGCGCCTTATCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTTACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatg                                                                                                                                      <  1:322898/150‑4 (MQ=255)
                  cATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttca                                                                                                                                  <  1:1084054/150‑1 (MQ=255)
                      tataGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                                                                                   <  2:1618514/145‑1 (MQ=255)
                         aGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATc                                                                                                                           <  2:616200/150‑1 (MQ=255)
                          gTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                                                                          >  1:1353731/1‑150 (MQ=255)
                          gTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                                                                          >  1:971506/1‑150 (MQ=255)
                          gTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                                                                          <  1:1023187/150‑1 (MQ=255)
                          atctCAACATTGATGGAGCCTTTTCAAAATGATATGTATTCAGAATATAATCCTATCTCACCTTACTATAAAATTGAGTTGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                                                                          <  2:1420215/147‑1 (MQ=255)
                              cAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcac                                                                                                                      >  1:351962/1‑150 (MQ=255)
                              cAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcac                                                                                                                      <  2:446643/150‑1 (MQ=255)
                               aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcaca                                                                                                                     <  1:842609/150‑1 (MQ=255)
                               aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATAATTTTATAATAGTATCTCGTTCCTTTCATTCAATCATcaca                                                                                                                     >  1:1807573/1‑150 (MQ=255)
                                aCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcacac                                                                                                                    >  1:814343/1‑150 (MQ=255)
                                          cgcCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAat                                                                                                          >  1:1283979/1‑150 (MQ=255)
                                            ccTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAat                                                                                                          >  2:459989/1‑148 (MQ=255)
                                                 tCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCt                                                                                                   <  1:235394/150‑1 (MQ=255)
                                                 tCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCt                                                                                                   >  1:286009/1‑150 (MQ=255)
                                                  cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTc                                                                                                  >  2:1504890/1‑150 (MQ=255)
                                                  cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTc                                                                                                  >  1:1588726/1‑150 (MQ=255)
                                                         atatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaaaa                                                                                           <  2:1060765/150‑1 (MQ=255)
                                                            tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGtt                                                                                        >  1:447802/1‑150 (MQ=255)
                                                             gTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTa                                                                                       >  1:83606/1‑150 (MQ=255)
                                                               aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                                                                          >  1:506650/1‑113 (MQ=255)
                                                                ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAAt                                                                                    >  2:370145/1‑150 (MQ=255)
                                                                  gagaagaTAATCCTATATGACCTTACTATAAAATTGAGTGGTAGAAGATTGAACAAATAAATAATAAAGAGATAAATTTATAAAAATATATAGTTAAATTCATTCAATCCTAACAAAAGAAAAAATAAGAGATAAACAAAACTTaaaata                                                                                  >  2:501882/2‑147 (MQ=255)
                                                                  cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                                                                                  >  2:1883490/1‑150 (MQ=255)
                                                                  cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                                                                                  >  1:449801/1‑150 (MQ=255)
                                                                   agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCt                                                                                 >  1:1972086/1‑150 (MQ=255)
                                                                     aagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTa                                                                               <  2:1360326/150‑1 (MQ=255)
                                                                       taTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                                                                             <  2:351962/149‑1 (MQ=255)
                                                                       gaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                                                                             <  2:252289/150‑1 (MQ=255)
                                                                       gaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                                                                             >  2:777245/1‑150 (MQ=255)
                                                                       gaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                                                                             <  1:1009010/150‑1 (MQ=255)
                                                                           atacTATCTAAAATTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATc                                                                         >  2:649040/4‑150 (MQ=255)
                                                                           aTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATc                                                                         >  2:694669/1‑150 (MQ=255)
                                                                              cTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTc                                                                      >  2:1231767/1‑150 (MQ=255)
                                                                                 tctcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg                                                                   <  1:1512085/150‑1 (MQ=255)
                                                                                 tctcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg                                                                   <  1:1504890/150‑1 (MQ=255)
                                                                                   tcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt                                                                 >  1:1708510/1‑150 (MQ=255)
                                                                                    cACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTGACGttt                                                                <  2:1353731/150‑1 (MQ=255)
                                                                                    cACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGttt                                                                >  2:1265595/1‑150 (MQ=255)
                                                                                      ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg                                                              >  2:1947917/1‑150 (MQ=255)
                                                                                      ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg                                                              >  1:1749687/1‑150 (MQ=255)
                                                                                      ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg                                                              >  1:1000476/1‑150 (MQ=255)
                                                                                      ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg                                                              <  2:841405/150‑1 (MQ=255)
                                                                                         tACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtat                                                           >  2:886283/1‑150 (MQ=255)
                                                                                           cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATATTAGAATGTCGTTCATTTCATTCAATCATCACACTTGAAAAAATATGCGCACACCAAAAGTTAAATTCTAAACAAGGTTCACGttttaaatcg                                                         >  2:1287778/1‑144 (MQ=255)
                                                                                           atataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTatat                                                         <  1:1947917/149‑1 (MQ=255)
                                                                                            tataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatt                                                        <  1:1412064/150‑1 (MQ=255)
                                                                                                aaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatttatt                                                    >  1:1332091/1‑150 (MQ=255)
                                                                                                  aTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTAc                                                  >  1:1579438/1‑150 (MQ=255)
                                                                                                         ggTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGAGCTCAGCAAAAGTTAAATTCTTATCAATCTTCACGTTTTGTATATTTATTAAAAAAATg                                           >  1:1530935/1‑150 (MQ=255)
                                                                                                           tAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg                                                                   <  2:186268/124‑1 (MQ=255)
                                                                                                           tAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg                                                                   >  1:186268/1‑124 (MQ=255)
                                                                                                             gACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCtt                                       >  1:1469841/1‑150 (MQ=255)
                                                                                                             gACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCtt                                       >  1:1767138/1‑150 (MQ=255)
                                                                                                               cGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTc                                     >  2:1837490/1‑150 (MQ=255)
                                                                                                                    gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAg                                >  1:467935/1‑150 (MQ=255)
                                                                                                                     aCCAAATAACTCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCCTCTTCACGGTTTGTATATTTATTACAAAAATGTCTTTCGCCAGc                               <  1:1396454/150‑1 (MQ=255)
                                                                                                                      ccAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCa                              <  1:41048/150‑1 (MQ=255)
                                                                                                                      ccAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCa                              <  2:832627/150‑1 (MQ=255)
                                                                                                                           taaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAActct                         <  1:1883490/150‑1 (MQ=255)
                                                                                                                            aaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCtt                        >  2:117651/1‑150 (MQ=255)
                                                                                                                              atCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTgg                      >  1:592180/1‑149 (MQ=255)
                                                                                                                              atCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGc                      >  1:1325026/1‑150 (MQ=255)
                                                                                                                              atCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGc                      >  1:1666804/1‑150 (MQ=255)
                                                                                                                              atCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTATTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGc                      >  1:1902070/1‑150 (MQ=255)
                                                                                                                                cATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGctc                    >  2:1547253/1‑150 (MQ=255)
                                                                                                                                 ggaaCGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGctct                   <  1:1493186/148‑1 (MQ=255)
                                                                                                                                 ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGctct                   >  2:1116788/3‑150 (MQ=255)
                                                                                                                                 aTAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGctct                   >  1:1382925/1‑150 (MQ=255)
                                                                                                                                  gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTa                  >  2:304855/2‑150 (MQ=255)
                                                                                                                                   aaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTaa                 >  1:1734630/1‑150 (MQ=255)
                                                                                                                                    aCGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAAt                <  1:777245/150‑1 (MQ=255)
                                                                                                                                    aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAAt                >  2:444330/1‑150 (MQ=255)
                                                                                                                                    aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAAt                >  2:1597299/1‑150 (MQ=255)
                                                                                                                                    aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAAt                <  1:1116788/150‑1 (MQ=255)
                                                                                                                                     cGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGATCTAAtt               >  1:1187394/1‑149 (MQ=255)
                                                                                                                                     cGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATg               >  2:1114943/1‑150 (MQ=255)
                                                                                                                                     cGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATg               <  1:1597299/150‑1 (MQ=255)
                                                                                                                                      gagaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGa              >  1:258436/1‑150 (MQ=255)
                                                                                                                                         aTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAGAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGc           <  2:814343/150‑1 (MQ=255)
                                                                                                                                            cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAgg        <  1:633871/150‑1 (MQ=255)
                                                                                                                                             ttttATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGa       <  1:886283/150‑1 (MQ=255)
                                                                                                                                              tttATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa                                                                     <  1:827062/87‑1 (MQ=255)
                                                                                                                                              tttATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGaa      <  1:649040/150‑1 (MQ=255)
                                                                                                                                              tttATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGCATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGaa      <  1:142279/150‑1 (MQ=255)
                                                                                                                                               ttATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAAc     <  2:1972086/150‑1 (MQ=255)
                                                                                                                                                 ataataGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAACaa   <  1:1265595/150‑1 (MQ=255)
                                                                                                                                                  taataGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAACAAt  >  1:1788201/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAACAAT  >  NC_000913/2993747‑2994042

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