Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,704,679 T→A 15.9% M1* (TTG→TAG) †‡ PP_1499 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,704,6790TA15.9% 222.4 / 35.2 88M1* (TTG→TAG) †‡PP_1499hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base A (9/5);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.48e-01

CGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGACGACAAC  >  NC_002947/1704531‑1704827
                                                                                                                                                    |                                                                                                                                                    
cGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTg                                                                                                                                                     >  2:364025/1‑150 (MQ=255)
        aaCGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                             >  2:1133256/1‑150 (MQ=255)
        aaCGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                             <  1:701005/150‑1 (MQ=255)
             cgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                        <  2:988699/150‑1 (MQ=255)
             cgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                        <  1:1136910/150‑1 (MQ=255)
               acggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                                                      <  1:980223/150‑1 (MQ=255)
                 ggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTCCCCAAACGACACCCGCACATTCGACCTGTACGACTATGTTTCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGt                                                                                                                                    >  1:610948/1‑150 (MQ=255)
                     tGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAg                                                                                                                                <  2:1377764/150‑1 (MQ=255)
                            cGAGAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCgt                                                                                                                         <  1:46453/150‑1 (MQ=255)
                            cGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCgt                                                                                                                         <  1:1217923/150‑1 (MQ=255)
                                agtacggcggCGCGGGTGCTGTCCTACCTGGACGACCAGGCCAACTCCGACCTGTTAGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGCCGTCTTCGCGGGCTTGCCCGCGAAGTGGCCTCTACAGGTTTTGTGTAc                                                                                                                     <  2:100669/146‑1 (MQ=255)
                                aGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTAc                                                                                                                     <  1:1133256/150‑1 (MQ=255)
                                   ccGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                                                      <  2:1240454/130‑1 (MQ=255)
                                   ccGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                                                      >  1:1240454/1‑130 (MQ=255)
                                         gcgcGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGa                                                                                                            <  2:916673/150‑1 (MQ=255)
                                         gcgcGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGa                                                                                                            >  2:603197/1‑150 (MQ=255)
                                                tGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCg                                                                                                     >  2:469062/1‑150 (MQ=255)
                                                 gCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCgg                                                                                                    >  2:1197572/1‑150 (MQ=255)
                                                   tGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTg                                                                                                  <  2:947228/150‑1 (MQ=255)
                                                   tGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTg                                                                                                  >  1:947228/1‑150 (MQ=255)
                                                   tGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTg                                                                                                  >  1:410925/1‑150 (MQ=255)
                                                      ccTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGAc                                                                                               <  1:1206368/150‑1 (MQ=255)
                                                       cTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACg                                                                                              >  1:1418296/1‑150 (MQ=255)
                                                       cTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCAGGCACTGTTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACg                                                                                              >  2:1206368/1‑150 (MQ=255)
                                                         aCCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGaa                                                                                            <  1:537986/150‑1 (MQ=255)
                                                           cTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATa                                                                                          >  2:1205264/1‑150 (MQ=255)
                                                            tGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGt                                                                                                          >  2:323145/1‑133 (MQ=255)
                                                            tGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGt                                                                                                          <  1:323145/133‑1 (MQ=255)
                                                               acgacCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAg                                                                                      >  1:252761/1‑150 (MQ=255)
                                                                  acCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                           <  2:870051/94‑1 (MQ=255)
                                                                  acCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                           >  1:870051/1‑94 (MQ=255)
                                                                  acCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGGGTTCCCGGTGGAc                                                                                               >  2:246820/1‑138 (MQ=255)
                                                                  acCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGGGTTCCCGGTGGAc                                                                                               <  1:246820/138‑1 (MQ=255)
                                                                   cccGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCaaa                                                                                  <  1:627682/150‑1 (MQ=255)
                                                                    ccGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCACGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAc                                                                                 >  2:51064/1‑150 (MQ=255)
                                                                       gCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCAt                                                                              <  2:410925/150‑1 (MQ=255)
                                                                        ccAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATc                                                                             <  2:252761/150‑1 (MQ=255)
                                                                         cAACTCCGACCTGTTCGTCGAAGTGGCGGGGTCCCAGTGCTCGGGCAAATGACCGGGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATcc                                                                            <  1:364025/150‑1 (MQ=255)
                                                                         cAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATAcc                                                                                        >  2:1303395/1‑138 (MQ=255)
                                                                         cAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATAcc                                                                                        <  1:1303395/138‑1 (MQ=255)
                                                                          aaCTCCAACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCcg                                                                           >  2:858418/1‑150 (MQ=255)
                                                                            cTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAc                                                                                 >  1:1029296/1‑142 (MQ=255)
                                                                            cTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAc                                                                                 <  2:1029296/142‑1 (MQ=255)
                                                                             tCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCgcgc                                                                        >  1:1370477/1‑150 (MQ=255)
                                                                               cGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                             <  1:908530/79‑1 (MQ=39)
                                                                               cGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                             >  2:908530/1‑79 (MQ=39)
                                                                                  ccTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCt                                                                     >  1:1243817/1‑148 (MQ=255)
                                                                                  ccTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCt                                                                     <  2:1243817/148‑1 (MQ=255)
                                                                                   cTGTGCGGCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCGGCAAAGAGGCCGGTATAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTc                                                                  <  2:660199/150‑1 (MQ=255)
                                                                                     gTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGa                                                                                                            >  1:65755/1‑106 (MQ=255)
                                                                                     gTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGa                                                                                                            <  2:65755/106‑1 (MQ=255)
                                                                                          tcgAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGGTGGGTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGc                                                           >  1:116236/1‑150 (MQ=255)
                                                                                            gAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                                                      >  1:1520153/1‑73 (MQ=255)
                                                                                            gAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                                                      <  2:1520153/73‑1 (MQ=255)
                                                                                            gAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAt                                                         <  2:1186078/150‑1 (MQ=255)
                                                                                            gAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAt                                                         >  2:576291/1‑150 (MQ=255)
                                                                                                 ggCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGGCTTCGCGGGCTTGCCGGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAtcgtcg                                                    <  2:478039/150‑1 (MQ=255)
                                                                                                  gcgcgcTTCCAGCGCTCGGGCATATGAGCGTGTTGGCGTCTTCGCAGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTcgc                                                   <  2:640995/150‑1 (MQ=255)
                                                                                                  gcgcGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTc                                                                                                                           >  2:232493/1‑78 (MQ=255)
                                                                                                  gcgcGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTc                                                                                                                           <  1:232493/78‑1 (MQ=255)
                                                                                                    gcgGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCgg                                                 >  1:382368/1‑150 (MQ=255)
                                                                                                        ttCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCgg                                                                                                    >  1:490620/1‑95 (MQ=255)
                                                                                                        ttCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCgg                                                                                                    <  2:490620/95‑1 (MQ=255)
                                                                                                          ccAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATCCCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGa                                           <  1:788208/150‑1 (MQ=255)
                                                                                                           cAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGTGTGTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGAc                                          >  1:515909/1‑150 (MQ=255)
                                                                                                           cAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGAc                                          <  2:1594251/150‑1 (MQ=255)
                                                                                                                         aTGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTcc                            >  1:1585646/1‑150 (MQ=255)
                                                                                                                             ccGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGa                        >  2:1254216/1‑150 (MQ=255)
                                                                                                                                  ttGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCcgc                                                                          >  2:337473/1‑95 (MQ=255)
                                                                                                                                  ttGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCcgc                                                                          <  1:337473/95‑1 (MQ=255)
                                                                                                                                    ggCCTCTTCGCGGGCAAGCCCGCGAAGACGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTg                 >  1:220131/1‑150 (MQ=255)
                                                                                                                                      cGTGTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGaa               <  2:551785/150‑1 (MQ=255)
                                                                                                                                        tCTTCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGt             >  1:1432310/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCCTGTCCACGCCACTACTTGTCAAAAAATCTATACCAGACAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTgg           <  1:603197/150‑1 (MQ=255)
                                                                                                                                          ttCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCg                                                                                                     <  1:1582641/60‑1 (MQ=25)
                                                                                                                                          ttCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCg                                                                                                     >  2:1582641/1‑60 (MQ=25)
                                                                                                                                          ttCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGGCGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTg            >  2:1492524/1‑149 (MQ=255)
                                                                                                                                          ttCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGGCGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTg            <  1:1492524/149‑1 (MQ=255)
                                                                                                                                          ttCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTgg           >  2:1614430/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTgg           >  1:1536818/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGGGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTgg           >  2:848148/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGCAAGCCCGCGAAGACGCCGGTACAGGTTTCGTGTACTAGAAAGGCGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTgg           >  1:293435/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGCAAGCCCGCGAAGACGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTgg           <  2:220131/150‑1 (MQ=255)
                                                                                                                                             gcgGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgac        >  2:1490434/1‑150 (MQ=255)
                                                                                                                                               gggCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgacga      <  1:1197572/150‑1 (MQ=255)
                                                                                                                                               gggCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgacga      >  2:1166668/1‑150 (MQ=255)
                                                                                                                                               gggCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgacga      <  1:576291/150‑1 (MQ=255)
                                                                                                                                                ggCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgacgac     <  2:382368/150‑1 (MQ=255)
                                                                                                                                                   ttGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGACGacaac  <  1:469062/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGACGACAAC  >  NC_002947/1704531‑1704827

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