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.C94.3% 293.3 / 11.3 88intergenic (‑77/+623)iraM/ycgXRpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base C (38/45);  total (40/48)
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.29e-01

ACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAAT  >  NC_000913/1211935‑1212227
                                                                                                                                                   |                                                                                                                                                      
aCCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑ccccc                                                                                                                                                  <  2:558685/151‑1 (MQ=255)
   ccATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑cccccccc                                                                                                                                               >  2:863811/1‑151 (MQ=255)
     atatGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑ccccccccct                                                                                                                                             <  1:226920/151‑3 (MQ=255)
      tatGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATC‑A‑‑ccccccccctc                                                                                                                                            >  1:1142877/1‑148 (MQ=255)
       atGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑ccccccccctc                                                                                                                                            <  1:26441/150‑4 (MQ=255)
        tGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑ccccccccctca                                                                                                                                           >  2:585767/1‑146 (MQ=255)
        tGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑ccccccccctc                                                                                                                                            >  2:1112333/1‑146 (MQ=255)
         ggCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑Caaa                                                                                                                                         <  1:98465/151‑1 (MQ=255)
         ggCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑ccccccccctcaa                                                                                                                                          <  2:591550/150‑6 (MQ=255)
              aaaaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑Caaaaaaa                                                                                                                                     <  1:235821/150‑1 (MQ=255)
              aaaaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAaaat                                                                                                                                    <  1:140305/151‑1 (MQ=255)
              aaaaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAaaat                                                                                                                                    <  2:1018569/151‑1 (MQ=255)
              aaaaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAaaat                                                                                                                                    <  2:175791/151‑1 (MQ=255)
               aaaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAaaat                                                                                                                                    >  1:818240/1‑150 (MQ=255)
                 aaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAaaataa                                                                                                                                  >  2:1030900/1‑150 (MQ=255)
                   gATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAaatt                                                                                                                               <  2:313517/151‑1 (MQ=255)
                    atatACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAaatt                                                                                                                               <  1:151235/150‑1 (MQ=255)
                    atatACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAaatt                                                                                                                               >  1:824166/1‑150 (MQ=255)
                        aCCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAaattaatt                                                                                                                           >  1:322195/1‑150 (MQ=255)
                        aCCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAaattaatt                                                                                                                           <  1:163049/150‑1 (MQ=255)
                         ccACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAattaatta                                                                                                                          <  2:221485/150‑1 (MQ=255)
                          cacaTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑ACCCCCCCCCCT‑CAAAAAAATAAAttaattaa                                                                                                                         <  1:1090820/151‑1 (MQ=255)
                          cacaTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAAttaattaa                                                                                                                         >  1:826692/1‑150 (MQ=255)
                          cacaTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTaaa                                                                                                                        <  1:658675/151‑1 (MQ=255)
                            caTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTaaaa                                                                                                                       <  2:972418/150‑1 (MQ=255)
                               gTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAatgat                                                                                                                   <  2:244568/151‑1 (MQ=255)
                               gTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAatga                                                                                                                    >  1:728523/1‑150 (MQ=255)
                                 aGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAatgatg                                                                                                                  <  2:1175343/150‑1 (MQ=255)
                                  gtttaaaTTACCGTTTCAATTACTATCCACTTCATTTTTTATGTCTTATCCCACGGTATTAAAAATGTTCCATTAGGATGTTAATTTTTTAATTTGGAATATGGGAAAATTCC‑‑CCCCCCCCT‑CAAAAAAATAAATTAATTaaaa                                                                                                                       <  1:957246/139‑1 (MQ=255)
                                    ttGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCtt                                                                                                              >  2:787503/1‑151 (MQ=255)
                                      gAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTt‑at                                                                                                           <  2:1194900/151‑1 (MQ=255)
                                      gAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTt‑at                                                                                                           <  2:377513/151‑1 (MQ=255)
                                          tACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑ACCCCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTt‑atat                                                                                                         >  2:454736/1‑150 (MQ=255)
                                          tACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTt‑atata                                                                                                        >  1:764449/1‑150 (MQ=255)
                                           aCCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGTTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATataa                                                                                                       >  1:556781/1‑150 (MQ=255)
                                           aCCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATataaa                                                                                                      <  1:1136952/151‑1 (MQ=255)
                                            ccGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATataaaa                                                                                                     <  2:230682/151‑1 (MQ=255)
                                            ccGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATataaa                                                                                                      <  1:129305/150‑1 (MQ=255)
                                             cGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATataaaat                                                                                                    <  1:897986/151‑1 (MQ=255)
                                                 tCAATTACTATCCACTTCGTTTGTTATGTCTTATCCCACGGTGTTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATataaaataaa                                                                                                 >  1:556468/1‑150 (MQ=255)
                                                 tCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATataaaataaa                                                                                                 >  2:591414/1‑150 (MQ=255)
                                                 tCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATataaaataaa                                                                                                 <  1:875238/150‑1 (MQ=255)
                                                       aCTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTTCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTaa                                                                                           <  1:406718/150‑1 (MQ=255)
                                                          aTCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGc                                                                                        >  1:817329/1‑150 (MQ=255)
                                                               cTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGa                                                                                   <  2:931666/150‑1 (MQ=255)
                                                                ttCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAAt                                                                                 >  2:838953/1‑151 (MQ=255)
                                                                 tCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGaa                                                                                  <  1:270727/149‑1 (MQ=255)
                                                                     ttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCa                                                                             <  2:447673/150‑1 (MQ=255)
                                                                       gTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATg                                                                          >  2:491218/1‑151 (MQ=255)
                                                                        ttATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATg                                                                          <  1:492832/150‑1 (MQ=255)
                                                                         tATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATgg                                                                         >  1:811457/1‑150 (MQ=255)
                                                                          aTGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGAt                                                                       >  1:375768/1‑151 (MQ=255)
                                                                           tGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGAt                                                                       >  1:485911/1‑150 (MQ=255)
                                                                           tGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGt                                                                     >  1:396355/1‑151 (MQ=255)
                                                                                tATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGtt                                                                    <  2:9621/148‑1 (MQ=255)
                                                                                 aTCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa                                                                 >  1:595297/1‑150 (MQ=255)
                                                                                 aTCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa                                                                 <  1:178022/150‑1 (MQ=255)
                                                                                 aTCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAAc                                                                >  2:489346/1‑151 (MQ=255)
                                                                                   cccACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaa                                                              <  1:170319/151‑1 (MQ=255)
                                                                                      aCGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAAt                                                           <  2:814236/151‑1 (MQ=255)
                                                                                            tttAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAtt                                                     >  2:329549/1‑151 (MQ=255)
                                                                                            tttAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAtt                                                     >  2:186690/1‑151 (MQ=255)
                                                                                               aaTATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTtgt                                                 <  2:870287/151‑1 (MQ=255)
                                                                                                atatGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTtgt                                                 <  2:940891/151‑1 (MQ=255)
                                                                                                atatGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTtgt                                                 >  2:702681/1‑151 (MQ=255)
                                                                                                 tatGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTtgtg                                                >  2:404144/1‑151 (MQ=255)
                                                                                                  atGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGa                                               <  1:1112079/151‑1 (MQ=255)
                                                                                                        cATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCaa                                         <  1:107729/151‑1 (MQ=255)
                                                                                                           tAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATa                                       <  1:484464/150‑1 (MQ=255)
                                                                                                              gATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAatt                                   <  1:927039/151‑1 (MQ=255)
                                                                                                               aTGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAatta                                  <  2:824353/151‑1 (MQ=255)
                                                                                                               aTGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAatt                                   >  2:1110292/1‑150 (MQ=255)
                                                                                                                tGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAattat                                 >  1:714212/1‑151 (MQ=255)
                                                                                                                  tttATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAattatt                                >  2:490374/1‑150 (MQ=255)
                                                                                                                   ttATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTg                               <  1:591272/150‑1 (MQ=255)
                                                                                                                    tATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAattatt                                <  1:1030657/148‑1 (MQ=255)
                                                                                                                    tATTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGAc                             >  1:566049/1‑151 (MQ=255)
                                                                                                                     aTTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACt                            >  2:84593/1‑151 (MQ=255)
                                                                                                                     aTTTCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACt                            >  1:120926/1‑151 (MQ=255)
                                                                                                                       ttCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTc                          >  2:1045840/1‑151 (MQ=255)
                                                                                                                        tCTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTc                          <  2:826880/150‑1 (MQ=255)
                                                                                                                         cTTGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGt                        >  2:1173916/1‑151 (MQ=255)
                                                                                                                           tGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCg                         <  1:368971/148‑1 (MQ=255)
                                                                                                                           tGATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGtt                       <  2:1147255/150‑1 (MQ=255)
                                                                                                                            gATTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGtttt                     <  1:863609/151‑1 (MQ=255)
                                                                                                                            gATTTTGCATATGAGTATATT‑A‑‑CCCCCCCCTCCAAAAAAAAAAATAATTTAAAATG‑TGGCTTAATTAAAAAAAAAATTTAAACCAAGAAAT‑TCCAAGGGTGTTTAAAAAAATGGGAATTTTTGGAAAGAATAAatttattaatttggttt                     >  2:1180860/1‑136 (MQ=255)
                                                                                                                             aTTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGtttt                     >  1:981912/1‑150 (MQ=255)
                                                                                                                             aTTTTGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGtttt                     >  2:280934/1‑150 (MQ=255)
                                                                                                                              ttttGCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGtttt                     >  2:1071234/1‑149 (MQ=255)
                                                                                                                                  gCATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGAGGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAAAttt               >  2:486711/1‑151 (MQ=255)
                                                                                                                                   cATATGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGAttt               >  1:836309/1‑150 (MQ=255)
                                                                                                                                    ataaGAGTATATT‑ACCCCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTg              <  1:454624/151‑1 (MQ=255)
                                                                                                                                      atGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGttt           <  1:571057/151‑1 (MQ=255)
                                                                                                                                       tGAGTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGttt           >  1:605499/1‑150 (MQ=255)
                                                                                                                                          gTATATT‑AC‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGc        <  2:425439/150‑1 (MQ=255)
                                                                                                                                               tt‑A‑‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAAt  <  2:396449/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑CCCCCCCCT‑CAAAAAAATAAATTAATTAAAATGATGGCTT‑ATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAAT  >  NC_000913/1211935‑1212227

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