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.

GGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGC  >  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)
        gCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑cc                                                                                                                                              >  4:522194/1‑149 (MQ=255)
        gCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑c                                                                                                                                               >  1:4399/1‑148 (MQ=255)
         cTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAACTGACCG‑ccc                                                                                                                                             <  4:5846/149‑2 (MQ=255)
           gTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtt                                                                                                                                            >  2:339948/1‑148 (MQ=255)
           gTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtt                                                                                                                                            >  2:567746/1‑148 (MQ=255)
            tACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtttt                                                                                                                                          <  4:760722/149‑1 (MQ=255)
            tACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtttt                                                                                                                                          <  4:186188/149‑1 (MQ=255)
             aCTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtttt                                                                                                                                          >  2:442775/1‑148 (MQ=255)
             aCTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtttt                                                                                                                                          >  2:416274/1‑148 (MQ=255)
              cTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGAAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtttttt                                                                                                                                        <  3:204279/149‑1 (MQ=255)
               tataTCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtttttt                                                                                                                                        <  2:818154/148‑1 (MQ=255)
                atatCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGATCG‑CCTTttttc                                                                                                                                       >  2:481568/1‑148 (MQ=255)
                atatCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTttttc                                                                                                                                       <  1:558581/148‑1 (MQ=255)
                 tatCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGTGCTGTCAGTTGACCG‑CCTTttttct                                                                                                                                      <  2:639736/148‑1 (MQ=255)
                 tatCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTttttct                                                                                                                                      <  1:194379/148‑1 (MQ=255)
                    cGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTttttcttttc                                                                                                                                  <  3:453429/149‑1 (MQ=255)
                      tCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTACATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCg                                                                                                                                 >  1:118960/1‑148 (MQ=255)
                       cAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGt                                                                                                                                <  2:24824/148‑1 (MQ=255)
                       cAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGt                                                                                                                                <  2:576091/148‑1 (MQ=255)
                        aGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTa                                                                                                                               <  1:128412/148‑1 (MQ=255)
                           aaCCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTttttct                                                                                                                                      <  2:183503/138‑1 (MQ=255)
                           aaCCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAggg                                                                                                                            <  2:4399/148‑1 (MQ=255)
                             ccTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCg                                                                                                                          <  2:948301/148‑1 (MQ=255)
                              cTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCgg                                                                                                                         <  1:128023/148‑1 (MQ=255)
                               tGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACGG‑CCTTTTTTCTTTTCGTAGGGCGGAt                                                                                                                       >  3:277332/1‑149 (MQ=255)
                                  tAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTa                                                                                                                               <  1:766962/138‑1 (MQ=255)
                                  tAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTa                                                                                                                               >  2:766962/1‑138 (MQ=255)
                                  tAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAg                                                                                                                    <  3:522194/149‑1 (MQ=255)
                                      cATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCAcc                                                                                                                >  4:121171/1‑149 (MQ=255)
                                       attattCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATaa                                                                                                                     <  1:374751/143‑1 (MQ=255)
                                       attattCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACcg                                                                                                               >  3:491722/1‑149 (MQ=255)
                                            tCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCa                                                                                                           >  1:678335/1‑148 (MQ=255)
                                             cagcagcagCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCAt                                                                                                          >  1:688915/1‑148 (MQ=255)
                                               gcagcagcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCg                                                                                                       >  4:118826/1‑149 (MQ=255)
                                                  gcagcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCca                                                                                                    <  3:129337/149‑1 (MQ=255)
                                                    agcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCcac                                                                                                   <  1:13364/148‑1 (MQ=255)
                                                    agcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCAAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCcac                                                                                                   <  2:414240/148‑1 (MQ=255)
                                                      cTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCcacac                                                                                                 >  1:344257/1‑148 (MQ=255)
                                                       tGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATaa                                                                                                                     >  4:111514/1‑127 (MQ=255)
                                                       tGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATaa                                                                                                                     <  3:111514/127‑1 (MQ=255)
                                                            tACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCa                                                                                           <  1:618545/148‑1 (MQ=255)
                                                             aCCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAAc                                                                                         <  3:21284/149‑1 (MQ=255)
                                                              ccGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                        <  3:724514/149‑1 (MQ=255)
                                                                    tCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                        >  2:851325/1‑143 (MQ=255)
                                                                    tCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                        <  1:851325/143‑1 (MQ=255)
                                                                      tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                        <  4:253356/141‑1 (MQ=255)
                                                                      tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                        >  3:253356/1‑141 (MQ=255)
                                                                      tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAcatc                                                                                >  3:288506/1‑149 (MQ=255)
                                                                      tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAcat                                                                                 <  1:527255/148‑1 (MQ=255)
                                                                        gAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAcatcat                                                                              <  3:121171/149‑1 (MQ=255)
                                                                          aaaaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATg                                                                             >  2:70017/1‑148 (MQ=255)
                                                                          aaaaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGa                                                                            <  4:166402/149‑1 (MQ=255)
                                                                            aaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAg                                                                           <  2:351111/148‑1 (MQ=255)
                                                                               gTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGAt                                                                       >  3:592776/1‑149 (MQ=255)
                                                                                  gCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAAt                                                                    <  4:326240/149‑1 (MQ=255)
                                                                                      gCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGac                                                                 <  1:202189/148‑1 (MQ=255)
                                                                                      gCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGac                                                                 >  1:82215/1‑148 (MQ=255)
                                                                                              gcgAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTa                                                        <  4:329003/149‑1 (MQ=255)
                                                                                               cgAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTa                                                        <  2:261527/148‑1 (MQ=255)
                                                                                               cgAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGAACAGCGCATCCGCCACACAAAGCGACAGGAACATCATGAGCGATAATGACACTATCGTa                                                        >  2:273872/1‑148 (MQ=255)
                                                                                                      aaaaaaTCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAgg                                                 >  1:750673/1‑148 (MQ=255)
                                                                                                      aaaaaaTCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAgg                                                 >  1:925712/1‑148 (MQ=255)
                                                                                                      aaaaaaGCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAgg                                                 >  2:531036/1‑148 (MQ=255)
                                                                                                         aaaTCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCa                                              >  2:192690/1‑148 (MQ=255)
                                                                                                            tCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGc                                           >  2:401882/1‑148 (MQ=255)
                                                                                                               tGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTcc                                       >  4:287060/1‑149 (MQ=255)
                                                                                                                gATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTcc                                       >  2:835370/1‑148 (MQ=255)
                                                                                                                      gTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACg                                 <  1:914071/148‑1 (MQ=255)
                                                                                                                       tGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTg                               >  4:752110/1‑149 (MQ=255)
                                                                                                                         aGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTggc                             <  4:59395/149‑1 (MQ=255)
                                                                                                                               cgttAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGtt                       <  3:634976/146‑1 (MQ=255)
                                                                                                                               cGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGt                        <  1:436481/148‑1 (MQ=255)
                                                                                                                               cGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGt                        <  1:413776/148‑1 (MQ=255)
                                                                                                                               cGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGt                        >  1:850744/1‑148 (MQ=255)
                                                                                                                                gCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTg                      <  3:259648/149‑1 (MQ=255)
                                                                                                                                 cTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTg                      <  1:70017/148‑1 (MQ=255)
                                                                                                                                  tAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGc                    >  3:306116/1‑149 (MQ=255)
                                                                                                                                   aaaaTAAGCGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATGCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACATGGCGGCGTTGGc                    >  2:736976/1‑148 (MQ=255)
                                                                                                                                      aTAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATcc                >  3:724458/1‑149 (MQ=255)
                                                                                                                                       tAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATcc                >  2:224758/1‑148 (MQ=255)
                                                                                                                                        aaaggCGGTCAACTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCg                                      <  1:156895/122‑1 (MQ=255)
                                                                                                                                        aaaggCGGTCAACTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCg                                      >  2:156895/4‑125 (MQ=255)
                                                                                                                                        aaGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCt               >  1:810668/1‑148 (MQ=255)
                                                                                                                                         aaggCGGTCAACTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGCTAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTg              >  1:850054/3‑148 (MQ=255)
                                                                                                                                         aGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTgc             <  3:939154/149‑1 (MQ=255)
                                                                                                                                          gggCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAg                                                                           <  3:248188/86‑1 (MQ=255)
                                                                                                                                          gggCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAg                                                                           >  4:248188/1‑86 (MQ=255)
                                                                                                                                          gggCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTgcg            >  3:143734/1‑149 (MQ=255)
                                                                                                                                          gggCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTgcg            >  4:813539/1‑149 (MQ=255)
                                                                                                                                          aggCGGTCAACTGACCGCCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTg              >  2:659690/2‑148 (MQ=255)
                                                                                                                                           ggCGGTCAACTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTAt                                                            <  1:543629/100‑1 (MQ=255)
                                                                                                                                              ggTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCAtc        <  4:307471/149‑1 (MQ=255)
                                                                                                                                                 cAGTTTACCG‑CCTTTTTTCGTTTCGTAGGGCGGATAAGCACCGCGCCTCCGCCACACAAAGCAACATGGACATCCTGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCAtctc      >  2:376492/1‑148 (MQ=255)
                                                                                                                                                  aGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGTGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCg    >  4:911905/1‑149 (MQ=255)
                                                                                                                                                   gTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCgg   >  4:929631/1‑149 (MQ=255)
                                                                                                                                                    ttGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGc  >  3:615918/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGC  >  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: ‑