Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,886,757 T→C 5.0% intergenic (+35/‑71) yidC → / → mnmE membrane protein insertase/tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,886,7570TC 5.0% 567.4 / 11.0 173intergenic (+35/‑71)yidC/mnmEmembrane protein insertase/tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase
Reads supporting (aligned to +/- strand):  ref base T (82/82);  new base C (4/5);  total (86/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 178 total aligned reads displayed.

GTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCG  >  NC_000913/3886610‑3886903
                                                                                                                                                   |                                                                                                                                                  
gTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGtt                                                                                                                                                   <  4:344384/149‑1 (MQ=255)
gTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAAGCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGtt                                                                                                                                                   <  1:1261509/149‑1 (MQ=255)
 tCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTg                                                                                                                                                  >  1:558841/1‑149 (MQ=255)
  cTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTGCCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGa                                                                                                                                                 <  1:1364309/149‑1 (MQ=255)
      tGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTTTGGAAAATCGTGGCATGCATAGCCGCGAGAAGACAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGc                                                                                                                                             <  3:57994/149‑1 (MQ=255)
        cTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAACTGACCGccc                                                                                                                                           <  1:8422/149‑2 (MQ=255)
         tGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAACTCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAACTGACCGccct                                                                                                                                          <  2:410926/149‑3 (MQ=255)
         tGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAACTGACCGccct                                                                                                                                          <  3:349963/149‑3 (MQ=255)
           tACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCtttt                                                                                                                                        >  1:801854/1‑149 (MQ=255)
            aCTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCttttt                                                                                                                                       <  1:543292/149‑1 (MQ=255)
             cTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCtttttt                                                                                                                                      >  3:601213/1‑149 (MQ=255)
             cTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCtttttt                                                                                                                                      >  3:1184174/1‑149 (MQ=255)
              tataTCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTttttc                                                                                                                                     >  2:98646/1‑149 (MQ=255)
                  tCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTttttctttt                                                                                                                                 <  1:353686/149‑1 (MQ=255)
                   cGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTttttcttttc                                                                                                                                <  3:770694/149‑1 (MQ=255)
                     tCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGt                                                                                                                              >  4:1183310/1‑149 (MQ=255)
                     tCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGt                                                                                                                              >  1:62141/1‑149 (MQ=255)
                      cAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAACCGGGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTa                                                                                                                             <  3:176085/149‑1 (MQ=255)
                       aGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAg                                                                                                                            <  1:752405/149‑1 (MQ=255)
                       aGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTATTCTATTCGTAg                                                                                                                            <  4:1254881/149‑1 (MQ=255)
                             cTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGa                                                                                                                      <  1:957502/149‑1 (MQ=255)
                                gTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATaa                                                                                                                   <  2:600550/149‑1 (MQ=255)
                                 tAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCCGTCAGTTGACCGCCTTTTTTCCTTTCGTAGGGCGGATAAg                                                                                                                  <  1:440410/149‑1 (MQ=255)
                                 tAACCATTATTCAGCAGCAGCTGATTTACCGTGGGCTGGAAAAACGTGGCCTGCATATCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAGATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGGGGATAAg                                                                                                                  >  1:1026825/1‑149 (MQ=255)
                                    ccATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCAc                                                                                                               <  2:669075/149‑1 (MQ=255)
                                     cATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCAcc                                                                                                              >  4:955275/1‑149 (MQ=255)
                                     cATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCAcc                                                                                                              >  3:885289/1‑149 (MQ=255)
                                     cATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCAcc                                                                                                              >  4:1030925/1‑149 (MQ=255)
                                     cATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCAcc                                                                                                              >  4:1288878/1‑149 (MQ=255)
                                       ttattCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACcgc                                                                                                            <  4:1019587/149‑1 (MQ=255)
                                       ttattCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACcgc                                                                                                            >  3:164658/1‑149 (MQ=255)
                                       ttattCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACcgc                                                                                                            <  2:889698/149‑1 (MQ=255)
                                       ttattCAGCAGAAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACcgc                                                                                                            >  1:789392/1‑149 (MQ=255)
                                            cagcagcagCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGGTAAGCACCGCGCATc                                                                                                       <  4:118725/149‑1 (MQ=255)
                                                  cagcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCcac                                                                                                 >  1:606348/1‑149 (MQ=255)
                                                   agcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCcaca                                                                                                >  3:229982/1‑149 (MQ=255)
                                                    gcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGc                                                                                                                         >  1:1369842/1‑123 (MQ=255)
                                                    gcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGc                                                                                                                         <  2:1369842/123‑1 (MQ=255)
                                                    gcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCcacac                                                                                               >  3:1123107/1‑149 (MQ=255)
                                                      tGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACaa                                                                                             >  2:802339/1‑149 (MQ=255)
                                                       gATTTACCGTGGTCTGGAAGAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACaaa                                                                                            <  4:1011497/149‑1 (MQ=255)
                                                       gATTTACCGTGGTCTGGAAAAACGTGGCCTTCATAGTCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTGCGCTAAAAGAAGGGCCGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACaaa                                                                                            <  1:1272228/149‑1 (MQ=255)
                                                       gATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACaaa                                                                                            >  1:1009721/1‑149 (MQ=255)
                                                          ttACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCa                                                                                         >  4:1377833/1‑149 (MQ=255)
                                                           tACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCaa                                                                                        <  3:1425263/149‑1 (MQ=255)
                                                             ccGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCtttt                                                                                                                                        <  2:1185740/99‑1 (MQ=255)
                                                             ccGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                      >  4:1299994/1‑149 (MQ=255)
                                                             ccGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                      <  4:713335/149‑1 (MQ=255)
                                                              cGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAg                                                                                     >  1:1199850/1‑149 (MQ=255)
                                                                  gTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAc                                                                                 <  1:892586/149‑1 (MQ=255)
                                                                  gTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAc                                                                                 <  1:245337/149‑1 (MQ=255)
                                                                   tCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAca                                                                                <  2:1199850/149‑1 (MQ=255)
                                                                     tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATaa                                                                                                                   <  1:1437378/112‑1 (MQ=255)
                                                                     tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATaa                                                                                                                   >  2:1437378/1‑112 (MQ=255)
                                                                     tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAcatc                                                                              <  2:400293/149‑1 (MQ=255)
                                                                     tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAcatc                                                                              >  2:517552/1‑149 (MQ=255)
                                                                     tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAcatc                                                                              >  4:713228/1‑149 (MQ=255)
                                                                      ggAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCgg                                                                                                                       >  2:1031951/1‑107 (MQ=255)
                                                                       gAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAcatcat                                                                            >  1:638492/1‑149 (MQ=255)
                                                                        aaaaaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATcc                                                                                                      <  3:1332282/122‑1 (MQ=255)
                                                                        aaaaaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATg                                                                           <  1:530112/149‑1 (MQ=255)
                                                                          aaaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAg                                                                         >  3:1369182/1‑149 (MQ=255)
                                                                          aaaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAg                                                                         <  3:1116606/149‑1 (MQ=255)
                                                                                 gCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAAt                                                                  <  3:1030925/149‑1 (MQ=255)
                                                                                 gCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAAt                                                                  <  4:1283094/149‑1 (MQ=255)
                                                                                   cTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGa                                                                >  4:221274/1‑149 (MQ=255)
                                                                                   cTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGa                                                                <  1:929555/149‑1 (MQ=255)
                                                                                   cTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGa                                                                <  2:1455020/149‑1 (MQ=255)
                                                                                   cTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGa                                                                >  4:737191/1‑149 (MQ=255)
                                                                                    tGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGac                                                               <  3:887631/149‑1 (MQ=255)
                                                                                    tGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGac                                                               >  2:36463/1‑149 (MQ=255)
                                                                                     gCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGaca                                                              <  3:815484/149‑1 (MQ=255)
                                                                                         aGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTAt                                                          <  4:601213/149‑1 (MQ=255)
                                                                                          gCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATc                                                         <  4:764960/149‑1 (MQ=255)
                                                                                                    gaaAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAgg                                               >  1:517037/1‑149 (MQ=255)
                                                                                                     aaaaaaTCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGc                                              >  3:725336/1‑149 (MQ=255)
                                                                                                       aaaaTCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCTTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTACCGTAGCCCAGGCCa                                            >  1:1071052/1‑149 (MQ=255)
                                                                                                          aTCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                      >  4:500359/1‑104 (MQ=255)
                                                                                                             cTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTc                                      <  3:395439/149‑1 (MQ=255)
                                                                                                              tGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTcc                                     >  1:1128443/1‑149 (MQ=255)
                                                                                                               gATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCg                                    <  2:977517/149‑1 (MQ=255)
                                                                                                                   cGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGAc                                <  3:1008776/149‑1 (MQ=255)
                                                                                                                     gTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGGACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGt                              >  3:1466564/1‑149 (MQ=255)
                                                                                                                     gTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGt                              <  1:998182/149‑1 (MQ=255)
                                                                                                                      tGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTg                             >  1:576162/1‑149 (MQ=255)
                                                                                                                       gAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTgg                            <  2:801854/149‑1 (MQ=255)
                                                                                                                          ttttCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTggcgg                         >  3:233677/1‑149 (MQ=255)
                                                                                                                           tttCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTggcggc                        <  1:1217901/149‑1 (MQ=255)
                                                                                                                              cGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGtt                     <  1:810391/149‑1 (MQ=255)
                                                                                                                               gCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTg                    <  3:1384460/149‑1 (MQ=255)
                                                                                                                               gCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTg                    >  2:872907/1‑149 (MQ=255)
                                                                                                                                 tAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAAt                                                                  >  2:840670/1‑101 (MQ=255)
                                                                                                                                    aaTAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATc               <  1:1476111/149‑1 (MQ=255)
                                                                                                                                      tAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACACCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTg                    <  1:24416/142‑1 (MQ=255)
                                                                                                                                       aaaggCGGTCAACTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTg            >  2:1095587/4‑149 (MQ=255)
                                                                                                                                        aaggCGGTCAACTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTgc           >  1:1331749/3‑149 (MQ=255)
                                                                                                                                          ggCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTgcgc         >  2:52374/1‑149 (MQ=255)
                                                                                                                                              gTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCAtct     >  2:1506869/1‑149 (MQ=255)
                                                                                                                                               tCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCAtctc    >  4:1142290/1‑149 (MQ=255)
                                                                                                                                                 aGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCg  <  4:156416/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCG  >  NC_000913/3886610‑3886903

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