Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,212,088 (C)8→9 intergenic (‑85/+615) iraM ← / ← ycgX RpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,212,0801.C96.7% 211.8 / 0.1 61intergenic (‑77/+623)iraM/ycgXRpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (26/33);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

CCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTAT  >  NC_000913/1211936‑1212224
                                                                                                                                                 |                                                                                                                                                  
ccccAAATGGCTGAAAAAGATATACCACATGAAGGTTGAATTACCTTGACAATTACGAGCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTAAATTA‑‑ccc‑ccc                                                                                                                                            >  2:308102/1‑151 (MQ=255)
  ccATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccc‑ccccc                                                                                                                                          >  2:539665/1‑151 (MQ=255)
    atatGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccc‑cccccct                                                                                                                                        <  1:819189/151‑3 (MQ=255)
      atGGCTGAAAACGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccc‑cccccctc                                                                                                                                       >  2:759920/1‑147 (MQ=255)
        ggCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccc‑cccccctcaa                                                                                                                                     <  2:1192461/150‑6 (MQ=255)
        ggCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccc‑cccccctcaa                                                                                                                                     <  2:69728/150‑6 (MQ=255)
        ggCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccc‑cccccctcaa                                                                                                                                     <  1:1072443/150‑6 (MQ=255)
          cTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCaaaa                                                                                                                                    >  1:809159/1‑150 (MQ=255)
            gAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGGGTATATTAC‑CCC‑CCCCCTCaaaaa                                                                                                                                   <  2:825349/149‑1 (MQ=255)
            gAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCaaaaaaa                                                                                                                                 <  2:1018777/151‑1 (MQ=255)
             aaaaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAaaat                                                                                                                                <  2:902512/151‑1 (MQ=255)
             aaaaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATAGGAGTATATTAC‑CCC‑CCCCCTCaaa                                                                                                                                     <  2:976436/146‑1 (MQ=255)
              aaaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAaaata                                                                                                                               <  2:526/151‑1 (MQ=255)
               aaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAaaataa                                                                                                                              <  2:1191707/151‑1 (MQ=255)
                    tataCCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAaatta                                                                                                                          >  1:841032/1‑150 (MQ=255)
                      taCCACATGTAGGTTGAATCACCGTCTCATTTACTCTCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATTTTTATTTCTTGATTTTGCACCTGATAATATTAC‑CCC‑CCCCCTCAAAAAAATAaattaatt                                                                                                                       >  1:210694/1‑151 (MQ=255)
                        ccACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAattaatta                                                                                                                      <  1:691147/150‑1 (MQ=255)
                        ccACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAattaatta                                                                                                                      <  1:923129/150‑1 (MQ=255)
                         cacaTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAAttaattaa                                                                                                                     <  2:660152/150‑1 (MQ=255)
                          acaTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTaaa                                                                                                                    >  1:439623/1‑150 (MQ=255)
                              gTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAatga                                                                                                                <  2:1214603/150‑1 (MQ=255)
                               tAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAatgatg                                                                                                              >  1:1098913/1‑151 (MQ=255)
                               tAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAatgatg                                                                                                              >  2:695781/1‑151 (MQ=255)
                                aGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATgg                                                                                                             >  2:553261/1‑151 (MQ=255)
                                     gAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtat                                                                                                        >  1:575071/1‑151 (MQ=255)
                                      aaTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtata                                                                                                       <  2:892569/151‑1 (MQ=255)
                                         tACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtatata                                                                                                     <  1:1181962/150‑1 (MQ=255)
                                                  aaTTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATaaaataaaat                                                                                            <  1:251294/150‑1 (MQ=255)
                                                  aaTTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAtt                                                                                           >  2:27660/1‑151 (MQ=255)
                                                  aaTTACTATCAACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCC‑CCCCCTCAAAACAATAAATTAATTAAAATGATGGCGAATataa                                                                                                    >  2:972513/1‑143 (MQ=255)
                                                   aTTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATAtaaaataaaact                                                                                           <  2:907506/150‑3 (MQ=255)
                                                        tATCCACTTCATTTTTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTAATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTaa                                                                                        <  2:680514/148‑1 (MQ=255)
                                                           ccACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAg                                                                                      <  2:1142041/147‑1 (MQ=255)
                                                             aCTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAgg                                                                                 >  1:575866/1‑150 (MQ=255)
                                                             aCTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGa                                                                                <  2:1124015/151‑1 (MQ=255)
                                                                       ttATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATgg                                                                      >  1:639846/1‑151 (MQ=255)
                                                                          tGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGAt                                                                    <  2:575841/150‑1 (MQ=255)
                                                                          tGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGAt                                                                    <  2:130029/150‑1 (MQ=255)
                                                                                aTCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa                                                              >  1:212372/1‑150 (MQ=255)
                                                                                aTCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAAc                                                             <  1:85386/151‑1 (MQ=255)
                                                                                     aCGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaaaa                                                         <  2:1098929/150‑1 (MQ=255)
                                                                                         tATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATCTAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTATACAAAAt                                                        >  2:676078/1‑147 (MQ=255)
                                                                                           tttAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAtt                                                  >  2:1102378/1‑151 (MQ=255)
                                                                                              aaTATGGTTCATTAGGATATTTATTTCTTGATTTTGCATATAAGTATATTA‑‑CCCACCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAGACAAGGAATCTCAATGGATCTTAAACAAAATGAGATTTtg                                               <  2:759677/151‑1 (MQ=255)
                                                                                                  tGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGa                                            >  1:1119109/1‑150 (MQ=255)
                                                                                                   ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGaa                                           <  1:977915/150‑1 (MQ=255)
                                                                                                    gTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAg                                         <  1:612952/151‑1 (MQ=255)
                                                                                                    gTTCATTAGGATGTTTAATTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATCGATGTTAAACAAAATGAGATTTTGTGAAAg                                         <  2:744728/151‑1 (MQ=255)
                                                                                                     ttCATTAGGAAGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGa                                            <  2:1023494/147‑1 (MQ=255)
                                                                                                      tCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGc                                        >  2:117394/1‑150 (MQ=255)
                                                                                                      tCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGc                                        <  1:1124439/150‑1 (MQ=255)
                                                                                                       cATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCaa                                      <  1:1144290/151‑1 (MQ=255)
                                                                                                       cATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCa                                       <  1:924346/150‑1 (MQ=255)
                                                                                                        aTTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAAt                                     >  1:625391/1‑151 (MQ=255)
                                                                                                        aTTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAAt                                     <  1:394492/151‑1 (MQ=255)
                                                                                                         ttAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAAt                                     >  1:1075661/1‑150 (MQ=255)
                                                                                                            ggATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAatt                                >  1:885366/1‑151 (MQ=255)
                                                                                                                       tCTTGATTTTGCATATGAGTATATTACCCCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACtt                        <  1:506066/150‑1 (MQ=255)
                                                                                                                        cTTGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGt                     >  1:514430/1‑151 (MQ=255)
                                                                                                                          tGATTTTGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCg                      <  1:211391/148‑1 (MQ=255)
                                                                                                                              tttGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGtttt                  >  2:640857/1‑148 (MQ=255)
                                                                                                                              tttGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAg                <  2:841026/150‑1 (MQ=255)
                                                                                                                               ttGCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGtttt                  >  1:115701/1‑147 (MQ=255)
                                                                                                                                 gCATATGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGAtt             <  2:743658/150‑1 (MQ=255)
                                                                                                                                      tGAGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGttt        >  2:901359/1‑150 (MQ=255)
                                                                                                                                        aGTATATTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGc     >  2:711170/1‑151 (MQ=255)
                                                                                                                                           atatTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCTAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCta   >  1:928988/1‑150 (MQ=255)
                                                                                                                                           atatTAC‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCtat  >  2:169590/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCC‑CCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTAT  >  NC_000913/1211936‑1212224

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