Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,932,965 G→A 6.5% R186W (CGG→TGG)  PP_3466 ← ABC transporter permease/ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,932,9650GA 6.5% 212.6 / 2.6 73R186W (CGG→TGG) PP_3466ABC transporter permease/ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base A (2/3);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TTGCCCAAGCTGGTCCAGGCGCCCTGCTTCCCGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGCCGATCACCGCCGACACCGCCGAGAACAGCGAGTCCAC  >  NC_002947/3932816‑3933103
                                                                                                                                                     |                                                                                                                                          
ttGCCCAAGCTGGTCCAGGCGCCCTGCTTCCCGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGcgccg                                                                                                                                            >  2:1850364/1‑150 (MQ=255)
   cccAAGCTGGTCCAGGCGCCCTGCTTCCCGGCCTTGCTGGACGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCAt                                                                                                                                         >  1:1694415/1‑150 (MQ=255)
      aaGCTGGTCCAGGCGCCCTGCTTCCCGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCa                                                                                                                                      <  2:1694415/150‑1 (MQ=255)
      aaGCTGGTCCAGGCGCCCTGCTTCCCGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCACATCCa                                                                                                                                      >  1:114268/1‑150 (MQ=255)
       aGCTGGTCCAGGCGCCCTGCTTCCCGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCACATCCAg                                                                                                                                     <  2:114268/150‑1 (MQ=255)
             tCCAGGCGCCCTGCTTGCCTGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGAAAGTTGCGAACTTCGCGCCCCACCTTGCGCCGCATCCAGCGCAGc                                                                                                                               >  2:1349210/1‑150 (MQ=255)
                       cTGCTTCCCGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACtt                                                                                                                     >  1:109921/1‑150 (MQ=255)
                       cTGCTTCCCGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACtt                                                                                                                     >  1:1715239/1‑150 (MQ=255)
                           ttCCCGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGAt                                                                                                                 <  1:1151288/150‑1 (MQ=255)
                           ttCCCGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGAt                                                                                                                 <  1:505078/150‑1 (MQ=255)
                            tCCCGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGc                                                                                                                               <  2:181235/135‑1 (MQ=255)
                            tCCCGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGc                                                                                                                               >  1:181235/1‑135 (MQ=255)
                            tCCCGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATc                                                                                                                <  1:207321/150‑1 (MQ=255)
                             cccGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATc                                                                                                                <  2:847910/149‑1 (MQ=255)
                             cccGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATc                                                                                                                >  1:847910/1‑149 (MQ=255)
                             cccGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCg                                                                                                               <  2:1715239/150‑1 (MQ=255)
                              ccGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCgg                                                                                                              >  1:1459521/1‑150 (MQ=255)
                                            cggccgCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCa                                                                                                >  1:301172/1‑150 (MQ=255)
                                                    tgaatgaaCTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGc                                                                                        <  2:1140933/150‑1 (MQ=255)
                                                      aatgaaCTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGa                                                                                      <  2:858724/150‑1 (MQ=255)
                                                       atgaaCTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGAt                                                                                                                 >  1:1631132/1‑122 (MQ=255)
                                                       atgaaCTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGAt                                                                                                                 <  2:1631132/122‑1 (MQ=255)
                                                           aCTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACt                                                                                 >  2:181210/1‑150 (MQ=255)
                                                           aCTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACt                                                                                 >  2:1304797/1‑150 (MQ=255)
                                                            cTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTg                                                                                >  2:341168/1‑150 (MQ=255)
                                                             ttCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGc                                                                               >  2:366119/1‑150 (MQ=255)
                                                             ttCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGc                                                                               <  1:366119/150‑1 (MQ=255)
                                                                 tCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACg                                                                           >  1:430244/1‑150 (MQ=255)
                                                                    ccGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGAca                                                                        <  2:1459521/150‑1 (MQ=255)
                                                                          gggTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGAcagcatca                                                                  <  1:181210/150‑1 (MQ=255)
                                                                          gggTCACGACCAGGAAAGTCGACACAGCCGCCGAGCGCAGCCGCAAGTTGAGCACTTCGAGCCCCACCTTGCCCCGCATCCAGCGCAGCCACAGCACTTCGAGCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGAcaccatca                                                                  >  2:102324/1‑150 (MQ=255)
                                                                            gTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGAt                                                                                                                 <  1:1267983/101‑1 (MQ=255)
                                                                            gTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGAt                                                                                                                 >  2:1267983/1‑101 (MQ=255)
                                                                            gTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGAcagcatcagc                                                                <  1:1050402/150‑1 (MQ=255)
                                                                             tctcGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGAcagcatcagca                                                               >  1:82767/1‑150 (MQ=255)
                                                                                 gACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCACATCCAGCGCAGCCACAGCACTTCGATCg                                                                                                               <  2:404657/98‑1 (MQ=255)
                                                                                 gACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCACATCCAGCGCAGCCACAGCACTTCGATCg                                                                                                               >  1:404657/1‑98 (MQ=255)
                                                                                  aCCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAAc                                                          <  2:430244/150‑1 (MQ=255)
                                                                                      ggAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACca                                                                                                         >  2:703618/1‑99 (MQ=255)
                                                                                      ggAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACca                                                                                                         <  1:703618/99‑1 (MQ=255)
                                                                                      ggAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCa                                                      <  2:1130835/150‑1 (MQ=255)
                                                                                            acgacACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCAcc                                                    <  1:828259/146‑1 (MQ=255)
                                                                                            acgacACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCAcc                                                    >  2:828259/1‑146 (MQ=255)
                                                                                            acgacACATCCGCCGAGCGCTCACGCAAGTCGCGCACTTCGCGCTCCACCTTGCCCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCACCAGCAACAGCAACAACTGCCACGACAGCCTCAGCATCAACGTCACCgcgc                                                >  2:115809/1‑150 (MQ=255)
                                                                                                   aTCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGAt                                                                                                                 >  1:790809/1‑78 (MQ=255)
                                                                                                   aTCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGAt                                                                                                                 <  2:790809/78‑1 (MQ=255)
                                                                                                   aTCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACa                                                                                          <  1:1258361/101‑1 (MQ=255)
                                                                                                   aTCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACa                                                                                          >  2:1258361/1‑101 (MQ=255)
                                                                                                     ccgccgAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGc                                       <  1:65099/150‑1 (MQ=255)
                                                                                                      cgccgAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACcagca                                                                                                      <  2:1634152/86‑1 (MQ=255)
                                                                                                      cgccgAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACcagca                                                                                                      >  1:1634152/1‑86 (MQ=255)
                                                                                                         cgAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGCCGAt                                   <  2:109921/150‑1 (MQ=255)
                                                                                                         cgAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGCCGAt                                   <  1:109541/150‑1 (MQ=255)
                                                                                                               cTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCGCAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCGTCAGCATCAACGTCACCGCGCCCACCAGGCCGATcaccgc                             <  1:1850364/150‑1 (MQ=255)
                                                                                                                  aCGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCCTCCAGCGCAGCCACAGCACTTCGACCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGCCGATcaccgccga                          >  2:265259/1‑150 (MQ=255)
                                                                                                                  aCGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCAc                                                                            <  1:544291/100‑1 (MQ=255)
                                                                                                                  aCGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCAc                                                                            >  2:544291/1‑100 (MQ=255)
                                                                                                                  aCGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGa                                                                          >  2:539628/1‑102 (MQ=255)
                                                                                                                  aCGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGa                                                                          <  1:539628/102‑1 (MQ=255)
                                                                                                                  aCGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGCCGATcaccgccga                          >  1:558252/1‑150 (MQ=255)
                                                                                                                  aCGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGCCGATcaccgccga                          >  1:1150502/1‑150 (MQ=255)
                                                                                                                   cGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGCCGATcaccgccgac                         >  1:937242/1‑150 (MQ=255)
                                                                                                                         ttGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGCCGATcaccgccga                          >  1:1055462/1‑143 (MQ=255)
                                                                                                                         ttGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGCCGATcaccgccga                          <  2:1055462/143‑1 (MQ=255)
                                                                                                                           gcgcACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGAc                                                                                     >  2:399440/1‑82 (MQ=255)
                                                                                                                           gcgcACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGAc                                                                                     <  1:399440/82‑1 (MQ=255)
                                                                                                                               aCTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGCCGATCACCGCCGACACCGCCGAGAAc             >  1:29599/1‑150 (MQ=255)
                                                                                                                                cTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGa                                                                          <  1:90109/88‑1 (MQ=255)
                                                                                                                                cTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGa                                                                          >  2:90109/1‑88 (MQ=255)
                                                                                                                                 ttCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCGACGTCACCGCGCCCACCAGGCCGATCACCGCCGACACCGCCGAGAACAg           >  1:1823620/1‑150 (MQ=255)
                                                                                                                                     cgcCCCTCCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGCCGATCACCGCCGACACCGCCGAGAACAACGAg       <  1:1349210/150‑1 (MQ=255)
                                                                                                                                       cTCCACCTTGCCCCACATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACCACTGCCACGACAGCATCCGCACCATTGTCACCGCGCCCACCAGGCCGATCACCGCCGACACCGCCGAGAACAGCGAGTc     <  1:830863/150‑1 (MQ=255)
                                                                                                                                          ccccTTCCGCCGCATCCAGCGCACCCACAGCACTCCGATCGGCCCCAGGCTCGCCTTCAATAGCGACCACTGCCACGACAGCCTCAGCATCAACGTCCCCGCGCCCACCAGGCCGATCACCGCCGACACCTCCGAGAACAGCGATTcccc  >  2:272170/3‑148 (MQ=255)
                                                                                                                                          cACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGCCGATCACCGCCGACACCGCCGAGAACAGCGAGTCCAc  >  1:553021/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
TTGCCCAAGCTGGTCCAGGCGCCCTGCTTCCCGGCCTTGCTGGCCGGCCGCCTGAATGAACTTCATCGCCGGCAGGGTCTCGACCAGGAAAGACGACACATCCGCCGAGCGCTCACGCAAGTTGCGCACTTCGCGCTCCACCTTGCGCCGCATCCAGCGCAGCCACAGCACTTCGATCGGCACCAGCAGCGCCAGCAACAGCGACAACTGCCACGACAGCATCAGCATCAACGTCACCGCGCCCACCAGGCCGATCACCGCCGACACCGCCGAGAACAGCGAGTCCAC  >  NC_002947/3932816‑3933103

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