Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,709,508 T→C 6.4% L15P (CTG→CCG) ‡ PP_2372 → sugar‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,709,5080TC 6.4% 181.4 / 2.4 62L15P (CTG→CCG) ‡PP_2372sugar‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (33/25);  new base C (2/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

CATATGAAAAGAACCAAATGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGACCCTGCAGGTGCAATATAGCGAGCGC  >  NC_002947/2709360‑2709650
                                                                                                                                                    |                                                                                                                                              
cATATGAAAAGAACCAAATGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACtg                                                                                                                                               <  2:942207/150‑1 (MQ=255)
   atGAAAAGAACCAAATGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACtggtg                                                                                                                                            >  1:409687/1‑150 (MQ=255)
    tGAAAAGAACCAAATGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACtggtgg                                                                                                                                           <  2:169950/150‑1 (MQ=255)
      aaaaGAACCAAATGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACCAGTGGCa                                                                                                                                         >  1:349722/1‑150 (MQ=255)
        aaGAACCAAATGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCTCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGTCTTGTCTTGCCACTGGTGGCaaa                                                                                                                                       >  1:453703/1‑150 (MQ=255)
         aGAACCAAATGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACCAGTGGCaaga                                                                                                                                      <  2:349722/150‑3 (MQ=255)
               aaaTGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCaa                                                                                                                                >  2:75723/1‑150 (MQ=255)
               aaaTGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCaa                                                                                                                                >  1:312394/1‑150 (MQ=255)
                 aTGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAAcc                                                                                                                              <  1:205639/150‑1 (MQ=255)
                      cTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCtt                                                                                                                         <  2:263083/150‑1 (MQ=255)
                             agaTTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCggtggt                                                                                                                  <  1:965035/150‑1 (MQ=255)
                                 tGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGttt                                                                                                              <  1:501033/150‑1 (MQ=255)
                                    cTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACCAGTGGCaagacaag                                                                                                                                  <  1:102754/127‑7 (MQ=255)
                                    cTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACCAGTGGCaagacaag                                                                                                                                  >  2:102754/1‑121 (MQ=255)
                                     tGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCa                                                                                                          >  1:967663/1‑150 (MQ=255)
                                      gCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCaa                                                                                                         >  2:1004302/1‑150 (MQ=255)
                                                   tGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATgcgc                                                                                                                                  <  2:523745/112‑1 (MQ=255)
                                                   tGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATgcgc                                                                                                                                  >  1:523745/1‑112 (MQ=255)
                                                   tGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCgg                                                                                                                            >  2:642411/1‑118 (MQ=255)
                                                   tGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCgg                                                                                                                            <  1:642411/118‑1 (MQ=255)
                                                   tGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTcc                                                                                            >  1:1343309/1‑150 (MQ=255)
                                                        aTCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCACATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCga                                                                                       <  1:983690/150‑1 (MQ=255)
                                                        aTCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCga                                                                                       <  1:345354/150‑1 (MQ=255)
                                                         tCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCCTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCCCGCGGGCGCCGTCCACCgag                                                                                      >  2:1105769/1‑150 (MQ=255)
                                                           gCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGAc                                                                                    >  2:709571/1‑150 (MQ=255)
                                                            ccAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCggtg                                                                                                                    <  2:463040/117‑1 (MQ=255)
                                                            ccAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCggtg                                                                                                                    >  1:463040/1‑117 (MQ=255)
                                                               ggCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACACGACGTTGTCCACCGAGACATTc                                                                                >  1:783858/1‑150 (MQ=255)
                                                                             ggTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATAc                                                                  <  1:75723/150‑1 (MQ=255)
                                                                             ggTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATAc                                                                  >  1:777851/1‑150 (MQ=255)
                                                                             ggTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATAc                                                                  >  1:1341282/1‑150 (MQ=255)
                                                                                 tttGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGc                                                              >  2:1109235/1‑150 (MQ=255)
                                                                                 tttGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGc                                                              >  2:33824/1‑150 (MQ=255)
                                                                                     cTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCa                                                          >  2:428480/1‑150 (MQ=255)
                                                                                     cTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCa                                                          >  2:777933/1‑150 (MQ=255)
                                                                                      ttCGCAGAGATGTTACCCCATGCTCCAGTCCGTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCAt                                                         >  1:750246/1‑150 (MQ=255)
                                                                                      ttCGCAGAGATGTTACCCCATGCCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATg                                                        >  1:829605/1‑150 (MQ=255)
                                                                                         gCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCa                                                      >  1:112238/1‑150 (MQ=255)
                                                                                          cAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCtata                                                                   <  1:1290695/136‑1 (MQ=255)
                                                                                          cAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCtata                                                                   >  2:1290695/1‑136 (MQ=255)
                                                                                            gagaTGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGCCCTCGGCGGTCTCTCTCAACCCGGGGCTGTCCACCGAGACATCCTCGGCCAGCTCTGCCCCCTCCaccgatgc                                                   >  2:190370/1‑143 (MQ=255)
                                                                                                     ccccATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCggtg                                                                                                                    >  1:714234/1‑76 (MQ=255)
                                                                                                     ccccATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCggtg                                                                                                                    <  2:714234/76‑1 (MQ=255)
                                                                                                      cccATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGAc                                         >  1:1077224/1‑150 (MQ=255)
                                                                                                       ccATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATTCAGGCCGCCGCGGACg                                        >  2:1383710/1‑150 (MQ=255)
                                                                                                        cATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGa                                       >  1:218001/1‑150 (MQ=255)
                                                                                                         aTGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAg                                      >  1:389601/1‑150 (MQ=255)
                                                                                                           gCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCt                                    <  1:777933/150‑1 (MQ=255)
                                                                                                            cTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTg                                   <  2:1341282/150‑1 (MQ=255)
                                                                                                               aaGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTggg                                 <  2:115690/149‑1 (MQ=255)
                                                                                                               aaGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTggg                                 >  1:115690/1‑149 (MQ=255)
                                                                                                               aaGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGc                                >  1:930172/1‑150 (MQ=255)
                                                                                                                   cccTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGa                            <  1:33824/150‑1 (MQ=255)
                                                                                                                   cccTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGa                            <  2:777851/150‑1 (MQ=255)
                                                                                                                   cccTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGa                            <  2:1343309/150‑1 (MQ=255)
                                                                                                                        tgCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGACCCTg                       <  1:709571/150‑1 (MQ=255)
                                                                                                                         gCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGACCCTGc                      <  2:829605/150‑1 (MQ=255)
                                                                                                                          ccAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGACCCTGCa                     <  1:1004302/150‑1 (MQ=255)
                                                                                                                           cAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTTGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGACCCTGCAg                    <  2:750246/150‑1 (MQ=255)
                                                                                                                            aaaGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGACCCTGCAgg                   <  2:389601/150‑1 (MQ=255)
                                                                                                                                      ccTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCAt                                                         <  1:816792/102‑1 (MQ=255)
                                                                                                                                      ccTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCAt                                                         >  2:816792/1‑102 (MQ=255)
                                                                                                                                      ccTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGACCCTGCAGGTGCAATATag         >  2:635241/1‑150 (MQ=255)
                                                                                                                                          gTCTTGCCACTGGTGGCAAATGCTCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGACCCTGCAGGTGCAATATagcgag     <  2:453703/150‑1 (MQ=255)
                                                                                                                                             ttGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGACCCTGCAGGTGCAATATAGCGAgcgc  >  1:83350/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CATATGAAAAGAACCAAATGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGACCCTGCAGGTGCAATATAGCGAGCGC  >  NC_002947/2709360‑2709650

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