Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 514,014 G→T 8.4% intergenic (+45/+32) 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,0140GT 8.4% 235.4 / 5.2 77intergenic (+45/+32)yhfA/PP_0426hypothetical protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (36/34);  new base T (3/4);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTGGAGCTGTCGGCGGAAAAGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGATAACCC  >  NC_002947/513868‑514158
                                                                                                                                                  |                                                                                                                                                
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                                                                                                                              >  2:797746/1‑150 (MQ=255)
                 ggAAAAGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGc                                                                                                                              >  1:1934546/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                  >  2:904128/1‑143 (MQ=255)
                                                                                                                                    ttGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAgg                  <  1:904128/143‑1 (MQ=255)
                                                                                                                                    ttGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGcc           <  2:2123155/150‑1 (MQ=255)
                                                                                                                                       tgGGAGCGGGTTCACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGa        >  2:236090/1‑150 (MQ=255)
                                                                                                                                         ggaGCGGGTTCACCCGCTTCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGctt                                                               <  1:122063/91‑1 (MQ=255)
                                                                                                                                         ggaGCGGGTTCACCCGCTTCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGctt                                                               >  2:122063/3‑93 (MQ=255)
                                                                                                                                         gAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGATa      >  2:426938/1‑150 (MQ=255)
                                                                                                                                         gAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGATa      >  1:1726352/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)
                                                                                                                                                  |                                                                                                                                                
GCGGTGGAGCTGTCGGCGGAAAAGTATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGCCAGATAACCC  >  NC_002947/513868‑514158

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: