Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,426,139 G→C 6.0% W209C (TGG→TGC) ‡ PP_4765 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,426,1390GC 6.0% 254.8 / 6.9 83W209C (TGG→TGC) ‡PP_4765hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (44/34);  new base C (3/2);  total (47/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

GTGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCACCGGTCC  >  NC_002947/5425990‑5426287
                                                                                                                                                     |                                                                                                                                                    
gTGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTgg                                                                                                                                                      >  2:33480/1‑150 (MQ=255)
 tGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGa                                                                                                                                                     >  2:1320083/1‑150 (MQ=255)
  gAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCGGGCAACCGCTGATGGCAGAACGGGCTGTGGAGATCGCTCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAg                                                                                                                                                    <  1:37164/150‑1 (MQ=255)
    aaGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcg                                                                                                                                                  >  2:322278/1‑150 (MQ=255)
         ttggTCGTGGTGCTGGCGTTGGGTGGTGCCAGCTGGGGGCGCCTGGGTTCTGACGCCAGCCGGCAACCGCTGCTGGCAGAACGGGCTGGGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTGGAAGTAGCGGGCTGGAgcgcgcaa                                                                                                                                             >  2:1882779/3‑148 (MQ=255)
           ggTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACt                                                                                                                                           <  1:40856/150‑1 (MQ=255)
              cGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCaa                                                                                                                                        <  2:603299/150‑1 (MQ=255)
                    gCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACaa                                                                                                                                  <  2:83478/150‑1 (MQ=255)
                          gcTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTcgc                                                                                                                            <  1:761066/150‑1 (MQ=255)
                          gcTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTACttcaaaa                                                                                                                                     <  1:649863/141‑7 (MQ=255)
                          gcTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTACttcaaaa                                                                                                                                     >  2:649863/1‑135 (MQ=255)
                           gtGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGGGCGCGCTTCTCAAGGTCAAGTTcgcg                                                                                                                           <  2:397033/149‑1 (MQ=255)
                           cTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTcgcg                                                                                                                           <  2:1077109/150‑1 (MQ=255)
                            tGGGTGGTGCCAGCTGGGCGCGCGTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGACATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCgg                                                                                                                          >  2:946488/1‑150 (MQ=255)
                            tGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCgg                                                                                                                          <  2:1395779/150‑1 (MQ=255)
                            tGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCgg                                                                                                                          <  1:33480/150‑1 (MQ=255)
                                  gtgCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCAc                                                                                                                    >  1:520637/1‑150 (MQ=255)
                                     ccAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGc                                                                                                                 >  2:1351781/1‑150 (MQ=255)
                                     ccAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGc                                                                                                                 >  1:157165/1‑150 (MQ=255)
                                      cAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCt                                                                                                                >  1:1049234/1‑150 (MQ=255)
                                      cAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCt                                                                                                                >  1:1320231/1‑150 (MQ=255)
                                        gCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCTCTACTCAAGGACAAGTTCGCGGGCGCACCGCTga                                                                                                              <  1:447646/150‑1 (MQ=255)
                                        gCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTga                                                                                                              <  2:622908/150‑1 (MQ=255)
                                         cTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACAGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTgaa                                                                                                             <  2:1001163/150‑1 (MQ=255)
                                                 gcCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCg                                                                                                     <  2:769366/150‑1 (MQ=255)
                                                  ccTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCt                                                                                                                                              >  1:85356/1‑108 (MQ=255)
                                                  ccTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCt                                                                                                                                              <  2:85356/108‑1 (MQ=255)
                                                   cTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGcc                                                                                                   <  1:932516/150‑1 (MQ=255)
                                                     gggATCTGACGGCAGCTGGCAACAGCTGCTGGCAGAACGGTCTGTGGATATCTAGCATTTGAAGCCGAGCAAATGTGACTCTGAATTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCt                                                                                                 >  1:232093/1‑150 (MQ=255)
                                                          cTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCg                                                                                            >  1:914371/1‑150 (MQ=255)
                                                           tGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGt                                                                                           >  1:829528/1‑150 (MQ=255)
                                                             aCGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTga                                                                                                              <  1:549458/129‑1 (MQ=255)
                                                             aCGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTga                                                                                                              >  2:549458/1‑129 (MQ=255)
                                                                     tGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTAGCGGGCGCACCGCTGAAGGACATCGCCCTCAGCGTCCCGGGCAgc                                                                                 >  1:126431/1‑150 (MQ=255)
                                                                         aaCCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAAc                                                                                                         >  1:1276691/1‑122 (MQ=255)
                                                                         aaCCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAAc                                                                                                         <  2:1276691/122‑1 (MQ=255)
                                                                         aaCCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCAc                                                                             >  2:1086923/1‑150 (MQ=255)
                                                                               tgctgGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCa                                                                       >  1:205051/1‑150 (MQ=255)
                                                                               tgctgGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCCACCGTGGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCa                                                                       >  1:511747/1‑150 (MQ=255)
                                                                                  tgGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAgt                                                                    <  1:1320083/150‑2 (MQ=255)
                                                                                   ggCAGAAGGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAggg                                                                   <  2:1423815/150‑1 (MQ=255)
                                                                                    gCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGc                                                                  >  2:1886934/1‑150 (MQ=255)
                                                                                    gCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGc                                                                  >  1:1442719/1‑150 (MQ=255)
                                                                                      aGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCCCCCTCAGCGTCCCTGTCAGCGCACCGCGCAAGGGCGa                                                                >  2:781577/1‑150 (MQ=255)
                                                                                             gCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCAtt                                                         >  1:1905590/1‑150 (MQ=255)
                                                                                             gCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCACACCGCGCAAGGGCGAGTTCAtt                                                         >  1:885718/1‑150 (MQ=255)
                                                                                               tgtgGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATtct                                                       >  1:1149130/1‑150 (MQ=255)
                                                                                               tgtgGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTCTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATtct                                                       >  2:1500577/1‑150 (MQ=255)
                                                                                                 tgGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGGTCGGGGGCCAACCGCTGAAGAACATCTCCCTCAGCGTCCCTGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCa                                                     >  1:224734/1‑150 (MQ=255)
                                                                                                   gATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCTCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCAcc                                                   >  2:401526/1‑150 (MQ=255)
                                                                                                       tctcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGa                                                                                                                                     >  2:578601/1‑64 (MQ=255)
                                                                                                       tctcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGa                                                                                                                                     <  1:578601/64‑1 (MQ=255)
                                                                                                       tctcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTc                                               >  1:466495/1‑150 (MQ=255)
                                                                                                         tcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCggg                                                                                                                         <  1:576222/74‑1 (MQ=255)
                                                                                                         tcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCggg                                                                                                                         >  2:576222/1‑74 (MQ=255)
                                                                                                         tcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCaa                                             >  1:1544324/1‑150 (MQ=255)
                                                                                                          cGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAg                                            <  2:466495/150‑1 (MQ=255)
                                                                                                            ccTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGc                                          >  1:406463/1‑150 (MQ=255)
                                                                                                                tGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGtt                                                                                                                               <  2:185823/61‑1 (MQ=255)
                                                                                                                tGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGtt                                                                                                                               >  1:185823/1‑61 (MQ=255)
                                                                                                                   aGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGaa                                    >  1:1401995/1‑149 (MQ=255)
                                                                                                                   aGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGaa                                    <  2:1401995/149‑1 (MQ=255)
                                                                                                                   aGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAACGGCGCGTTCATTCTCACCGCTCAAGGCGGGGAAg                                   >  2:658153/1‑150 (MQ=255)
                                                                                                                   aGCCGAGCAACTGCGGGGTTGAAGTAGCGGGCTGGAGCGCGCTACTCAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAgg                                  >  2:1368844/1‑150 (MQ=255)
                                                                                                                           aaCTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTg                           <  1:1351781/150‑1 (MQ=255)
                                                                                                                           aaCTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTg                           <  1:946488/150‑1 (MQ=255)
                                                                                                                             cTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGgt                         <  1:42514/150‑1 (MQ=255)
                                                                                                                             cTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGgt                         <  2:1228828/150‑1 (MQ=255)
                                                                                                                              tGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGgtg                        <  2:157165/150‑1 (MQ=255)
                                                                                                                                 ggTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTAt                     <  1:781577/150‑1 (MQ=255)
                                                                                                                                 ggTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTAt                     >  1:1922825/1‑150 (MQ=255)
                                                                                                                                   ttttGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACAGCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGc                   >  2:1698882/1‑150 (MQ=255)
                                                                                                                                   tccttgaGTAGCGCGCTCCAGCCCGCTACTAAAGGACAAGTTAGCGGGCGCACAGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCTCACCGCGCAAGGGCGAGTTCATTCTCACCGATCAAGGCGTGGAGGGCCGCCTGGTGGATgt                   >  2:408271/7‑149 (MQ=255)
                                                                                                                                     ttGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGGGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCGTTCTCACCGTTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCt                 <  1:1882779/150‑1 (MQ=255)
                                                                                                                                     ttGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCt                 <  1:1698882/150‑1 (MQ=255)
                                                                                                                                     ttGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCt                 >  2:1528577/1‑150 (MQ=255)
                                                                                                                                     cttgaGTAGCGCGCTCCAGCCCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAg                                                               <  2:673687/100‑1 (MQ=38)
                                                                                                                                     cttgaGTAGCGCGCTCCAGCCCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAg                                                               >  1:673687/5‑104 (MQ=38)
                                                                                                                                       gAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTgg               >  2:1581285/1‑150 (MQ=255)
                                                                                                                                        aaGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGt              >  1:702331/1‑150 (MQ=255)
                                                                                                                                        aaGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGt              >  2:1372178/1‑150 (MQ=255)
                                                                                                                                        aaGTAGCGGGCTGGAGCGCGCTACGCAAGGACAAGTTCGCGGGCGGACGGCTGAACAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGt              >  1:521823/1‑150 (MQ=255)
                                                                                                                                         aGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTc             <  2:1905590/150‑1 (MQ=255)
                                                                                                                                                ggCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCACCg      <  2:520637/150‑1 (MQ=255)
                                                                                                                                                    ggAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCACCGGTcc  <  1:775494/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCACCGGTCC  >  NC_002947/5425990‑5426287

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: