Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,169,846 T→A 6.1% S169C (AGC→TGC)  PP_5408 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,169,8460TA 6.1% 239.4 / 6.4 82S169C (AGC→TGC) PP_5408hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/42);  new base A (3/2);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.60e-01

GTCTATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGACGTCGCCAATAAGAGTAAC  >  NC_002947/6169700‑6169995
                                                                                                                                                  |                                                                                                                                                     
gTCTATGTTTTTTTTTTTTCTTTAGAGTACACGTAGAACTTATCAAGGTCCAAGAGGTGCTGAAAGAGGGAGAACCGTCCGCAGAATTGTTCCAGCCGCAATTCCAATTGGTTCCTGTTTGCCGTAGGCCTTGCGCTGCCTATGCTgcgc                                                                                                                                                   >  2:1834631/1‑150 (MQ=255)
 tCTATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTgcgc                                                                                                                                                   >  2:2216135/1‑150 (MQ=255)
 tCTATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTgcgc                                                                                                                                                   >  1:1189816/1‑150 (MQ=255)
   tATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGccc                                                                                                                                                 >  1:274902/1‑150 (MQ=255)
   tATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGccc                                                                                                                                                 >  2:1047906/1‑150 (MQ=255)
    attttttttAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCg                                                                                                                                                <  1:645500/147‑1 (MQ=255)
       ttttttAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCaa                                                                                                                                             <  1:662622/150‑1 (MQ=255)
       ttttttAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCaa                                                                                                                                             >  2:77755/1‑150 (MQ=255)
         ttttAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACa                                                                                                                                           <  2:1344718/150‑1 (MQ=255)
              gTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGc                                                                                                                                      <  1:1547693/150‑1 (MQ=255)
              gTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGc                                                                                                                                      <  2:1859717/150‑1 (MQ=255)
                  gCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCca                                                                                                                                  <  2:1333906/150‑1 (MQ=255)
                    tttCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCt                                                                                                                                                       <  2:1846621/127‑1 (MQ=255)
                    tttCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCt                                                                                                                                                       >  1:1846621/1‑127 (MQ=255)
                    tttCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCcaca                                                                                                                                <  1:77755/150‑1 (MQ=255)
                    tttCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCcaca                                                                                                                                >  2:1313603/1‑150 (MQ=255)
                    tttCGAGCACCCGTAGAACCGTTTCAAGTGAAAGGTGGTGATCAAGGAGGCAGATACTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTAACAGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCcaca                                                                                                                                >  2:535315/1‑150 (MQ=255)
                                 tAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACaga                                                                                                                   <  1:1047906/150‑1 (MQ=255)
                                    aaCCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTg                                                                                                                >  2:1522143/1‑150 (MQ=255)
                                    aaCCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTg                                                                                                                <  2:2374670/150‑1 (MQ=255)
                                    aaCCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTg                                                                                                                <  2:651009/150‑1 (MQ=255)
                                     aCCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTgg                                                                                                               <  1:1924706/150‑1 (MQ=255)
                                      ccGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGc                                                                                                              >  2:1302982/1‑150 (MQ=255)
                                         tttCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCAt                                                                                                           >  2:2288301/1‑150 (MQ=255)
                                             aaGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATcgcg                                                                                                       >  1:1041958/1‑150 (MQ=255)
                                                        ggtgCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCAGCCTAGGCAgcgcaa                                                                                                                                                 <  2:325682/97‑3 (MQ=255)
                                                        ggtgCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCAGCCTAGGCAgcgcaa                                                                                                                                                 >  1:325682/1‑95 (MQ=255)
                                                          tgCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGgctgc                                                                                          <  2:886622/150‑1 (MQ=255)
                                                             tCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTca                                                                                       <  2:1615428/150‑1 (MQ=255)
                                                               aaGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTcaca                                                                                     <  2:1041958/150‑1 (MQ=255)
                                                               aaGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATCGACCCGGGAACCCGGTAGTTTTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTcaca                                                                                     <  2:63748/150‑1 (MQ=255)
                                                                   aggCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAAc                                                                                 <  2:335472/150‑1 (MQ=255)
                                                                    ggCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACt                                                                                <  1:1302982/150‑1 (MQ=255)
                                                                     gCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCt                                                                                                                          <  2:1674657/107‑1 (MQ=255)
                                                                     gCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCt                                                                                                                          >  1:1674657/1‑107 (MQ=255)
                                                                     gCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTg                                                                               >  1:1984129/1‑150 (MQ=255)
                                                                     gCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTg                                                                               >  1:358026/1‑150 (MQ=255)
                                                                      cAGATTCTTCCGTAGAATTGTGCCAGCCGCAATGCTAATTGTTCCCTGTTGGCTGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGTGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTgg                                                                              <  1:535315/150‑1 (MQ=255)
                                                                           ccTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTa                                                                         <  1:1110950/150‑1 (MQ=255)
                                                                            ctgccGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCACGAAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTAc                                                                        <  1:1741728/147‑1 (MQ=255)
                                                                             ttCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACg                                                                       >  2:772905/1‑150 (MQ=255)
                                                                              tCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACga                                                                      <  1:2216135/150‑1 (MQ=255)
                                                                               ccGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACgag                                                                     <  2:866736/150‑1 (MQ=255)
                                                                                cGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACgaga                                                                    <  2:1189816/150‑1 (MQ=255)
                                                                                      atcgtgCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCAGGACAAAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAAt                                                              <  1:1313603/147‑1 (MQ=255)
                                                                                      aTTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGAATCGTCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAAt                                                              <  1:1834631/150‑1 (MQ=255)
                                                                                            ccAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGt                                                        >  1:392001/1‑150 (MQ=255)
                                                                                            ccAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGt                                                        >  2:2205372/1‑150 (MQ=255)
                                                                                             cAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCTCGCAACATCCGCTCCACATTACCTGGACAGAGTg                                                                                                                <  2:871542/93‑1 (MQ=255)
                                                                                             cAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTg                                                                                                                >  1:871542/1‑93 (MQ=255)
                                                                                             cAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGtt                                                       >  2:2229879/1‑150 (MQ=255)
                                                                                              aGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTtg                                                      <  2:2054227/150‑1 (MQ=255)
                                                                                               gCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTtgt                                                     >  1:238815/1‑150 (MQ=255)
                                                                                                    aaTGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTc                                                >  2:1164424/1‑150 (MQ=255)
                                                                                                        ccAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCgcgc                                            >  1:227528/1‑150 (MQ=255)
                                                                                                          aaTGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGcac                                          >  1:1227035/1‑150 (MQ=255)
                                                                                                           aTGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGcaca                                         >  1:2328186/1‑150 (MQ=255)
                                                                                                           aTGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGACCGCAACATCCGCTCCACATAACCTGGTCGCAGTGGCCACTGTGATAAAATGGCTGATGACGGAACGGCCATTACGAGAGCAAATCGACGTTGTGCTTCGCGcaca                                         <  2:274902/150‑1 (MQ=255)
                                                                                                                acccTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATcgc                                                                                                        >  2:1044415/2‑82 (MQ=255)
                                                                                                                acccTGTTGGCCGTAGGCCTTGCGCTGCATAGGCTTCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATcgc                                                                                                        <  1:1044415/81‑1 (MQ=255)
                                                                                                                  ccTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTgg                                                                              <  1:309636/106‑1 (MQ=255)
                                                                                                                  ccTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTgg                                                                              >  2:309636/1‑106 (MQ=255)
                                                                                                                      ttGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGgct                                                                                            >  1:1822706/1‑88 (MQ=255)
                                                                                                                      ttGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGgct                                                                                            <  2:1822706/88‑1 (MQ=255)
                                                                                                                      ttGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTAACTACAACGAGTCGCAATGGCTAGATCATTGGTTCGGACAGTTGTGGGAGTACGAGTGTAAATTGACGGTGTGCTTCGCGCACACAGTTACAtgt                              <  2:275084/150‑1 (MQ=255)
                                                                                                                            gTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGa                        <  1:2205372/150‑1 (MQ=255)
                                                                                                                              aGGCCTTGCGCTGACTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGGGGGTGAAg                      <  2:437267/150‑1 (MQ=255)
                                                                                                                               ttgcgggcgcAGCCTAGGCAGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGa                     >  1:2408473/7‑150 (MQ=255)
                                                                                                                                gCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGAc                    >  2:2149871/1‑150 (MQ=255)
                                                                                                                                 gcgggcgcAGCCTAGGCAGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTTCTCACAGAACTgg                                                                              <  2:593816/87‑1 (MQ=255)
                                                                                                                                 gcgggcgcAGCCTAGGCAGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTgg                                                                              >  1:593816/5‑91 (MQ=255)
                                                                                                                                  cTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTAc                                                                                                                            >  2:1873579/1‑44 (MQ=255)
                                                                                                                                  cTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTAc                                                                                                                            <  1:1873579/44‑1 (MQ=255)
                                                                                                                                  cTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGACGt                  >  2:1082684/1‑150 (MQ=255)
                                                                                                                                   ttGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGACGTc                 <  1:2149871/150‑1 (MQ=255)
                                                                                                                                           ccTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCgag                                                                                                     >  1:880927/1‑58 (MQ=255)
                                                                                                                                           ccTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCgag                                                                                                     <  2:880927/58‑1 (MQ=255)
                                                                                                                                           ccTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGcaca                                         <  2:169029/118‑1 (MQ=255)
                                                                                                                                           ccTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGcaca                                         >  1:169029/1‑118 (MQ=255)
                                                                                                                                             tAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCGGTACAAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGACGTCGCCAATAaga       >  1:1620497/1‑150 (MQ=255)
                                                                                                                                             tAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGACGTCGCCAATAaga       <  2:238815/150‑1 (MQ=255)
                                                                                                                                              aGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGACGTCGCCAATAagag      <  2:358026/150‑1 (MQ=255)
                                                                                                                                                gCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTAc                                                                        >  2:679147/1‑82 (MQ=255)
                                                                                                                                                gCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTAc                                                                        <  1:679147/82‑1 (MQ=255)
                                                                                                                                                  tGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGACGTCGCCAATAAGAGTAAc  >  1:1733101/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTCTATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGACGTCGCCAATAAGAGTAAC  >  NC_002947/6169700‑6169995

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: