Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,491,540 C→G 5.6% intergenic (‑28/+24) PP_5715 ← / ← PP_5716 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,491,5400CG 5.6% 221.6 / 4.2 71intergenic (‑28/+24)PP_5715/PP_5716hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/27);  new base G (2/2);  total (42/29)
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.45e-01

GGTCAGGTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCAAAATTG  >  NC_002947/5491396‑5491684
                                                                                                                                                |                                                                                                                                                
ggTCAGGTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCATGCTATCGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTCTGCCCTCCAGGCCCGTCGACGCCGTGTc                                                                                                                                             <  1:50433/150‑1 (MQ=255)
        aGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCTGTCGACACGGCGTCGACGGGcc                                                                                                                                     <  1:378762/150‑1 (MQ=255)
           cTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAAc                                                                                                                                  >  2:91943/1‑150 (MQ=255)
             gCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACtt                                                                                                                                >  2:1075175/1‑150 (MQ=255)
                    gAACATCCAACAGTGCTGCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTc                                                                                                                         >  2:1013555/1‑150 (MQ=255)
                     aaCATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCa                                                                                                                        >  2:1198776/1‑150 (MQ=255)
                         tCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGt                                                                                                                    >  1:350611/1‑150 (MQ=255)
                                        gCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCTGTCGACACGGCGTCGACGGGcct                                                                                                                                    <  2:525553/119‑2 (MQ=255)
                                        gCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCTGTCGACACGGCGTCGACGGGcct                                                                                                                                    >  1:525553/1‑118 (MQ=255)
                                        gCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGt                                                                                                              <  1:1014982/141‑1 (MQ=255)
                                        gCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGt                                                                                                              >  2:1014982/1‑141 (MQ=255)
                                         cTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGa                                                                                                    >  2:239501/1‑150 (MQ=255)
                                           cGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGAcc                                                                                                  <  1:1075175/150‑1 (MQ=255)
                                                aGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAg                                                                                             <  1:1013555/150‑1 (MQ=255)
                                                aGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAg                                                                                             >  1:1216298/1‑150 (MQ=255)
                                                     tGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCa                                                                                                                >  1:203303/1‑126 (MQ=255)
                                                     tGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCa                                                                                                                <  2:203303/126‑1 (MQ=255)
                                                     tGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAAATGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCa                                                                                                                >  1:203637/1‑126 (MQ=255)
                                                     tGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAAATGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCa                                                                                                                <  2:203637/126‑1 (MQ=255)
                                                                  gATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGttt                                                                           >  1:728103/1‑150 (MQ=255)
                                                                   aTCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGAGAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTg                                                                          >  1:313251/1‑150 (MQ=255)
                                                                   aTCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTg                                                                          >  2:1286190/1‑150 (MQ=255)
                                                                       aTACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGc                                                                                                            <  1:982/112‑1 (MQ=255)
                                                                       aTACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGc                                                                                                            >  2:982/1‑112 (MQ=255)
                                                                       aTACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTgtcgtc                                                                                 >  2:321485/1‑139 (MQ=255)
                                                                       aTACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTgtcgtc                                                                                 <  1:321485/139‑1 (MQ=255)
                                                                       aTACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGttt                                                                      <  1:91943/150‑1 (MQ=255)
                                                                        tACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCTGAGTTCCAGTGCca                                                                                                          >  1:116624/1‑113 (MQ=255)
                                                                        tACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCTGAGTTCCAGTGCca                                                                                                          <  2:116624/113‑1 (MQ=255)
                                                                        tACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTc                                                                     <  1:1006256/150‑1 (MQ=255)
                                                                           cAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTAc                                                                  >  1:1283142/1‑150 (MQ=255)
                                                                             ggCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACaa                                                                <  1:905556/150‑1 (MQ=255)
                                                                                 aTTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGACGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGgc                                                            <  2:238007/150‑1 (MQ=255)
                                                                                 aTTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGACGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGgc                                                            >  1:238007/1‑150 (MQ=255)
                                                                                             ccACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGACCTTTGTc                                                >  2:142665/1‑150 (MQ=255)
                                                                                             ccACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTc                                                >  1:300216/1‑150 (MQ=255)
                                                                                             ccACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTc                                                <  2:1216298/150‑1 (MQ=255)
                                                                                               acaAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTc                                                >  2:717718/1‑148 (MQ=255)
                                                                                               acaAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTc                                                <  1:717718/148‑1 (MQ=255)
                                                                                               acaAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGt                                              >  1:183778/1‑150 (MQ=255)
                                                                                                caAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTgtc                                                                                    <  2:1106863/111‑1 (MQ=255)
                                                                                                caAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTgtc                                                                                    >  1:1106863/1‑111 (MQ=255)
                                                                                                    gtgtAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTgtcgtc                                                                                 <  2:1025708/110‑1 (MQ=255)
                                                                                                    gtgtAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTgtcgtc                                                                                 >  1:1025708/1‑110 (MQ=255)
                                                                                                    gtgtAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAAGGGCGCGAACTTTGTCGTATcgg                                         <  1:239501/150‑1 (MQ=255)
                                                                                                          ttGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTa                                   <  2:728103/150‑1 (MQ=255)
                                                                                                               ttGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCa                                                                                         >  1:528944/1‑91 (MQ=255)
                                                                                                               ttGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCa                                                                                         <  2:528944/91‑1 (MQ=255)
                                                                                                               ttGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGttt                                                                           <  1:365252/105‑1 (MQ=255)
                                                                                                               ttGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGttt                                                                           >  2:365252/1‑105 (MQ=255)
                                                                                                                   aTCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAg                                                                                                                                        >  2:428961/1‑40 (MQ=255)
                                                                                                                   aTCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAg                                                                                                                                        <  1:428961/40‑1 (MQ=255)
                                                                                                                   aTCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTc                          >  2:1260149/1‑150 (MQ=255)
                                                                                                                     cAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCaa                        <  2:313251/150‑1 (MQ=255)
                                                                                                                       gTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAaca                      <  1:1198776/150‑1 (MQ=255)
                                                                                                                       gTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAaca                      >  1:816081/1‑150 (MQ=255)
                                                                                                                          tGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATGGAAATTGAACAcct                   >  2:1055965/1‑149 (MQ=255)
                                                                                                                            cccTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGCAATTTAACACCGca                 >  2:732874/1‑150 (MQ=255)
                                                                                                                             ccTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGcac                >  1:1041592/1‑150 (MQ=255)
                                                                                                                             ccTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTGCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGGTGAGTTCGGGCGCATTTTCTAACTTTGTCGTATCGGGGGCTATCGAAATTCAACACCGcac                >  2:358904/1‑150 (MQ=255)
                                                                                                                                 cAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACcaca            >  2:352792/1‑150 (MQ=255)
                                                                                                                                   ggCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAgg          >  2:928239/1‑150 (MQ=255)
                                                                                                                                   ggCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAgg          <  2:183778/150‑1 (MQ=255)
                                                                                                                                    gCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCaa                        <  1:776535/135‑1 (MQ=255)
                                                                                                                                    gCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCaa                        >  2:776535/1‑135 (MQ=255)
                                                                                                                                     cctgTCGACACGGCGTCGACGGGCCGACTTCATGGTCAGAGTTCCAATGCCAGCTGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGGACTTTGTCGTATCTGCGTGTCTCGAAATTCAACACCGCACCACAGGcc        >  2:346339/4‑149 (MQ=255)
                                                                                                                                      ccGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCaa       >  2:868112/1‑150 (MQ=255)
                                                                                                                                       cGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCaaa      >  1:319890/1‑150 (MQ=255)
                                                                                                                                         tCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGt                                              <  2:555238/108‑1 (MQ=255)
                                                                                                                                         tCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGt                                              >  1:555238/1‑108 (MQ=255)
                                                                                                                                           gACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCAAAATTg  >  1:393277/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GGTCAGGTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCAAAATTG  >  NC_002947/5491396‑5491684

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: