Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,601,000 T→A 6.1% intergenic (‑217/+28) 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,601,0000TA 6.1% 384.0 / 9.6 127intergenic (‑217/+28)rpoH/ftsXRNA polymerase, sigma 32 (sigma H) factor/putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (59/60);  new base A (6/2);  total (65/62)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 129 total aligned reads displayed.

CAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCA  >  NC_000913/3600851‑3601148
                                                                                                                                                       |                                                                                                                                                    
cAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCAtt                                                                                                                                                      >  1:1069437/1‑150 (MQ=255)
 aaCGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACACACAAAATGTGACATAGAGATG‑A‑ATTACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTg                                                                                                                                                     <  1:297685/150‑1 (MQ=255)
   cGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCa                                                                                                                                                   <  2:1787621/150‑1 (MQ=255)
      ggTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAggg                                                                                                                                                <  2:1157476/150‑1 (MQ=255)
       gTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAAACCATTGCAGGGa                                                                                                                                               >  1:1735452/1‑150 (MQ=255)
         tATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGaaa                                                                                                                                             >  2:375832/1‑150 (MQ=255)
         tATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGaaa                                                                                                                                             <  2:919472/150‑1 (MQ=255)
         tATCATTCACTGTTTTATCAGACCGTGATTTTATCCAAAAGTTCTATGCAAGCTTGTGAATTAATTTCGCACAAAATGTGTCATAGAGATG‑A‑ATTACAGGGAATAGACAACGGGGTCGCTTTCCCTGCAATGGGTTCCGTAGCAGGGaaa                                                                                                                                             <  1:346020/150‑1 (MQ=255)
          aTCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAag                                                                                                                                            <  2:1320772/150‑1 (MQ=255)
          aTCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCAATGGGTTCCGTAGCAGGGAAag                                                                                                                                            <  2:1666355/150‑1 (MQ=255)
            cATTCACTGTTTTATCAGACAGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGAAATAGACATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAagag                                                                                                                                          >  1:497458/1‑150 (MQ=255)
             agtCACTGTTTAATCAGACCCTGATTTTATCCACAAGTTCAATGCAAGCTTCTTAATAAATAACGCACAAAATTTGAAATAGAGATG‑ATTATA‑AGGGAAGCGAAAACGGGGTCTCTTACCATTCTACGGCACCCATTGCAGGGAAagaga                                                                                                                                         >  2:1223171/3‑149 (MQ=255)
             aTTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCGATTGCAGGGAAAGAGt                                                                                                                                         >  2:1191178/1‑150 (MQ=255)
               tCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtat                                                                                                                                       <  1:1431211/150‑1 (MQ=255)
                cACTGTTTTATCAGACCATGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAA‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtat                                                                                                                                       >  2:1482076/1‑150 (MQ=255)
                cACTGTTATATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGAGG‑A‑AATACCGGGAAGAGACAACGGTGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtata                                                                                                                                      >  1:286383/1‑150 (MQ=255)
                 aCTGTTTTATCAGACCGTGATTTTATCCACAAGATCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGATACGGAACACATTGCAGGGAAAGAGTATaa                                                                                                                                     >  2:1191433/1‑150 (MQ=255)
                  cTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAac                                                                                                                                    >  1:723189/1‑150 (MQ=255)
                   tGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAaca                                                                                                                                   <  1:1514814/150‑1 (MQ=255)
                    gTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAacac                                                                                                                                  >  2:450736/1‑150 (MQ=255)
                    gTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAacac                                                                                                                                  >  1:1194600/1‑150 (MQ=255)
                        tATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCtt                                                                                                                              >  2:1273546/1‑150 (MQ=255)
                            aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                                          >  2:701504/1‑150 (MQ=255)
                               ccGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGACCCCATTGCAGGGAAAGAGTAAAACCCGCTTTtattat                                                                                                                       >  2:1362004/1‑150 (MQ=255)
                               ccGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTtattat                                                                                                                       <  2:1351581/150‑1 (MQ=255)
                                cGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTattatt                                                                                                                      <  1:973978/150‑1 (MQ=255)
                                 gTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTattatta                                                                                                                     >  1:1295521/1‑149 (MQ=255)
                                 gTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTc                                                                                                                     <  1:1217484/150‑1 (MQ=255)
                                    aTTTTATACACAAGTTCCATGCAAGCTTGTGAATAAATTACGCACAAAAGGTGACATAGAGATG‑A‑AATACCGTGAAGAGACAACGGGGGCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                                  <  2:1268567/150‑1 (MQ=255)
                                       ttATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGt                                                                                                               <  2:839179/150‑1 (MQ=255)
                                        tATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTa                                                                                                              >  2:1026339/1‑150 (MQ=255)
                                          tCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTaaa                                                                                                            >  2:1287858/1‑150 (MQ=255)
                                           ccACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAg                                                                                                           >  1:25132/1‑150 (MQ=255)
                                            cacaAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGt                                                                                                          >  2:863053/1‑150 (MQ=255)
                                              caAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTgg                                                                                                        <  1:750792/150‑1 (MQ=255)
                                                 gTTCAATGCAAGCTTGTGAATAAATTCCGCACAAAATGGGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGt                                                                                                     <  2:734788/150‑1 (MQ=255)
                                                 gTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGt                                                                                                     >  1:1754093/1‑150 (MQ=255)
                                                 gTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGt                                                                                                     >  1:1099294/1‑150 (MQ=255)
                                                      aTGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGTAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtg                                                                                                <  2:137702/150‑1 (MQ=255)
                                                       tGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgt                                                                                               >  1:1770781/1‑150 (MQ=255)
                                                        gCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgtt                                                                                              >  2:786594/1‑150 (MQ=255)
                                                        gCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgtt                                                                                              <  2:944668/150‑1 (MQ=255)
                                                         cAAGATTGTGAATAAATAACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGTCAACGGGGTCTCTTTCCCTGCTACGGATCCCATTGCAGTGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAATGTGGCGTAAAtgttg                                                                                             <  1:1191433/150‑1 (MQ=255)
                                                          aaGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgttgt                                                                                            >  2:802411/1‑150 (MQ=255)
                                                          aaGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgttgt                                                                                            <  1:1183434/150‑1 (MQ=255)
                                                              ttGTGAACAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTATTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTAAGGCGTAAAGTTTCGTAAATGTTGTACCg                                                                                        >  1:773396/1‑150 (MQ=255)
                                                                 tgAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTgg                                                                                     >  2:1619857/1‑150 (MQ=255)
                                                                 tgAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGACACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCATTGAAATAGTAAAACACGCTTTTATTATTCAGGCGGAAAGTGGCGTAAATGTTGTACCGTgg                                                                                     <  2:280389/150‑1 (MQ=255)
                                                                   aaTAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTTCGGAACCCATTGCAGGGAAAGAGTATAACAAGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTTTACCGTGGCa                                                                                   >  1:1623565/1‑150 (MQ=255)
                                                                         ttACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCAc                                                                             <  1:739841/150‑1 (MQ=255)
                                                                           aCGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGc                                                                           <  1:600995/150‑1 (MQ=255)
                                                                           aCGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGc                                                                           <  1:199339/150‑1 (MQ=255)
                                                                            cGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCt                                                                          >  1:12694/1‑150 (MQ=255)
                                                                             ccacaAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTg                                                                         <  2:463936/149‑1 (MQ=255)
                                                                              cacaAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTTAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGc                                                                        >  1:763208/1‑150 (MQ=255)
                                                                              cacaAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGc                                                                        >  2:1199078/1‑150 (MQ=255)
                                                                               acaAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGcc                                                                       <  2:1308100/150‑1 (MQ=255)
                                                                                   aaTGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCAccc                                                                   <  1:360836/150‑1 (MQ=255)
                                                                                     tgtgACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAg                                                                 <  2:135591/150‑1 (MQ=255)
                                                                                     tgggACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAg                                                                 <  2:1031921/147‑1 (MQ=255)
                                                                                        gACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGAGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCa                                                              >  2:1065384/1‑150 (MQ=255)
                                                                                          cATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAat                                                            <  1:500180/150‑1 (MQ=255)
                                                                                               agaTG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCg                                                       <  2:746082/150‑1 (MQ=255)
                                                                                                 aTG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAg                                                     <  1:1512804/150‑1 (MQ=255)
                                                                                                  tG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGc                                                    >  1:384940/1‑150 (MQ=255)
                                                                                                  tG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGc                                                    >  1:1374838/1‑150 (MQ=255)
                                                                                                  tG‑A‑AATACCGGGAAGAGACAAAGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCCAGCCACGCTGCCACCCAGCCAATCATCCAGc                                                    >  2:136482/1‑150 (MQ=255)
                                                                                                          aCCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACCCGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATCTTGTACCGTGGCCAGCCACGCTGACCCCCAGCCAATCATCGAGCATACca                                              >  2:330032/1‑150 (MQ=255)
                                                                                                          aCCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACca                                              <  1:735987/150‑1 (MQ=255)
                                                                                                           ccGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACcag                                             >  2:1199300/1‑150 (MQ=255)
                                                                                                             gggAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACcagca                                           >  2:1485121/1‑150 (MQ=255)
                                                                                                               gAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCagcagc                                         <  1:1089580/150‑1 (MQ=255)
                                                                                                               gAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGAATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCagcagc                                         >  1:522089/1‑150 (MQ=255)
                                                                                                                 agagACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCaa                                       <  2:1304712/150‑1 (MQ=255)
                                                                                                                  gagaCAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCATGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGATGCCACCCATCCAATAAGAGAGAATACCAGCAGCAAt                                      >  2:93156/1‑150 (MQ=255)
                                                                                                                   agaCAACGGTGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATa                                     <  2:1736762/150‑1 (MQ=255)
                                                                                                                     aCAACTGGGTCTCTTTCCCTGCTACTGAACCCATTGCAGGGAAAGAGTATAACACGCTTTAATAATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGc                                   <  1:294406/150‑1 (MQ=255)
                                                                                                                      cAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCa                                  <  2:689066/150‑1 (MQ=255)
                                                                                                                       aaCGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAg                                 <  2:204497/150‑1 (MQ=255)
                                                                                                                            ggTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGTCGTAAAGTTGCGTAAATGTTGTACCGTGGCAAGCCAAGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCAtt                            <  2:286383/150‑1 (MQ=255)
                                                                                                                            atactctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCAtt                            >  2:394828/4‑150 (MQ=38)
                                                                                                                            acactctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCAtt                            >  2:1569875/4‑150 (MQ=255)
                                                                                                                              tctcTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCa                          <  2:885453/150‑1 (MQ=255)
                                                                                                                                tctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATc                        <  2:513865/150‑1 (MQ=255)
                                                                                                                                tctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATc                        >  1:1332484/1‑150 (MQ=39)
                                                                                                                                 ctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTCCCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCGTACCTGCAGCAATAGCAGGCATTCATCg                       >  1:1132172/1‑150 (MQ=255)
                                                                                                                                 ctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCg                       >  2:1508648/1‑150 (MQ=255)
                                                                                                                                  tttCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGa                      >  1:372465/1‑150 (MQ=255)
                                                                                                                                      ccTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGa                  >  2:862323/1‑150 (MQ=255)
                                                                                                                                       cTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGAt                 >  1:1703161/1‑150 (MQ=255)
                                                                                                                                       cTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGAt                 >  1:1092682/1‑150 (MQ=255)
                                                                                                                                           tACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGACACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGc             >  1:223792/1‑150 (MQ=255)
                                                                                                                                              ggAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCAt          <  2:1770781/150‑1 (MQ=255)
                                                                                                                                              ggAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCAAt          >  1:1100355/1‑150 (MQ=255)
                                                                                                                                               gAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGAGTAAATGTTGTACCGTGGCAAGCCACGCAGCCACACAGCCAATCATCGAGCATACCAGCAGCAAAATCATGCATTCATCCAATGATAAGCCaat         >  1:516182/1‑148 (MQ=255)
                                                                                                                                                aaCCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTg        <  2:367180/150‑1 (MQ=255)
                                                                                                                                                aaCCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGACAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTg        >  2:1227422/1‑150 (MQ=255)
                                                                                                                                                 aCCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGa       <  1:802411/150‑1 (MQ=255)
                                                                                                                                                 aCCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGa       <  2:666854/150‑1 (MQ=255)
                                                                                                                                                      atGCAGGGAAAGAGTATAACACGCTTTTATAATTCAGGCGTAAAGTGGCGTAAATGTTGTAACGTGGAAAGCCACGCTGCCAACCAGCCAATCATCGAGCATACCAGCAGCAATATCAGGCATTCATCGAATGATAAGCCATTGATATCa  <  1:876706/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
CAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCA  >  NC_000913/3600851‑3601148

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: ‑