Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,147,548 T→G 5.9% intergenic (+8/‑73) rpmF → / → plsX 50S ribosomal subunit protein L32/putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,147,5480TG 5.9% 274.2 / 3.9 90intergenic (+8/‑73)rpmF/plsX50S ribosomal subunit protein L32/putative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (47/37);  new base G (2/4);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATGCGTCGTTCCCATGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGC  >  NC_000913/1147406‑1147696
                                                                                                                                              |                                                                                                                                                    
ggCATGCGTCGTTCCCATGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGt                                                                                                                                                <  1:525012/149‑1 (MQ=255)
  cATGCGTCGTTCCCATGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGtga                                                                                                                                              <  2:195550/149‑1 (MQ=255)
   aTGCGTCGTTCCCATGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGtgat                                                                                                                                             <  1:1302007/149‑1 (MQ=255)
   aTGCGTCGTTCCCATGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGtgat                                                                                                                                             <  1:454475/149‑1 (MQ=255)
      cgtcgtTCCCATGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGaa                                                                                                                                          >  2:382359/1‑149 (MQ=255)
      cgtcgtTCCCATGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCAGATGCGtgattac                                                                                                                                          <  1:1432954/149‑4 (MQ=255)
               cATGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGa                                                                                                                                 >  2:1592720/1‑149 (MQ=255)
                aTGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGGCAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAg                                                                                                                                <  1:769792/149‑1 (MQ=255)
                aTGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGGCAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAg                                                                                                                                <  1:1315511/149‑1 (MQ=255)
                aTGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAg                                                                                                                                <  2:1443113/149‑1 (MQ=255)
                  gACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGCCGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGa                                                                                                                              >  1:1048931/1‑149 (MQ=255)
                   aCGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGAt                                                                                                                             <  2:1302583/149‑1 (MQ=255)
                      cgcTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGAtttt                                                                                                                          >  2:229906/1‑149 (MQ=255)
                        cTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTcc                                                                                                                        >  2:448191/1‑149 (MQ=255)
                         tGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTccc                                                                                                                       <  2:1372056/149‑1 (MQ=255)
                           aCCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCCCCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCa                                                                                                                     >  2:1432217/1‑149 (MQ=255)
                            ccGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAg                                                                                                                    <  2:1263953/149‑1 (MQ=255)
                            ccGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGACTTTCCCCAg                                                                                                                    >  2:462512/1‑149 (MQ=255)
                             cGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAgg                                                                                                                   >  1:727309/1‑149 (MQ=255)
                                   cACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTg                                                                                                             <  1:160960/149‑1 (MQ=255)
                                    aCCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACCTCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTgg                                                                                                            >  2:233902/1‑149 (MQ=255)
                                      cAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTgggg                                                                                                          >  1:1567017/1‑149 (MQ=255)
                                      cAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTgggg                                                                                                          >  2:772910/1‑149 (MQ=255)
                                      cAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGAGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTgggg                                                                                                          >  2:1275704/1‑149 (MQ=255)
                                         ccTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGTTGCCTGGGGaaa                                                                                                       <  2:864193/149‑1 (MQ=255)
                                         ccTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGaaa                                                                                                       >  2:1134128/1‑149 (MQ=255)
                                               tgtAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCac                                                                                                 >  1:785953/1‑148 (MQ=255)
                                                 tAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAAc                                                                                               >  2:847937/1‑149 (MQ=255)
                                                      aaaaCTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGggc                                                                                          >  1:294510/1‑149 (MQ=255)
                                                      aaaaCTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGggc                                                                                          >  1:1466083/1‑149 (MQ=255)
                                                      aaaaCTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTAACGCGGCCGCAAGGTCTTCTCTAAGTAATCACGCATCTGAGTGATGAAGGTTTGTGAGGATTTTCCAGAGGGATCTGGGGAAAGACCAAACCgggg                                                                                          >  2:308931/1‑148 (MQ=255)
                                                        aaCTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCCAGGTCATCGCTACGTAATCACGCAGCTGCGTGATGAAGCTTAGTGAGGACTTTCCCCAGGCAACTGGGGAAAGACCAAACCGggcgg                                                                                        <  1:62617/149‑1 (MQ=255)
                                                             cTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGcgacg                                                                                   >  2:576969/1‑149 (MQ=255)
                                                                gTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATa                                                                                >  1:827253/1‑149 (MQ=255)
                                                                gTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTAATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATa                                                                                >  2:1006632/1‑149 (MQ=255)
                                                                  gAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATAcc                                                                              <  2:294510/149‑1 (MQ=255)
                                                                    aaaaCACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCtt                                                                            <  1:100744/149‑1 (MQ=255)
                                                                     aaaCACCTGCGTCACCACATCACTGCCGACGGTTACTCCCGCGGCAGCAAGGTCATCGCTAAGTAATCCCGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCtag                                                                           >  1:66231/1‑147 (MQ=255)
                                                                     aaaCACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGCTACCTTg                                                                           >  2:1418208/1‑149 (MQ=255)
                                                                      cacagCTGCGTCACCACATCAGTGCCGACGGTTACTACCGCGGGCGCAAGGTCATCGCTAAGTATTCACGCATCTGCGTGATGAAGCTTAGTGAGGTTTTTCCCCAGGCAACTGGGGAAAGACCATTCCGGGCGGCGACGATACCTTgt                                                                          <  1:1080624/148‑2 (MQ=255)
                                                                      aaCACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGa                                                                          <  1:1592720/149‑1 (MQ=255)
                                                                         acCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGacac                                                                       <  1:772910/149‑1 (MQ=255)
                                                                           cTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGt                                                                     <  2:1033787/149‑1 (MQ=255)
                                                                                tCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAAc                                                                >  1:747971/1‑149 (MQ=255)
                                                                                tCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAAc                                                                <  1:960070/149‑1 (MQ=255)
                                                                                tCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAAc                                                                <  1:1309895/149‑1 (MQ=255)
                                                                                  accacATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAAccc                                                              >  2:1193352/1‑149 (MQ=255)
                                                                                    cacaTCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTg                                                            <  1:542684/149‑1 (MQ=255)
                                                                                    cacaTCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTg                                                            >  2:396137/1‑149 (MQ=255)
                                                                                    cacaTCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTg                                                            >  2:964111/1‑149 (MQ=255)
                                                                                     acaTCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCATCTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTgg                                                           >  1:688968/1‑149 (MQ=255)
                                                                                      caTCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCACACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGc                                                          >  1:1164783/1‑149 (MQ=255)
                                                                                        tCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGt                                                        <  2:24692/149‑1 (MQ=255)
                                                                                        tCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGt                                                        <  2:1311858/149‑1 (MQ=255)
                                                                                        tCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGt                                                        <  2:1317833/149‑1 (MQ=255)
                                                                                             gCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGAt                                                   <  1:576969/149‑1 (MQ=255)
                                                                                              ccGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATg                                                  >  2:1570862/1‑149 (MQ=255)
                                                                                              ccGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATg                                                  >  1:521695/1‑149 (MQ=255)
                                                                                              ccGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATg                                                  <  2:5740/149‑1 (MQ=255)
                                                                                                gACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTc                                                >  2:1001123/1‑149 (MQ=255)
                                                                                                 aCGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCa                                               >  2:180451/1‑149 (MQ=255)
                                                                                                   ggTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATg                                             <  1:382359/149‑1 (MQ=255)
                                                                                                    gtaactacCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATgg                                            <  1:233902/146‑1 (MQ=255)
                                                                                                    gTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATgg                                            <  2:827253/149‑1 (MQ=255)
                                                                                                     ttACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGCCCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATggg                                           >  2:1539171/1‑149 (MQ=255)
                                                                                                        ctacCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAgg                                        >  2:64844/1‑149 (MQ=255)
                                                                                                           ccGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGa                                     >  1:1209810/1‑149 (MQ=255)
                                                                                                                  cGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGc                              >  1:397634/1‑149 (MQ=255)
                                                                                                                   gCAAGGGCATCGCTAACTAATCACGCCTCTGCGTGATGAAGCTTAGTGCGGAGTGTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGCTACCTTGACACGGCTAACCCTGGCGTTAGATGGCATGGGAGGGGATTTTGGcc                             <  1:98362/149‑1 (MQ=255)
                                                                                                                     aaGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGAGTTTGGCCCt                           >  1:275232/1‑149 (MQ=255)
                                                                                                                            tCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGa                    <  1:448191/149‑1 (MQ=255)
                                                                                                                            tCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGa                    <  1:1462770/149‑1 (MQ=255)
                                                                                                                              gCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACa                  >  2:868141/1‑149 (MQ=255)
                                                                                                                                tAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGt                >  2:465080/1‑149 (MQ=255)
                                                                                                                                  aGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGc              <  1:1470259/149‑1 (MQ=255)
                                                                                                                                  aGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGc              <  2:1341146/149‑1 (MQ=255)
                                                                                                                                   ttcaTCACGCAGATGCGTGATTAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGa                                                                          <  2:197113/85‑1 (MQ=39)
                                                                                                                                   ttcaTCACGCAGATGCGTGATTAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGa                                                                          >  1:197113/4‑88 (MQ=39)
                                                                                                                                   gTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCACACCGGGCGGCGACGATGCCTTGACCCGTCAAACCCTGGCGTTAGATGTCTTGGGAGGGGATTTTGGCCCTTCCGTGACAGTGcg             >  2:122508/1‑148 (MQ=255)
                                                                                                                                   gTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATCTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGcc             >  2:314015/1‑149 (MQ=255)
                                                                                                                                    tAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCt            <  1:229906/149‑1 (MQ=255)
                                                                                                                                    tAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCt            >  1:1026513/1‑149 (MQ=255)
                                                                                                                                     caTCACGCAGATGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTg           >  1:1179066/2‑149 (MQ=255)
                                                                                                                                      aTCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgc          <  2:877261/149‑1 (MQ=255)
                                                                                                                                       tCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgca         <  2:1209810/149‑1 (MQ=255)
                                                                                                                                        cACGCAGCTGCGTGATGAAGCTTAGGGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGGTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgcag        <  1:810169/149‑1 (MQ=255)
                                                                                                                                         aCGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgcagc       >  2:571789/1‑149 (MQ=255)
                                                                                                                                         aCGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgcagc       >  1:1392624/1‑149 (MQ=255)
                                                                                                                                         aCGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgcagc       >  1:775056/1‑149 (MQ=255)
                                                                                                                                         aCGCATCTGCGTGATGAAGCTTAGTGAGGATATTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgcagc       >  2:315019/1‑149 (MQ=255)
                                                                                                                                              tCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGc  >  1:1348485/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GGCATGCGTCGTTCCCATGACGCGCTGACCGCAGTCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGC  >  NC_000913/1147406‑1147696

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