Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,187,090 G→C 7.2% intergenic (+20/+29) pepT → / ← roxA peptidase T/50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,187,0900GC 7.2% 335.0 / 13.0 105intergenic (+20/+29)pepT/roxApeptidase T/50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase
Reads supporting (aligned to +/- strand):  ref base G (47/50);  new base C (4/3);  total (52/53)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 111 total aligned reads displayed.

ACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAGTGC  >  NC_000913/1186941‑1187238
                                                                                                                                                     |                                                                                                                                                    
aCCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtg                                                                                                                                                      <  1:1326200/150‑1 (MQ=255)
aCCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtg                                                                                                                                                      <  2:143726/150‑1 (MQ=255)
  cGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgtg                                                                                                                                                    >  2:1720797/1‑150 (MQ=255)
  cGTGACCGAACCTGTTCACTGTCGGTTAGAACTATCATGTGAAGCATGAGTTTGTGACTCTGGAATGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGAGGCGGCAtgtg                                                                                                                                                    <  2:1533520/150‑1 (MQ=255)
    tGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGAt                                                                                                                                                  <  1:168071/150‑1 (MQ=255)
        cGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATTGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCg                                                                                                                                              <  1:1327572/150‑1 (MQ=255)
        cGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGTATGTGATGCCg                                                                                                                                              <  2:630420/150‑1 (MQ=255)
        cGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCg                                                                                                                                              <  1:1282276/150‑1 (MQ=255)
         gAACCTTTTCACTGGCGGTTAAAACTTTCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCTACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGc                                                                                                                                             <  2:741245/150‑1 (MQ=255)
         gAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGc                                                                                                                                             >  2:40121/1‑150 (MQ=255)
            ccTGTTCACTGGCGGTTACAACTATCATGGTAATCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                                                                          <  2:1290277/150‑1 (MQ=255)
            ccTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATc                                                                                                                                          >  2:1380270/1‑150 (MQ=255)
            ccTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                                                                          <  2:1664204/150‑1 (MQ=255)
             cTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATcc                                                                                                                                         >  1:659931/1‑150 (MQ=255)
              tGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCg                                                                                                                                        >  1:585863/1‑150 (MQ=255)
               gTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCgg                                                                                                                                       >  1:538678/1‑150 (MQ=255)
                ttCACTGTCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGACGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGc                                                                                                                                      <  2:1054229/150‑1 (MQ=255)
                ttCACTGGCGGTTACAACTATGATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGGCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGc                                                                                                                                      <  2:356833/150‑1 (MQ=255)
                ttCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGc                                                                                                                                      <  1:1288346/150‑1 (MQ=255)
                       gCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATc                                                                                                                               >  1:1176421/1‑150 (MQ=255)
                       gCGGTAACAACTATCATGGTATGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATc                                                                                                                               <  2:283906/150‑1 (MQ=255)
                         ggTTACAACTATCATGGTAAGCATGAGATTGTGACTCTGGAAGGTATGGAAAAAGCGGGGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGTATGCGGCATGTGATGCCGCATCCGGCTTAAATCCa                                                                                                                             <  2:887894/150‑1 (MQ=255)
                          gTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAACGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCaa                                                                                                                            <  1:1230071/150‑1 (MQ=255)
                            tACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAAc                                                                                                                          >  1:911675/1‑150 (MQ=255)
                            tACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAAc                                                                                                                          >  2:1271768/1‑150 (MQ=255)
                              cAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACtt                                                                                                                        >  1:82617/1‑150 (MQ=255)
                              cAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACtt                                                                                                                        >  2:460527/1‑150 (MQ=255)
                                 cTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTTCAGGTGATCGTCCGTATTGCCGAGATAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAcc                                                                                                                     <  1:335787/150‑1 (MQ=255)
                                      aTGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCg                                                                                                                <  1:897129/150‑1 (MQ=255)
                                      aTGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCg                                                                                                                <  1:1113351/150‑1 (MQ=255)
                                       tGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCga                                                                                                               >  1:1680825/1‑150 (MQ=255)
                                       tGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCga                                                                                                               >  1:394986/1‑150 (MQ=255)
                                             gCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAcc                                                                                                         <  1:1383686/150‑1 (MQ=255)
                                             gCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAcc                                                                                                         <  1:1045457/150‑1 (MQ=255)
                                                tGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAAt                                                                                                      >  1:40684/1‑150 (MQ=255)
                                                 gAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATa                                                                                                     >  1:399777/1‑150 (MQ=255)
                                                   gTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAcc                                                                                                   >  1:96588/1‑150 (MQ=255)
                                                       gtgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCt                                                                                               >  1:1715127/1‑150 (MQ=255)
                                                        tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAAGCCAAACTTACCCTTCGAAGAACCAAGACCCGCTa                                                                                              >  2:318855/1‑150 (MQ=255)
                                                          aCTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGTAGAACCAATACCCGCTAtt                                                                                            <  2:1540018/150‑1 (MQ=255)
                                                           ctctGGAAGCTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTg                                                                                           <  2:1453893/150‑1 (MQ=255)
                                                            tctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGa                                                                                          <  2:1731148/150‑1 (MQ=255)
                                                               ggAAGGTTTGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCa                                                                                       <  1:633767/150‑1 (MQ=255)
                                                               ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCa                                                                                       >  1:612012/1‑150 (MQ=255)
                                                                   ggTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAgcgc                                                                                   >  1:288406/1‑150 (MQ=255)
                                                                   ggTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATATGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAACAgcgc                                                                                   >  2:361347/1‑150 (MQ=255)
                                                                     tATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGcgccg                                                                                 >  2:1109754/1‑150 (MQ=255)
                                                                      aTGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGcgccgc                                                                                >  2:1695116/1‑150 (MQ=255)
                                                                       tGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGAATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCcgcg                                                                               >  2:727542/1‑150 (MQ=255)
                                                                        ggAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGa                                                                              >  1:387501/1‑150 (MQ=255)
                                                                        cgAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGa                                                                              >  1:248447/2‑150 (MQ=255)
                                                                         gAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCTAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAg                                                                             >  2:1298500/1‑150 (MQ=255)
                                                                             aaGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATc                                                                         >  1:27329/1‑150 (MQ=255)
                                                                              aGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTAACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATcg                                                                        <  1:1083917/150‑1 (MQ=255)
                                                                                cGGTTCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATcgcg                                                                      <  2:1294171/150‑1 (MQ=255)
                                                                                  gTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAGCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAg                                                                    <  1:334017/150‑1 (MQ=255)
                                                                                  gTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAg                                                                    >  1:104205/1‑150 (MQ=255)
                                                                                    gCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGa                                                                  <  1:503046/150‑1 (MQ=255)
                                                                                     cAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGaa                                                                 <  2:82617/150‑1 (MQ=255)
                                                                                       ggTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAAtc                                                               >  1:649874/1‑149 (MQ=255)
                                                                                          gATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACg                                                            <  1:1380270/150‑1 (MQ=255)
                                                                                          gATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACg                                                            >  1:1446568/1‑150 (MQ=255)
                                                                                           aTCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGCCCAGCGCCGCGAGCATCGCGAGGAATGACgg                                                           >  1:593718/1‑150 (MQ=255)
                                                                                           aTCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACgg                                                           <  2:1709748/150‑1 (MQ=255)
                                                                                                ccGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCtt                                                      >  1:1064139/1‑150 (MQ=255)
                                                                                                  gTATTGCCGAGTTAACGGCGCAACGGACGTCCGCGAAACGGGCTGCGGCATGGGCTGCCGCAGCCGGCTTAAATCCAAACTTACCCTTCGCAGACCCAATACCCGCTATTGACCAGCGCCGCGAGCAGCGCGAGGAATGACGGATCTTcc                                                    <  2:192848/150‑1 (MQ=255)
                                                                                                  gTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGCCTGATCTTcc                                                    >  2:1234747/1‑150 (MQ=255)
                                                                                                    aTTGCCGAGTTAACGGCGCAACGGAAGTAAGCGATAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAg                                                  <  2:82894/150‑1 (MQ=255)
                                                                                                     ttGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAgc                                                 <  2:1413032/150‑1 (MQ=255)
                                                                                                      tGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAgcg                                                <  2:561033/150‑1 (MQ=255)
                                                                                                      tGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAgcg                                                <  2:1205261/150‑1 (MQ=255)
                                                                                                      tGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAgcg                                                <  1:40121/150‑1 (MQ=255)
                                                                                                        ccGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCa                                              >  2:484249/1‑150 (MQ=255)
                                                                                                         cGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGAATCCGGCTTAAATCCAACATTACCATTCGAAGAACAAATACCCGCTATTGCCCAGCGCCGCGAGCAACGCGAGGAATGACGGATATACCAGCGCaa                                             >  1:757466/1‑149 (MQ=255)
                                                                                                               aaCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGCCCAGCGCCGCGAGCATCGCGAGTAATGACGGATCTTCAAGCGCATCGCCaa                                       >  2:725525/1‑150 (MQ=255)
                                                                                                                aCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCaaa                                      <  1:1109754/150‑1 (MQ=255)
                                                                                                                 cGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCGTCCGGCTTAAATCCAAACTTACCATTCGAAAAACCACTACCCGCTATTGACCAGCGCCGCGAGATACCCGAGTAATGACGTATCTTCCAGCGCATAGCCaaaa                                     >  2:93531/1‑150 (MQ=255)
                                                                                                                  ggCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAAt                                    <  2:109004/150‑1 (MQ=255)
                                                                                                                   gcgcAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAAtt                                   >  1:1592121/1‑150 (MQ=255)
                                                                                                                     gcAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATtct                                 <  2:245427/150‑1 (MQ=255)
                                                                                                                     gcAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATtct                                 <  2:399777/150‑1 (MQ=255)
                                                                                                                      cAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATtctc                                >  2:1251636/1‑150 (MQ=255)
                                                                                                                         cGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGc                             <  1:460527/150‑1 (MQ=255)
                                                                                                                          ggAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCa                            <  2:911675/150‑1 (MQ=255)
                                                                                                                             aGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTc                         <  1:993170/150‑1 (MQ=255)
                                                                                                                               tAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAg                       >  2:664119/1‑150 (MQ=255)
                                                                                                                               tAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAg                       <  2:14849/150‑1 (MQ=255)
                                                                                                                                 aGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAgcg                     <  1:878838/150‑1 (MQ=255)
                                                                                                                                 aGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAgcg                     <  1:1113605/150‑1 (MQ=255)
                                                                                                                                     aaaaGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAAt                 >  2:1225266/1‑150 (MQ=255)
                                                                                                                                      aaaGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATg                >  1:1424949/1‑150 (MQ=255)
                                                                                                                                       gccggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGt               >  2:440509/4‑150 (MQ=255)
                                                                                                                                       gccggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGt               >  1:1089857/4‑150 (MQ=255)
                                                                                                                                       gccggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGt               >  1:880552/4‑150 (MQ=255)
                                                                                                                                       aaGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGt               <  1:727542/150‑1 (MQ=255)
                                                                                                                                         gggATGCGGCATGTGATGTCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTg             <  2:1382756/150‑1 (MQ=255)
                                                                                                                                          ggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGc            >  2:1371577/1‑150 (MQ=255)
                                                                                                                                           gATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCt           >  1:1029990/1‑150 (MQ=255)
                                                                                                                                                ggCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGa      >  1:501307/1‑150 (MQ=255)
                                                                                                                                                    tgtgATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAGTGc  >  2:1190259/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAGTGC  >  NC_000913/1186941‑1187238

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