Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,579,833 C→G 5.9% V361V (GTG→GTC serB ← phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,579,8330CG 5.9% 272.2 / 5.6 85V361V (GTG→GTCserBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (39/41);  new base G (3/2);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CGATCACGCAGGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCAAGTTCGTTGGCGAACACATAG  >  NC_002947/5579684‑5579980
                                                                                                                                                     |                                                                                                                                                   
cgATCACGCAGGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGccc                                                                                                                                                     >  2:1838454/1‑150 (MQ=255)
       gCAGGCCCAGCATGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCCGCCAGCGACAGCATCGGCAGGTCATTGGCGCCAGCGCCCACGGCGa                                                                                                                                              <  2:242176/150‑1 (MQ=255)
           gCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCGGGGGCGATggc                                                                                                                                          <  2:746540/150‑2 (MQ=255)
            cccAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGAACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTc                                                                                                                                         <  1:1352814/150‑1 (MQ=255)
              cagcagGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTg                                                                                                                                       >  2:1850934/1‑150 (MQ=255)
                  agGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATg                                                                                                                                            >  1:1922286/1‑141 (MQ=255)
                  agGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATg                                                                                                                                            <  2:1922286/141‑1 (MQ=255)
                    gTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGATCTGCTCCAg                                                                                                                                 <  2:2402735/150‑1 (MQ=255)
                      aCAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCATGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGtt                                                                                                                               >  2:1412165/1‑150 (MQ=255)
                        aGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGc                                                                                                                             >  1:2306634/1‑150 (MQ=255)
                          cACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCaa                                                                                                                           >  1:1325343/1‑150 (MQ=255)
                            ccccGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGc                                                                                                                         <  2:2320816/150‑1 (MQ=255)
                            ccccGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGc                                                                                                                         >  2:262875/1‑150 (MQ=255)
                             cccGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATg                                                                                                                                            <  1:2268550/130‑1 (MQ=255)
                             cccGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATg                                                                                                                                            >  2:2268550/1‑130 (MQ=255)
                             cccGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGcc                                                                                                                        >  2:2437968/1‑150 (MQ=255)
                             cccGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCACCATCGCCGTGGGCGATggcgccaa                                                                                                                                     >  2:492180/1‑131 (MQ=255)
                                gTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACGAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTtc                                                                                                                     <  1:135577/150‑1 (MQ=255)
                                     gCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGc                                                                                                                >  2:1530313/1‑150 (MQ=255)
                                       cccAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTg                                                                                                              >  2:2188457/1‑150 (MQ=255)
                                        ccAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTgg                                                                                                             >  2:1040352/1‑150 (MQ=255)
                                                   aTCGCCTGCTTGGCCGACTGGCGCACCAGCGGGTCGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCttcttc                                                                                                  <  1:1163323/150‑1 (MQ=255)
                                                   aTCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGGCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGTCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCttcttc                                                                                                  <  1:2367518/150‑1 (MQ=255)
                                                   aTCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCttcttc                                                                                                  <  1:2015000/150‑1 (MQ=255)
                                                   aTCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCttcttc                                                                                                  <  2:1586668/150‑1 (MQ=255)
                                                       ccTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTcaaca                                                                                              >  1:590824/1‑150 (MQ=255)
                                                           cTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATc                                                                                          <  1:1765809/150‑1 (MQ=255)
                                                               gCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGcc                                                                                      <  2:1521584/150‑1 (MQ=255)
                                                                ccGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCt                                                                                     >  1:790750/1‑150 (MQ=255)
                                                                  gACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTg                                                                                   >  2:1523592/1‑150 (MQ=255)
                                                                        cgcACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTgg                                                                             <  2:286822/150‑1 (MQ=255)
                                                                            ccAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTg                                                                                                              >  2:2319382/1‑113 (MQ=255)
                                                                            ccAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTg                                                                                                              <  1:2319382/113‑1 (MQ=255)
                                                                             cAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATc                                                                        <  1:1040352/150‑1 (MQ=255)
                                                                             cAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATc                                                                        >  1:1290632/1‑150 (MQ=255)
                                                                               gCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTcaa                                                                                                >  2:400000/1‑124 (MQ=255)
                                                                               gCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTcaa                                                                                                <  1:400000/124‑1 (MQ=255)
                                                                                  gCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACaa                                                                   <  1:1850934/150‑1 (MQ=255)
                                                                                   cTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAAt                                                                  <  1:1412165/150‑1 (MQ=255)
                                                                                    ttGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCGACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGa                                                                      <  2:1937095/145‑1 (MQ=255)
                                                                                    ttGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCGACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGa                                                                      >  1:1937095/1‑145 (MQ=255)
                                                                                    ttGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATc                                                                 >  2:1987795/1‑150 (MQ=255)
                                                                                       gCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTc                                                                                                                                    <  2:2074271/80‑1 (MQ=255)
                                                                                       gCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTc                                                                                                                                    >  1:2074271/1‑80 (MQ=255)
                                                                                       gCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGc                                                              >  2:2008837/1‑150 (MQ=255)
                                                                                        cACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCt                                                             <  2:877094/150‑1 (MQ=255)
                                                                                        cACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCCTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCt                                                             >  1:1996472/1‑150 (MQ=255)
                                                                                           gAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCaa                                                          <  1:2437968/150‑1 (MQ=255)
                                                                                           gAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCaa                                                          >  1:393966/1‑150 (MQ=255)
                                                                                            aaaCGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAAc                                                         >  1:1177787/1‑150 (MQ=255)
                                                                                            aaaCGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAAc                                                         <  1:1044337/150‑1 (MQ=255)
                                                                                                gccacgccCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGcc                                                     <  1:1530313/150‑1 (MQ=255)
                                                                                                  cacnccCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCGTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCAc                                                   >  2:1520375/1‑150 (MQ=255)
                                                                                                  cacgccCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCAc                                                   >  1:1394220/1‑150 (MQ=255)
                                                                                                    cgccCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGc                                                 <  1:2188457/150‑1 (MQ=255)
                                                                                                       ccAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGa                                                                      <  2:2351056/126‑1 (MQ=255)
                                                                                                       ccAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGa                                                                      >  1:2351056/1‑126 (MQ=255)
                                                                                                       ccAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCg                                              >  2:961225/1‑150 (MQ=255)
                                                                                                         aaCCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGc                                                                           <  1:1491340/119‑1 (MQ=255)
                                                                                                         aaCCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGc                                                                           >  2:1491340/1‑119 (MQ=255)
                                                                                                         aaCCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTc                                            >  1:1572283/1‑150 (MQ=255)
                                                                                                         aaCCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCAGTc                                            <  2:590824/150‑1 (MQ=255)
                                                                                                             cagccagcGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGTCATCGACAATCGGCTCAACAGCCACGCCCGTCACCt                                        >  1:2440541/1‑150 (MQ=255)
                                                                                                                  agcGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCa                                   >  1:1032436/1‑150 (MQ=255)
                                                                                                                   gcGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCAt                                  >  2:1065848/1‑150 (MQ=255)
                                                                                                                    cGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATc                                 <  1:240017/150‑1 (MQ=255)
                                                                                                                     gACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGACCTTTCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATc                                 <  2:1327237/149‑1 (MQ=255)
                                                                                                                     gACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGACCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATc                                 >  1:1327237/1‑149 (MQ=255)
                                                                                                                       cAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGAc                              <  2:790750/150‑1 (MQ=255)
                                                                                                                         gCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCa                            <  1:961225/150‑1 (MQ=255)
                                                                                                                          cATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCAc                           >  2:1522381/1‑150 (MQ=255)
                                                                                                                           aTCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTtc                                                                                                                     >  1:552153/1‑59 (MQ=255)
                                                                                                                           aTCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTtc                                                                                                                     <  2:552153/59‑1 (MQ=255)
                                                                                                                               gCAGGTCATTGGCGCCATCGCCCACGGGGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATc                                                                        <  2:678560/100‑1 (MQ=255)
                                                                                                                               gCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATc                                                                        >  1:678560/1‑100 (MQ=255)
                                                                                                                               gCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCCTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCaa                      <  2:1996472/150‑1 (MQ=255)
                                                                                                                                   gTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCAAGTTc                  <  1:1520375/150‑1 (MQ=255)
                                                                                                                                    tCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCAAGTTCg                 >  1:1152356/1‑150 (MQ=255)
                                                                                                                                      aTTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCAAGTTCGtt               <  2:1032436/150‑1 (MQ=255)
                                                                                                                                       ttGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCgg                                                                                        >  1:1883948/1‑76 (MQ=255)
                                                                                                                                       ttGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCgg                                                                                        <  2:1883948/76‑1 (MQ=255)
                                                                                                                                       ttGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCt                                                             >  1:698899/1‑103 (MQ=255)
                                                                                                                                       ttGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCt                                                             <  2:698899/103‑1 (MQ=255)
                                                                                                                                          agaccATCGCCGTGGGCGATGGCCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAAc                                                         >  1:1790963/4‑104 (MQ=255)
                                                                                                                                             ccATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCAAGTTCGTTGGCGAac        <  2:1549780/150‑1 (MQ=255)
                                                                                                                                             ccATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTAACGGTCGACCACTTCAAGTTCGTTGGCGAac        <  1:665943/150‑1 (MQ=255)
                                                                                                                                                   cccACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCACGTTCGTTGGCGAACACATAg  >  1:2283995/1‑150 (MQ=255)
                                                                                                                                                    ccACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCaa                      <  1:1711603/129‑1 (MQ=255)
                                                                                                                                                    ccACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCaa                      >  2:1711603/1‑129 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGATCACGCAGGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCAAGTTCGTTGGCGAACACATAG  >  NC_002947/5579684‑5579980

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