Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 514,015 A→C 7.9% intergenic (+46/+31) yhfA → / ← PP_0426 hypothetical protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947514,0150AC 7.9% 207.4 / 9.1 77intergenic (+46/+31)yhfA/PP_0426hypothetical protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base C (3/3);  total (39/38)
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.61e-01

GCGGTGGAGCTGTCGGCGGAAAAGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGATAACCCCCAGCC  >  NC_002947/513868‑514164
                                                                                                                                                   |                                                                                                                                                     
gcggTGGAGCTGTCGGCGGAAAAGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAAc                                                                                                                                                     <  2:1831638/150‑1 (MQ=255)
    gggAGCTGTCGGCGGAAAAGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGc                                                                                                                                                 <  2:187659/149‑1 (MQ=255)
     ggAGCTGTCGGCGGAAAAGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCt                                                                                                                                                <  2:2387349/150‑1 (MQ=255)
         cTGTCGGCGGAAAAGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGGAGCGGGTTCACCCGCTTCCa                                                                                                                                            <  1:479422/150‑1 (MQ=255)
         cTGTCGGCGGAAAAGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCa                                                                                                                                            >  2:885990/1‑150 (MQ=255)
                 ggAAAAGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGc                                                                                                                                    >  1:1934546/1‑150 (MQ=255)
                 ggAAAAGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGc                                                                                                                                    >  2:797746/1‑150 (MQ=255)
                  gAAAAGTATTGCTCGGCTTCGATCATGATCGCGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGTGTGAACCCGCTCCCACCCGCGGcc                                                                                                                                   >  2:374584/1‑150 (MQ=255)
                    aaaGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCAc                                                                                                                                           >  1:566394/1‑140 (MQ=255)
                    aaaGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCAc                                                                                                                                           <  2:566394/140‑1 (MQ=255)
                     aaGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCgg                                                                                                                                <  2:407140/150‑1 (MQ=255)
                       gTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGc                                                                                                                              <  2:1283001/150‑1 (MQ=255)
                         aTTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGa                                                                                                                            <  1:885990/150‑1 (MQ=255)
                          ttGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACACGCGGCCCGGGCGAg                                                                                                                           >  1:668599/1‑150 (MQ=255)
                             cTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGttt                                                                                                                        >  2:2031583/1‑150 (MQ=255)
                             cTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGAGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCGCCCACCCGCGGCCCGGGCGAGttt                                                                                                                        >  2:1210939/1‑150 (MQ=255)
                              ccGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCAAAGGGTATAAAGTAAAAGCTTTGGATAAAAGCGGAAAGGTTGGGGTGGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTc                                                                                                                       >  1:2123155/2‑150 (MQ=255)
                                  cTTCGACCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGtt                                                                                                                   <  1:116491/150‑1 (MQ=255)
                                   ttCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTg                                                                                                                  >  2:869724/1‑150 (MQ=255)
                                          aTGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAGGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAggg                                                                                                           <  2:2173935/150‑1 (MQ=255)
                                            gCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGcc                                                                                                         <  1:566917/150‑1 (MQ=255)
                                              tCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCgg                                                                                                       >  1:1702188/1‑150 (MQ=255)
                                                 aGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGGTTTGCAGCCCTTGTGGAAGCGGGTGAACCGGCTGCCGCCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAg                                                                                                    <  2:1387208/150‑1 (MQ=255)
                                                   cGTGCCGGGGTTGAGACTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAggg                                                                                                  <  1:1633726/150‑1 (MQ=255)
                                                        cGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGt                                                                                                                          >  2:745759/1‑121 (MQ=255)
                                                        cGGGGTTGAGATTACCCATGGCTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGt                                                                                                                          <  1:745759/121‑1 (MQ=255)
                                                            gTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCa                                                                                         >  2:991649/1‑150 (MQ=255)
                                                            gTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCa                                                                                         >  1:637284/1‑150 (MQ=255)
                                                             ttGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCaa                                                                                        >  1:1027583/1‑150 (MQ=255)
                                                                   ttACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAAc                                                                                                                                                     >  1:1746250/1‑83 (MQ=255)
                                                                   ttACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAAc                                                                                                                                                     <  2:1746250/83‑1 (MQ=255)
                                                                    tACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGc                                                                                 <  2:1934546/150‑1 (MQ=255)
                                                                       ccATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGa                                                                              >  1:401106/1‑150 (MQ=255)
                                                                       ccATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGa                                                                              >  2:1334258/1‑150 (MQ=255)
                                                                         aTGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAAcc                                                                                      <  1:1382646/140‑1 (MQ=255)
                                                                         aTGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAAcc                                                                                      >  2:1382646/1‑140 (MQ=255)
                                                                          tGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAgg                                                                           >  2:291507/1‑149 (MQ=255)
                                                                             gTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACgg                                                                        <  1:803084/150‑1 (MQ=255)
                                                                                   aaTCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAGCGGGTTCACCCGCTTCCAcaa                                                                                                                                         >  1:1824716/1‑76 (MQ=25)
                                                                                   aaTCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAGCGGGTTCACCCGCTTCCAcaa                                                                                                                                         <  2:1824716/78‑3 (MQ=25)
                                                                                    aTCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTa                                                                 >  1:381162/1‑150 (MQ=255)
                                                                                        tAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAgcgc                                                             >  1:764663/1‑150 (MQ=255)
                                                                                          gagCTGCGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCt                                                           <  1:991649/150‑1 (MQ=255)
                                                                                            gCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCGGGGCGAGTTTCAGTGGCTCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTgg                                                         <  1:291507/150‑1 (MQ=255)
                                                                                             cTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGc                                                        >  1:1425794/1‑150 (MQ=255)
                                                                                                 gTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCaa                                                    <  1:1466044/150‑1 (MQ=255)
                                                                                                  ttAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAAc                                                   <  1:554365/150‑1 (MQ=255)
                                                                                                    aaGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCg                                                 <  2:637284/150‑1 (MQ=255)
                                                                                                       gCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAATTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTc                                              >  2:873163/1‑150 (MQ=255)
                                                                                                       gCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGAGAGGTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTc                                              >  2:697742/1‑150 (MQ=255)
                                                                                                         tGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTccc                                            >  2:1447265/1‑150 (MQ=255)
                                                                                                                 cGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCttt                                    <  1:869724/150‑1 (MQ=255)
                                                                                                                   gggCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTg                                  <  1:2031583/150‑1 (MQ=255)
                                                                                                                   gggCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTg                                  <  1:544910/150‑1 (MQ=255)
                                                                                                                      ccGCTTTGCAGCCCTTGTGGAAGCGGGTGAANCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGa                                                                                                      <  1:597087/79‑1 (MQ=255)
                                                                                                                      ccGCTGTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGa                                                                                                      >  2:597087/1‑79 (MQ=255)
                                                                                                                       cGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTaa                              <  2:2373747/150‑1 (MQ=255)
                                                                                                                         cTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTgg                                                         >  2:339364/1‑121 (MQ=255)
                                                                                                                         cTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTgg                                                         <  1:339364/121‑1 (MQ=255)
                                                                                                                         cTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAAc                            >  1:1686441/1‑150 (MQ=255)
                                                                                                                          tttGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTg                                  <  1:1049327/143‑1 (MQ=255)
                                                                                                                          tttGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTg                                  >  2:1049327/1‑143 (MQ=255)
                                                                                                                           ttGCAGCCCTTGTGGAAGCGGGTTTACTCGCTTCCACCCGCGGCCGGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTTAACCGCTCCCGGCGCTTTTGGTAACCCa                          <  2:416798/150‑1 (MQ=255)
                                                                                                                             gCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAgg                        <  2:83427/150‑1 (MQ=255)
                                                                                                                              cAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGc                       >  1:1284215/1‑150 (MQ=255)
                                                                                                                                    ttGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAgg                        <  1:904128/143‑1 (MQ=255)
                                                                                                                                    ttGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAgg                        >  2:904128/1‑143 (MQ=255)
                                                                                                                                    ttGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGcc                 <  2:2123155/150‑1 (MQ=255)
                                                                                                                                       tgGGAGCGGGTTCACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGa              >  2:236090/1‑150 (MQ=255)
                                                                                                                                         ggaGCGGGTTCACCCGCTTCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGctt                                                                     >  2:122063/3‑93 (MQ=255)
                                                                                                                                         ggaGCGGGTTCACCCGCTTCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGctt                                                                     <  1:122063/91‑1 (MQ=255)
                                                                                                                                         gAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGATa            >  1:1726352/1‑150 (MQ=255)
                                                                                                                                         gAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGATa            >  2:426938/1‑150 (MQ=255)
                                                                                                                                          aaGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTACACCAACCCGAGCCAGATaa           >  2:2274798/1‑150 (MQ=255)
                                                                                                                                            gCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCGGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCCACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGATAAcc         >  2:122633/1‑150 (MQ=255)
                                                                                                                                            gCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGATAAcc         <  2:1686441/150‑1 (MQ=255)
                                                                                                                                             cGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGATAAccc        <  2:1027583/150‑1 (MQ=255)
                                                                                                                                                   aaCCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGATAACCCCCAGcc  <  1:2290390/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGGTGGAGCTGTCGGCGGAAAAGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGATAACCCCCAGCC  >  NC_002947/513868‑514164

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