Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,886,757 T→C 6.9% 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 6.9% 367.1 / 12.8 116intergenic (+35/‑71)yidC/mnmEmembrane protein insertase/tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase
Reads supporting (aligned to +/- strand):  ref base T (53/55);  new base C (4/4);  total (57/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

Only 100 of 122 total aligned reads displayed.

GGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGC  >  NC_000913/3886609‑3886905
                                                                                                                                                    |                                                                                                                                                     
ggTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTAAGt                                                                                                                                                       <  3:454754/149‑1 (MQ=255)
 gTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGtt                                                                                                                                                      >  4:816407/1‑149 (MQ=255)
  tCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTg                                                                                                                                                     >  4:259648/1‑149 (MQ=255)
     ggTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGAc                                                                                                                                                   <  2:516066/148‑1 (MQ=255)
        gCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGcc                                                                                                                                              >  4:522194/1‑149 (MQ=255)
        gCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGc                                                                                                                                               >  1:4399/1‑148 (MQ=255)
         cTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAACTGACCGccc                                                                                                                                             <  4:5846/149‑2 (MQ=255)
           gTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCtt                                                                                                                                            >  2:339948/1‑148 (MQ=255)
           gTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCtt                                                                                                                                            >  2:567746/1‑148 (MQ=255)
            tACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCtttt                                                                                                                                          <  4:760722/149‑1 (MQ=255)
            tACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCtttt                                                                                                                                          <  4:186188/149‑1 (MQ=255)
             aCTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCtttt                                                                                                                                          >  2:442775/1‑148 (MQ=255)
             aCTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCtttt                                                                                                                                          >  2:416274/1‑148 (MQ=255)
              cTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGAAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCtttttt                                                                                                                                        <  3:204279/149‑1 (MQ=255)
               tataTCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCtttttt                                                                                                                                        <  2:818154/148‑1 (MQ=255)
                atatCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGATCGCCTTttttc                                                                                                                                       >  2:481568/1‑148 (MQ=255)
                atatCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTttttc                                                                                                                                       <  1:558581/148‑1 (MQ=255)
                 tatCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGTGCTGTCAGTTGACCGCCTTttttct                                                                                                                                      <  2:639736/148‑1 (MQ=255)
                 tatCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTttttct                                                                                                                                      <  1:194379/148‑1 (MQ=255)
                    cGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTttttcttttc                                                                                                                                  <  3:453429/149‑1 (MQ=255)
                      tCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTACATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCg                                                                                                                                 >  1:118960/1‑148 (MQ=255)
                       cAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGt                                                                                                                                <  2:24824/148‑1 (MQ=255)
                       cAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGt                                                                                                                                <  2:576091/148‑1 (MQ=255)
                        aGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTa                                                                                                                               <  1:128412/148‑1 (MQ=255)
                           aaCCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTttttct                                                                                                                                      <  2:183503/138‑1 (MQ=255)
                           aaCCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAggg                                                                                                                            <  2:4399/148‑1 (MQ=255)
                             ccTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCg                                                                                                                          <  2:948301/148‑1 (MQ=255)
                              cTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCgg                                                                                                                         <  1:128023/148‑1 (MQ=255)
                               tGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACGGCCTTTTTTCTTTTCGTAGGGCGGAt                                                                                                                       >  3:277332/1‑149 (MQ=255)
                                  tAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTa                                                                                                                               <  1:766962/138‑1 (MQ=255)
                                  tAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTa                                                                                                                               >  2:766962/1‑138 (MQ=255)
                                  tAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAg                                                                                                                    <  3:522194/149‑1 (MQ=255)
                                      cATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCAcc                                                                                                                >  4:121171/1‑149 (MQ=255)
                                       attattCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATaa                                                                                                                     <  1:374751/143‑1 (MQ=255)
                                       attattCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACcg                                                                                                               >  3:491722/1‑149 (MQ=255)
                                            tCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCa                                                                                                           >  1:678335/1‑148 (MQ=255)
                                             cagcagcagCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCAt                                                                                                          >  1:688915/1‑148 (MQ=255)
                                               gcagcagcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCg                                                                                                       >  4:118826/1‑149 (MQ=255)
                                                  gcagcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCca                                                                                                    <  3:129337/149‑1 (MQ=255)
                                                    agcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCcac                                                                                                   <  1:13364/148‑1 (MQ=255)
                                                    agcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCAAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCcac                                                                                                   <  2:414240/148‑1 (MQ=255)
                                                      cTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCcacac                                                                                                 >  1:344257/1‑148 (MQ=255)
                                                       tGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATaa                                                                                                                     >  4:111514/1‑127 (MQ=255)
                                                       tGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATaa                                                                                                                     <  3:111514/127‑1 (MQ=255)
                                                            tACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCa                                                                                           <  1:618545/148‑1 (MQ=255)
                                                             aCCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAAc                                                                                         <  3:21284/149‑1 (MQ=255)
                                                              ccGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                        <  3:724514/149‑1 (MQ=255)
                                                                    tCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                        >  2:851325/1‑143 (MQ=255)
                                                                    tCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                        <  1:851325/143‑1 (MQ=255)
                                                                      tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                        <  4:253356/141‑1 (MQ=255)
                                                                      tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                        >  3:253356/1‑141 (MQ=255)
                                                                      tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAcatc                                                                                >  3:288506/1‑149 (MQ=255)
                                                                      tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAcat                                                                                 <  1:527255/148‑1 (MQ=255)
                                                                        gAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAcatcat                                                                              <  3:121171/149‑1 (MQ=255)
                                                                          aaaaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATg                                                                             >  2:70017/1‑148 (MQ=255)
                                                                          aaaaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGa                                                                            <  4:166402/149‑1 (MQ=255)
                                                                            aaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAg                                                                           <  2:351111/148‑1 (MQ=255)
                                                                               gTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGAt                                                                       >  3:592776/1‑149 (MQ=255)
                                                                                  gCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAAt                                                                    <  4:326240/149‑1 (MQ=255)
                                                                                      gCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGac                                                                 <  1:202189/148‑1 (MQ=255)
                                                                                      gCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGac                                                                 >  1:82215/1‑148 (MQ=255)
                                                                                              gcgAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTa                                                        <  4:329003/149‑1 (MQ=255)
                                                                                               cgAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTa                                                        <  2:261527/148‑1 (MQ=255)
                                                                                               cgAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGAACAGCGCATCCGCCACACAAAGCGACAGGAACATCATGAGCGATAATGACACTATCGTa                                                        >  2:273872/1‑148 (MQ=255)
                                                                                                      aaaaaaTCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAgg                                                 >  1:750673/1‑148 (MQ=255)
                                                                                                      aaaaaaTCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAgg                                                 >  1:925712/1‑148 (MQ=255)
                                                                                                      aaaaaaGCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAgg                                                 >  2:531036/1‑148 (MQ=255)
                                                                                                         aaaTCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCa                                              >  2:192690/1‑148 (MQ=255)
                                                                                                            tCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGc                                           >  2:401882/1‑148 (MQ=255)
                                                                                                               tGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTcc                                       >  4:287060/1‑149 (MQ=255)
                                                                                                                gATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTcc                                       >  2:835370/1‑148 (MQ=255)
                                                                                                                      gTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACg                                 <  1:914071/148‑1 (MQ=255)
                                                                                                                       tGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTg                               >  4:752110/1‑149 (MQ=255)
                                                                                                                         aGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTggc                             <  4:59395/149‑1 (MQ=255)
                                                                                                                               cgttAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGtt                       <  3:634976/146‑1 (MQ=255)
                                                                                                                               cGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGt                        <  1:436481/148‑1 (MQ=255)
                                                                                                                               cGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGt                        <  1:413776/148‑1 (MQ=255)
                                                                                                                               cGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGt                        >  1:850744/1‑148 (MQ=255)
                                                                                                                                gCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTg                      <  3:259648/149‑1 (MQ=255)
                                                                                                                                 cTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTg                      <  1:70017/148‑1 (MQ=255)
                                                                                                                                  tAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGc                    >  3:306116/1‑149 (MQ=255)
                                                                                                                                   aaaaTAAGCGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATGCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACATGGCGGCGTTGGc                    >  2:736976/1‑148 (MQ=255)
                                                                                                                                      aTAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATcc                >  3:724458/1‑149 (MQ=255)
                                                                                                                                       tAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATcc                >  2:224758/1‑148 (MQ=255)
                                                                                                                                        aaaggCGGTCAACTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCg                                      <  1:156895/122‑1 (MQ=255)
                                                                                                                                        aaaggCGGTCAACTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCg                                      >  2:156895/4‑125 (MQ=255)
                                                                                                                                        aaGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCt               >  1:810668/1‑148 (MQ=255)
                                                                                                                                         aaggCGGTCAACTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGCTAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTg              >  1:850054/3‑148 (MQ=255)
                                                                                                                                         aGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTgc             <  3:939154/149‑1 (MQ=255)
                                                                                                                                          gggCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAg                                                                           <  3:248188/86‑1 (MQ=255)
                                                                                                                                          gggCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAg                                                                           >  4:248188/1‑86 (MQ=255)
                                                                                                                                          gggCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTgcg            >  3:143734/1‑149 (MQ=255)
                                                                                                                                          gggCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTgcg            >  4:813539/1‑149 (MQ=255)
                                                                                                                                          aggCGGTCAACTGACCGCCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTg              >  2:659690/2‑148 (MQ=255)
                                                                                                                                           ggCGGTCAACTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTAt                                                            <  1:543629/100‑1 (MQ=255)
                                                                                                                                              ggTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCAtc        <  4:307471/149‑1 (MQ=255)
                                                                                                                                                 cAGTTTACCGCCTTTTTTCGTTTCGTAGGGCGGATAAGCACCGCGCCTCCGCCACACAAAGCAACATGGACATCCTGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCAtctc      >  2:376492/1‑148 (MQ=255)
                                                                                                                                                  aGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGTGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCg    >  4:911905/1‑149 (MQ=255)
                                                                                                                                                   gTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCgg   >  4:929631/1‑149 (MQ=255)
                                                                                                                                                    ttGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGc  >  3:615918/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGC  >  NC_000913/3886609‑3886905

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