Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,258,247 A→G F41L (TTC→CTC)  psuG ← pseudouridine 5'‑phosphate glycosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,258,2470AG96.5% 373.0 / 2.9 114F41L (TTC→CTC) psuGpseudouridine 5'‑phosphate glycosidase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (49/61);  total (51/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

Only 100 of 114 total aligned reads displayed.

TATGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGAT  >  NC_000913/2258100‑2258388
                                                                                                                                                   |                                                                                                                                             
tATGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGTATTTTGTGGGAGCg                                                                                                                                             <  1:928626/150‑1 (MQ=255)
tATGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCgg                                                                                                                                            >  2:982647/1‑151 (MQ=255)
   gCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTTTTTACGAATAGTTTTTTCAACTTCAATTGCGGTCTGGGCATTTTTTTGGAGCGGCAt                                                                                                                                         >  1:829629/1‑151 (MQ=255)
   gCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCAt                                                                                                                                         <  1:704572/151‑1 (MQ=255)
     ccTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATcc                                                                                                                                       <  2:1596295/151‑1 (MQ=255)
       ttCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAACGCCCCGTTTACGAATAGTTTTTTCAACTTCAATTGGGGGCTGGGCATTTTTTGGGAGCGGCATccc                                                                                                                                      >  1:221829/1‑150 (MQ=255)
        tCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATccc                                                                                                                                      <  1:353771/149‑1 (MQ=255)
         cACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGt                                                                                                                                    <  1:810863/150‑1 (MQ=255)
              cccAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAACGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGaa                                                                                                                               >  1:257849/1‑150 (MQ=255)
               ccAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGaaa                                                                                                                              >  1:999080/1‑150 (MQ=255)
                cAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAaata                                                                                                                            <  1:160828/151‑1 (MQ=255)
                cAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAaat                                                                                                                             >  1:112568/1‑150 (MQ=255)
                       tCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTc                                                                                                                      <  2:69876/150‑1 (MQ=255)
                        cAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCg                                                                                                                     >  1:1463351/1‑150 (MQ=255)
                          aTTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGAtt                                                                                                                  <  1:1138449/151‑1 (MQ=255)
                          aTTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGAtt                                                                                                                  <  1:285911/151‑1 (MQ=255)
                            ttcttcTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGGCTGGGCATTTTTTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTcc                                                                                                                >  1:423014/1‑151 (MQ=255)
                            ttcttcTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGGACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTcc                                                                                                                >  2:969851/1‑151 (MQ=255)
                             tcttctTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTc                                                                                                                 <  1:52040/149‑1 (MQ=255)
                             tcttctCTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACCGCGCCCTTGTTACGAATAGGTTCTTCTACATTAATTGAGGGCTGGGCATTTTTTGGGAGGGGGCTCCCGTGAGACATAATGGTCGATTcc                                                                                                                >  1:453479/1‑150 (MQ=255)
                                tcttTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAgc                                                                                                             <  1:637776/150‑1 (MQ=255)
                                tcttTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAg                                                                                                              <  2:11145/149‑1 (MQ=255)
                                   ttGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAgcgc                                                                                                           <  1:1021350/149‑1 (MQ=255)
                                   ttGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAgcg                                                                                                            <  1:610900/148‑1 (MQ=255)
                                     gCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCcac                                                                                                        <  2:1226880/150‑1 (MQ=255)
                                      cTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTTCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCCGCGCcac                                                                                                        >  2:139277/1‑149 (MQ=255)
                                       ttAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACaa                                                                                                      >  2:385071/1‑150 (MQ=255)
                                        tAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAAc                                                                                                     <  1:1529061/150‑1 (MQ=255)
                                         aaCCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCg                                                                                                   <  2:238940/151‑1 (MQ=255)
                                         aaCCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCg                                                                                                   <  2:1508499/151‑1 (MQ=255)
                                          aCCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTAAAGCGCCCTGTTTACGAAAAGTTTCTTCAAATTCAATTGCGGTTTGGGCATTTTTTGGGAGCGGCATCCCGTGGGAAATAATGGGGGATTTCCCCGCCCCAACCgg                                                                                                  >  1:1121454/1‑151 (MQ=255)
                                           ccccACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGAGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgt                                                                                                 <  2:139210/151‑1 (MQ=255)
                                            cccACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgt                                                                                                 >  1:376699/1‑150 (MQ=255)
                                             ccACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgttt                                                                                               <  2:843945/151‑1 (MQ=255)
                                              cACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAAACGgtttt                                                                                              >  1:989926/1‑151 (MQ=255)
                                                  ttCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTAGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgttttttg                                                                                           <  2:1586002/150‑1 (MQ=255)
                                                   tcatcaCGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACGACCGgttttttg                                                                                           <  1:454790/149‑1 (MQ=255)
                                                   tcatcaCGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTAAGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGGATTTTGTGGGAGCGGCATCCCGTGGGGAATAATGGTCGATTCCAGCGCCACAACCGgttttttgtt                                                                                         >  2:798227/1‑151 (MQ=255)
                                                    catcaCGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCCGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACCACCGggtttttgttt                                                                                        >  2:1098479/1‑151 (MQ=255)
                                                      tcaCGCCGCCAATAATGGCGATCGTGGCAGGGACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgttttttgttttt                                                                                      >  1:1222289/1‑151 (MQ=255)
                                                       caCGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgttttttgtttttt                                                                                     <  1:38592/151‑1 (MQ=255)
                                                       caCGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgttttttgttttt                                                                                      <  2:4390/150‑1 (MQ=255)
                                                           ccgccAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTaaa                                                                                  <  1:243283/150‑1 (MQ=255)
                                                            cgccAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGc                                                                                <  1:1346063/151‑1 (MQ=255)
                                                             gccAATAATGGCGATCGTGGCAGGGACAGCGCCCTGTTTTCGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCCCCACCGGTTTTTTGTTTTTTAAAGCg                                                                               >  2:1458923/1‑151 (MQ=255)
                                                             gccAATAATGGCGATCGTGGCAGGGACAGCGCCCTGTTTACGAATAGTTTTTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGGGCGGCATCCCGTGAGAAATAATTGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTaa                                                                                   >  1:599639/1‑147 (MQ=255)
                                                               cAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTc                                                                             <  2:1512641/151‑1 (MQ=255)
                                                                   aatGGCGATCGTGGCAGGTACCGCGCCCTGTTTACGAATAGTTTTTTCCACTTCAATTGCGGTCTGGGCCTTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTGGATTCCAGCGCCCCAACCGGTTTTTTGTTTTTTAAAGCGGCCTGc                                                                         >  1:381176/1‑151 (MQ=255)
                                                                   aatGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGc                                                                         <  2:79036/151‑1 (MQ=255)
                                                                            cGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCGCAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCg                                                                  <  2:1577452/149‑1 (MQ=255)
                                                                            cGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCgg                                                                 <  2:291841/150‑1 (MQ=255)
                                                                            cGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCgg                                                                 <  1:1642746/150‑1 (MQ=255)
                                                                              tGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCCCAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCgggg                                                               >  1:889245/1‑150 (MQ=255)
                                                                              tGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCgggg                                                               >  1:1641407/1‑150 (MQ=255)
                                                                               ggCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGaa                                                             <  2:863261/151‑1 (MQ=255)
                                                                               ggCAGGTACAGCGCCCTGCTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGa                                                              <  1:897276/150‑1 (MQ=255)
                                                                                gCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGaaa                                                            <  2:1648252/151‑1 (MQ=255)
                                                                                gCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGaa                                                             <  2:557534/150‑1 (MQ=255)
                                                                                     tACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAAttt                                                         <  1:1099675/149‑1 (MQ=255)
                                                                                     tACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAAttt                                                         >  2:1099761/1‑149 (MQ=255)
                                                                                     tACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTg                                                        <  2:990003/150‑1 (MQ=255)
                                                                                      aCAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGta                                                      <  1:1060272/151‑1 (MQ=255)
                                                                                      aCAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGt                                                       >  2:719966/1‑150 (MQ=255)
                                                                                       cAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGta                                                      >  1:650105/1‑150 (MQ=255)
                                                                                         gcgcCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGtaata                                                   <  2:529221/151‑1 (MQ=255)
                                                                                         gcgcCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGtaata                                                   <  1:1516133/151‑1 (MQ=255)
                                                                                          cgcCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGtaata                                                   >  2:1649473/1‑150 (MQ=255)
                                                                                           gcCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCGGCACTTCCGGGGAAATTTGTaataat                                                 >  2:353604/1‑151 (MQ=255)
                                                                                             ccTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTc                                               <  2:919127/151‑1 (MQ=255)
                                                                                             ccTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTc                                               >  2:419411/1‑151 (MQ=255)
                                                                                                gTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAg                                             <  2:1560551/150‑1 (MQ=255)
                                                                                                gTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAg                                             <  1:1182566/150‑1 (MQ=255)
                                                                                                 tttACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAgg                                            <  1:206058/150‑1 (MQ=255)
                                                                                                  ttACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAggg                                           <  2:1001267/150‑1 (MQ=255)
                                                                                                    aCGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGaa                                        <  2:614294/151‑1 (MQ=255)
                                                                                                     cGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGaaa                                       <  1:268863/151‑1 (MQ=255)
                                                                                                      gAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAAt                                      >  2:1236881/1‑151 (MQ=255)
                                                                                                       aaTAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAAtt                                     <  2:1463450/151‑1 (MQ=255)
                                                                                                        aTAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAAttt                                    <  1:559373/151‑1 (MQ=255)
                                                                                                         tAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAAttt                                    <  2:1525525/150‑1 (MQ=255)
                                                                                                          aGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAAttttt                                  >  1:240120/1‑151 (MQ=255)
                                                                                                               cttcaacttcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAAtt                              >  1:1251789/1‑150 (MQ=255)
                                                                                                                 tcaacttcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAAtt                              <  2:437885/148‑1 (MQ=255)
                                                                                                                 tcaacttcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTc                             >  2:1543250/1‑149 (MQ=255)
                                                                                                                 tcaacttcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTc                             >  1:1170475/1‑149 (MQ=255)
                                                                                                                 tcaacttcaaTTGCGGTCTGGGCATCTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTc                             >  2:1084915/1‑149 (MQ=255)
                                                                                                                   aacttcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTGATAATTCAGGGGAAATTTTTAATTCAGAc                         >  2:1390002/1‑151 (MQ=255)
                                                                                                                   aacttcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGAc                         >  1:237760/1‑151 (MQ=255)
                                                                                                                   aacttcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGAc                         <  2:1679497/151‑1 (MQ=255)
                                                                                                                     cttcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACAt                       >  1:106961/1‑151 (MQ=255)
                                                                                                                       tcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTc                     >  2:1326969/1‑151 (MQ=255)
                                                                                                                        caaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCt                    >  1:780170/1‑151 (MQ=255)
                                                                                                                          aTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGc                  <  1:1112382/151‑1 (MQ=255)
                                                                                                                               ggTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTtct             >  1:1570512/1‑151 (MQ=255)
                                                                                                                               ggTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTtct             >  2:1716220/1‑151 (MQ=255)
                                                                                                                                 tCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTcc           >  2:910190/1‑151 (MQ=255)
                                                                                                                                    gggCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACt        <  1:982572/151‑1 (MQ=255)
                                                                                                                                      gCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTaa      <  1:1573121/151‑1 (MQ=255)
                                                                                                                                        aTTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACg    <  1:39090/151‑1 (MQ=255)
                                                                                                                                           ttGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGat  <  1:1277368/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TATGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGAT  >  NC_000913/2258100‑2258388

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