Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,600,988 T→A 5.3% intergenic (‑205/+40) rpoH ← / ← ftsX RNA polymerase, sigma 32 (sigma H) factor/putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,600,9880TA 5.3% 377.1 / 5.2 126intergenic (‑205/+40)rpoH/ftsXRNA polymerase, sigma 32 (sigma H) factor/putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (60/58);  new base A (6/2);  total (66/60)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

AGAGCTTAAGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAA  >  NC_000913/3600839‑3601135
                                                                                                                                                     |                                                                                                                                                   
    cTTAAGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATTTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACgg                                                                                                                                                 >  1:972977/1‑150 (MQ=255)
      tAAGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGaa                                                                                                                                               >  2:1180794/1‑150 (MQ=255)
      tAAGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGaa                                                                                                                                               >  2:1574260/1‑150 (MQ=255)
       aaGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAAc                                                                                                                                              >  1:308507/1‑150 (MQ=255)
           gCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCAt                                                                                                                                          >  2:1344899/1‑150 (MQ=255)
           gCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCAt                                                                                                                                          >  1:1216945/1‑150 (MQ=255)
             aaCGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTg                                                                                                                                        <  2:486410/150‑1 (MQ=255)
                gAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAg                                                                                                                                     >  1:1188709/1‑150 (MQ=255)
                gAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAg                                                                                                                                     <  2:1531820/150‑1 (MQ=255)
                 aGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAgg                                                                                                                                    <  1:253512/150‑1 (MQ=255)
                 aGGTTATCAATCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAgg                                                                                                                                    <  1:1292002/150‑1 (MQ=255)
                  ggTTATCATTCACTGTTTTATCAGACAGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAATAGACAACGGTGTCTCTTTCCATGCAACGGAACCCATTGCAggg                                                                                                                                   >  1:799489/1‑150 (MQ=255)
                   gTTATCCTTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATACATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGa                                                                                                                                  <  2:252637/150‑1 (MQ=255)
                    ttATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACAAAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGCGaa                                                                                                                                 >  1:387738/1‑150 (MQ=255)
                      aTCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCAATGGGTTCCGTAGCAGGGAAag                                                                                                                               <  2:1450333/150‑1 (MQ=255)
                      aTCATTCACGGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCGGCTACGGAACCCATTGCAGGGAAag                                                                                                                               <  2:58405/150‑1 (MQ=255)
                       ccATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAaga                                                                                                                              <  1:338091/149‑1 (MQ=255)
                         aTTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACACCGGGGTCTCTTTCCCTGCTACGGAACCGAGTGCATGGAATGAGt                                                                                                                            >  2:1187698/1‑150 (MQ=255)
                         aTTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGt                                                                                                                            >  2:1647741/1‑150 (MQ=255)
                         aTTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGt                                                                                                                            <  1:1699964/150‑1 (MQ=255)
                          ttCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCACGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCAATGGGTTCCGTCGCAGGGAAagagac                                                                                                                           <  1:899899/150‑3 (MQ=38)
                              cTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTGGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAac                                                                                                                       <  1:1697164/150‑1 (MQ=255)
                                gTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAacac                                                                                                                     >  2:1008298/1‑150 (MQ=255)
                                  tttATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGc                                                                                                                   <  1:1825026/150‑1 (MQ=255)
                                     aTCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCttt                                                                                                                <  2:730891/150‑1 (MQ=255)
                                       cAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACAAGCTTtta                                                                                                              >  2:1670303/1‑150 (MQ=255)
                                       cAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTATCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTtta                                                                                                              <  1:1575300/150‑1 (MQ=255)
                                        aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                             <  1:643900/150‑1 (MQ=255)
                                        aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                             >  1:1148880/1‑150 (MQ=255)
                                        aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                             <  2:1758717/150‑1 (MQ=255)
                                        aGACCCTGATATTATCCACAAGTTCAATGCAAGCTTGTGAATCAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                             <  2:898852/150‑1 (MQ=255)
                                            cGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTattatt                                                                                                         <  1:836991/150‑1 (MQ=255)
                                               gATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATCACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAg                                                                                                      <  1:636720/150‑1 (MQ=255)
                                                aTTTTATCCACTAGTTCAATGCAAGCTTGTGAATAAATTACGCACATAATGTGACATAGAGATGATATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                     <  2:421947/150‑1 (MQ=255)
                                                aTTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAATTACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGCACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTAATCAgg                                                                                                     >  1:1478331/1‑150 (MQ=255)
                                                aTTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                     >  1:1521084/1‑150 (MQ=255)
                                                      tCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGACCCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTaaa                                                                                               >  2:961116/1‑150 (MQ=255)
                                                       ccACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTaaat                                                                                              >  2:674247/1‑149 (MQ=255)
                                                       ccACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAg                                                                                              >  2:690230/1‑150 (MQ=255)
                                                        cacaAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGt                                                                                             >  2:10686/1‑150 (MQ=255)
                                                        cacaAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTAAGGCGTAAAGt                                                                                             >  1:520121/1‑150 (MQ=255)
                                                         acaAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTg                                                                                            >  2:1761731/1‑150 (MQ=255)
                                                            aGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCg                                                                                         <  2:207024/150‑1 (MQ=255)
                                                             gTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGt                                                                                        >  1:1311820/1‑150 (MQ=255)
                                                                 aaTGCAAGCGTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAt                                                                                    <  1:1457637/150‑1 (MQ=255)
                                                                  aTGCAAGCTTGTTAATAAATTACGCACAATATGTGACATAGAGATTAAATACCGGGAAGAGACAACGGGGTATCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtg                                                                                   <  1:275608/150‑1 (MQ=255)
                                                                  aTGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtt                                                                                   >  2:1464939/1‑149 (MQ=255)
                                                                    gCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgtt                                                                                 >  2:955592/1‑150 (MQ=255)
                                                                    gCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgtt                                                                                 <  2:62580/150‑1 (MQ=255)
                                                                    gCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgtt                                                                                 <  1:1344899/150‑1 (MQ=255)
                                                                           tgtgAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGt                                                                          <  2:1216945/150‑1 (MQ=255)
                                                                           tgtgAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCATGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTAATATTATTCAGACGTAAAGTGACGTAAATGTTGTACCGt                                                                          >  1:1324425/1‑150 (MQ=255)
                                                                             tgAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTgg                                                                        >  1:424116/1‑150 (MQ=255)
                                                                             tgAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTgg                                                                        >  1:1736157/1‑150 (MQ=255)
                                                                             tgAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGAAAATGTTGTACCGTgg                                                                        >  2:593108/1‑150 (MQ=255)
                                                                             tgAAGAAATTACGCACAAAATGTGACATAGAGATGAACTACCGGCCTGAGACAACGGGGGCTCTTTCCCTGCTACGGAACCCCTTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTgg                                                                        <  2:372383/150‑1 (MQ=255)
                                                                              gAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGc                                                                       >  1:870863/1‑150 (MQ=255)
                                                                                aTAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCaa                                                                     <  2:1006949/150‑1 (MQ=255)
                                                                                 tAAAATACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAg                                                                    <  2:1386485/150‑1 (MQ=255)
                                                                                   aaTTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGcc                                                                  >  2:820757/1‑150 (MQ=255)
                                                                                      tACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACg                                                               <  1:1083392/150‑1 (MQ=255)
                                                                                          cacaAAATGTGACATAGAGAAGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGTAAGAGTATAACACGCTGTTATTATTCTGGCGTAAAGAGGCGTAAATGTTGTCCCTTTGCAAGCCCAGCTGc                                                           >  2:350153/1‑150 (MQ=255)
                                                                                             aaaaTGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCAc                                                        >  1:620172/1‑150 (MQ=255)
                                                                                                 tgtgACATAGAGATGAAATACCGGGAAGAGACAACGGTGTCTCTTTCCCTGCTACGGAACCCTTTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAACTGTTGTACCGTGGCAAGCCACGCTGCCACCCAg                                                    <  1:354559/150‑1 (MQ=255)
                                                                                                 tgtgACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAg                                                    <  2:896702/150‑1 (MQ=255)
                                                                                                  gtgACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGc                                                   <  2:1086418/150‑1 (MQ=255)
                                                                                                   ttaCATAGAGATGAAATACCGTGAAGAGACAACGGGCTCTCATTCCCAGCTACGGACCCCCTAGCAGGGAAAGAGTAGACCACGCTTTTAGTCTTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTCGCAAGCCACGCTGCCACCCAGcc                                                  <  1:279378/148‑1 (MQ=255)
                                                                                                   tgACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGcc                                                  <  2:1449690/150‑1 (MQ=255)
                                                                                                      cATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAat                                               >  1:1300008/1‑150 (MQ=255)
                                                                                                        tAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAatca                                             >  1:1276213/1‑150 (MQ=255)
                                                                                                         agagATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAatcat                                            >  2:24101/1‑150 (MQ=255)
                                                                                                          gagaTGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAatcatc                                           <  1:1003243/150‑1 (MQ=255)
                                                                                                            gaTGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCgg                                         >  2:557571/1‑149 (MQ=255)
                                                                                                             aTGAAATACCGGGAAGCGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGGATAACACGCTTTTATTATTCAGGCGTACAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAg                                        <  1:992945/150‑1 (MQ=255)
                                                                                                                  aTACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATAc                                   <  2:973101/150‑1 (MQ=255)
                                                                                                                   tACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATAcc                                  <  2:943388/150‑1 (MQ=255)
                                                                                                                   gcagggaaaGAGACCCCGTTGTCTCTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACca                                 >  1:1575801/8‑150 (MQ=255)
                                                                                                                    aCCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACca                                 >  2:1757381/1‑150 (MQ=255)
                                                                                                                     ccGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTAAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACcac                                >  2:1401211/1‑149 (MQ=255)
                                                                                                                      cGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATAAcagc                               >  2:746924/1‑150 (MQ=255)
                                                                                                                         gAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCagcagc                            >  2:112401/1‑150 (MQ=255)
                                                                                                                         gAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCagcagc                            <  1:239058/150‑1 (MQ=255)
                                                                                                                         gAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCagcagc                            <  2:586756/150‑1 (MQ=255)
                                                                                                                           agagACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCagcagcca                          >  1:1332306/1‑148 (MQ=255)
                                                                                                                             agaCAACGGGGTCTCTTTACCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATa                        <  2:646566/150‑1 (MQ=255)
                                                                                                                              gaCAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAg                       >  1:480625/1‑150 (MQ=255)
                                                                                                                              gaCAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAg                       >  1:11775/1‑150 (MQ=255)
                                                                                                                                  aCGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAgg                   <  1:1552788/150‑1 (MQ=255)
                                                                                                                                       gTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTc              >  1:1108282/1‑150 (MQ=255)
                                                                                                                                       gTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTc              >  1:1661156/1‑150 (MQ=255)
                                                                                                                                          tctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATc           <  1:198205/150‑1 (MQ=255)
                                                                                                                                           ctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCATTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCAGTCATCg          >  1:959862/1‑150 (MQ=255)
                                                                                                                                           ctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATCGCAGGCATTCATCg          >  1:1625734/1‑150 (MQ=255)
                                                                                                                                            tttCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGa         >  1:1526503/1‑150 (MQ=255)
                                                                                                                                            tttCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGa         >  2:761963/1‑150 (MQ=255)
                                                                                                                                              tCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAAt       <  1:1180794/150‑1 (MQ=255)
                                                                                                                                               cccTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATg      <  1:1761731/150‑1 (MQ=255)
                                                                                                                                               cccTGCTACGGAACCCATTGCAGCGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGAGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAGTCATCGAGAATCCCAGCAGCAATAGCAGGAATTCATCGAATg      >  2:150323/1‑150 (MQ=255)
                                                                                                                                                 cTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGAt    >  1:651408/1‑150 (MQ=255)
                                                                                                                                                  tGCTACGGAACCCATTGCATGGAAAGAGTATAACACGCTTTTATTATTCATGCTTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATa   <  2:1533995/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGAGCTTAAGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAA  >  NC_000913/3600839‑3601135

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