Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,378,716 C→A 6.0% L273L (CTG→CTT dnaK ← chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,378,7160CA 6.0% 839.6 / 27.3 249L273L (CTG→CTTdnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (100/134);  new base A (10/5);  total (110/139)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 262 total aligned reads displayed.

GCGATGCGGCAAGGCTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTC  >  NC_002947/5378567‑5378865
                                                                                                                                                     |                                                                                                                                                     
    tGCGGCAAGGCTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCt                                                                                                                                                   <  1:4350870/150‑1 (MQ=255)
      cGGCAAGGCTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCg                                                                                                                                                 <  1:2658492/150‑1 (MQ=255)
        gCAAGGCTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGAt                                                                                                                                               <  1:3657292/150‑1 (MQ=255)
         cAAGGCTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATc                                                                                                                                              <  1:4245517/150‑1 (MQ=255)
         cAAGGCTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATc                                                                                                                                              <  1:4245362/150‑1 (MQ=255)
           aGGCTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCtt                                                                                                                                            <  1:8506820/150‑1 (MQ=255)
            ggCTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCTGTCGTCCGCTcagtcgac                                                                                                                                           <  2:6570211/150‑8 (MQ=255)
            ggCTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTg                                                                                                                                           <  1:4886957/150‑1 (MQ=255)
            ggCTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTg                                                                                                                                           <  2:6898207/150‑1 (MQ=255)
            ggCTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTg                                                                                                                                           <  2:6006051/150‑1 (MQ=255)
              cTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGc                                                                                                                                         <  1:510521/150‑1 (MQ=255)
              cTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGc                                                                                                                                         <  1:284648/150‑1 (MQ=255)
               tCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCt                                                                                                                                        >  2:6173765/1‑150 (MQ=255)
               tCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCt                                                                                                                                        >  1:785675/1‑150 (MQ=255)
                 gATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCttt                                                                                                                                      <  2:8422177/150‑1 (MQ=255)
                 gATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCttt                                                                                                                                      >  2:1622417/1‑150 (MQ=255)
                      tACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTg                                                                                                                                 <  2:5325958/150‑1 (MQ=255)
                      tACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTg                                                                                                                                 <  1:3439212/150‑1 (MQ=255)
                       aCGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGc                                                                                                                                <  2:7537432/150‑1 (MQ=255)
                        cGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCg                                                                                                                               <  2:5525701/150‑1 (MQ=255)
                         gCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCgg                                                                                                                              >  2:157772/1‑150 (MQ=255)
                           ttGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGct                                                                                                                            >  1:1952538/1‑150 (MQ=255)
                             gACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGcttc                                                                                                                          <  1:658817/150‑1 (MQ=255)
                               ccAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGcttctt                                                                                                                        <  2:7863745/150‑1 (MQ=255)
                                 aGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTc                                                                                                                      <  2:6081238/150‑1 (MQ=255)
                                 aGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTc                                                                                                                      <  1:5430858/150‑1 (MQ=255)
                                  ggTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCa                                                                                                                     <  1:9120133/150‑1 (MQ=255)
                                    tCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGa                                                                                                                   <  2:1731234/150‑1 (MQ=255)
                                         gACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTg                                                                                                              >  2:4134214/1‑150 (MQ=255)
                                                  cTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCag                                                                                                     >  2:644176/1‑150 (MQ=255)
                                                   tCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCaga                                                                                                    <  2:6086081/150‑1 (MQ=255)
                                                   tCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCaga                                                                                                    <  2:3482172/150‑1 (MQ=255)
                                                   tCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCaga                                                                                                    >  2:246245/1‑150 (MQ=255)
                                                   gccAGCGTGGCGCGGGAGATCGTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCaga                                                                                                    <  2:3261311/149‑1 (MQ=255)
                                                    ccAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGTGGCTTCTTTCAGACGCTGCAGGGCCagag                                                                                                   >  2:2416494/1‑150 (MQ=255)
                                                      aGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGa                                                                                                 >  1:174949/1‑150 (MQ=255)
                                                       gCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGAt                                                                                                >  2:8880306/1‑150 (MQ=255)
                                                             cgcgGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGttct                                                                                          <  1:1046654/150‑1 (MQ=255)
                                                               cgGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGttcttc                                                                                        <  1:6537531/150‑1 (MQ=255)
                                                               cgGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGttcttc                                                                                        <  2:2743019/150‑1 (MQ=255)
                                                                     aTCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTcc                                                                                  <  2:5715918/150‑1 (MQ=255)
                                                                     aTCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTcc                                                                                  >  2:818083/1‑150 (MQ=255)
                                                                        ttCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATg                                                                               >  2:1462018/1‑150 (MQ=255)
                                                                        ttCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATg                                                                               >  1:6237715/1‑150 (MQ=255)
                                                                        ttCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATg                                                                               <  2:8728592/150‑1 (MQ=255)
                                                                        ttCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATg                                                                               >  1:3708284/1‑150 (MQ=255)
                                                                        ttCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATg                                                                               <  1:5181214/150‑1 (MQ=255)
                                                                         tCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGc                                                                              >  2:6320263/1‑150 (MQ=255)
                                                                          cACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGcc                                                                             >  2:4262028/1‑150 (MQ=255)
                                                                           aCGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGccc                                                                            <  2:2521906/150‑1 (MQ=255)
                                                                             gTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGa                                                                          >  1:1762412/1‑150 (MQ=255)
                                                                             gTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGa                                                                          <  1:2020934/150‑1 (MQ=255)
                                                                                    tGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttc                                                                   <  2:1148816/150‑1 (MQ=255)
                                                                                    tGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttc                                                                   <  1:2234918/150‑1 (MQ=255)
                                                                                     gCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttct                                                                  <  2:5832112/150‑1 (MQ=255)
                                                                                     gCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttct                                                                  >  2:8610004/1‑150 (MQ=255)
                                                                                     gCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttct                                                                  <  2:4615936/150‑1 (MQ=255)
                                                                                     gCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttct                                                                  >  1:1095359/1‑150 (MQ=255)
                                                                                      cTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttctt                                                                 <  1:3939187/150‑1 (MQ=255)
                                                                                       ttAGGACCGGTCTGGTCAGCCGTGATGAATGGCAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCAGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTg                                                                >  1:796676/1‑150 (MQ=255)
                                                                                           gACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTg                                                                >  2:5143285/1‑146 (MQ=255)
                                                                                           gACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTg                                                                <  1:569258/146‑1 (MQ=255)
                                                                                               ggTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcgtc                                                        <  2:777677/150‑1 (MQ=255)
                                                                                                gTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcgtcg                                                       >  2:4944764/1‑150 (MQ=255)
                                                                                                     aTCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAg                                                  >  2:8927512/1‑150 (MQ=255)
                                                                                                        tGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTa                                               <  1:6787380/150‑1 (MQ=255)
                                                                                                        tGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTa                                               <  1:8880306/150‑1 (MQ=255)
                                                                                                        tGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTa                                               >  1:7446571/1‑150 (MQ=255)
                                                                                                           aGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTc                                            >  1:6827978/1‑150 (MQ=255)
                                                                                                           aGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTc                                            >  2:8349439/1‑150 (MQ=255)
                                                                                                             tGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGa                                          <  2:4446767/150‑1 (MQ=255)
                                                                                                                  tACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAgg                                     >  2:1279916/1‑150 (MQ=255)
                                                                                                                  tACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAgg                                     >  2:1681344/1‑150 (MQ=255)
                                                                                                                   aCGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGgc                                    >  1:2574906/1‑150 (MQ=255)
                                                                                                                   aCGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCATGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTAGAACTCGTCGACAAGGTAGACGATCAGgc                                    >  2:3877363/1‑150 (MQ=255)
                                                                                                                    cGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGgcg                                   <  1:9059641/150‑1 (MQ=255)
                                                                                                                      gCAGGTTGAATTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTATCTAAAGATCGATGCCCTGTTCCTTCTTGAACTCGCCAACCAGATATTCGACCAGACGCa                                 >  2:3631120/1‑150 (MQ=255)
                                                                                                                       cAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCATTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCAt                                <  1:4554274/150‑1 (MQ=255)
                                                                                                                          gTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTc                             >  1:7697178/1‑150 (MQ=255)
                                                                                                                           ttgatttcggtctGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTca                            >  2:6519457/7‑149 (MQ=255)
                                                                                                                            tGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGa                           >  2:9036881/1‑150 (MQ=255)
                                                                                                                             gACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGaa                          <  2:2541111/150‑1 (MQ=255)
                                                                                                                             gACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGaa                          <  1:1622416/150‑1 (MQ=255)
                                                                                                                             gACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGaa                          <  1:2489850/150‑1 (MQ=255)
                                                                                                                               tttcggtctGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGt                        >  1:4099803/3‑150 (MQ=255)
                                                                                                                                gtcggtcgACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTc                       >  1:8244597/1‑150 (MQ=255)
                                                                                                                                       gACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATc                                                                                               <  2:3669369/71‑1 (MQ=255)
                                                                                                                                       gACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCa                <  2:5442057/150‑1 (MQ=255)
                                                                                                                                       gACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCa                >  1:4025569/1‑150 (MQ=255)
                                                                                                                                          tGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGccg             >  2:4083978/1‑149 (MQ=255)
                                                                                                                                           gAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGccga            <  1:1500389/150‑3 (MQ=255)
                                                                                                                                               ggACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGaa        <  1:7540294/150‑1 (MQ=255)
                                                                                                                                               ggACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGaa        >  1:574934/1‑150 (MQ=255)
                                                                                                                                               ggACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGaa        <  1:246245/150‑1 (MQ=255)
                                                                                                                                               ggACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGCCGAAGTCTTCGCCACCCAGGaa        <  1:7276940/150‑1 (MQ=255)
                                                                                                                                               ggACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGac        >  2:1992170/1‑149 (MQ=255)
                                                                                                                                                  cgacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGgt     <  1:301483/150‑1 (MQ=255)
                                                                                                                                                   gacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGgtg    >  1:3152471/1‑150 (MQ=255)
                                                                                                                                                   gacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGgtg    >  2:2806260/1‑150 (MQ=255)
                                                                                                                                                     cAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTc  <  1:157772/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGATGCGGCAAGGCTCGATGGTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTC  >  NC_002947/5378567‑5378865

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