Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 337,831 Δ1 bp 5.7% intergenic (‑78/+51) PP_0279 ← / ← PP_0280 hypothetical protein/amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947337,8310T. 5.7% 216.2 / 6.2 70intergenic (‑78/+51)PP_0279/PP_0280hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (39/26);  new base . (2/2);  total (42/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCG  >  NC_002947/337687‑337978
                                                                                                                                                 |                                                                                                                                                   
acctggtgatctcgttcatgctggtgcaaggcttcaagctgctggagcgctggttgcgcgtggatgcctgccaaggccgctgagacctcactgactgCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGgcgc                                                                                                                                                >  1:1911131‑M1/98‑150 (MQ=255)
 gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCt                                                                                                                                               >  2:805684/1‑150 (MQ=255)
   tggtgatctcgttcatgctggtgcaaggcttcaagctgctggagcgctggttgcgcgtggatgcctgccaaggccgctgagacctcactgactgCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTac                                                                                                                                             >  2:1917866‑M1/95‑150 (MQ=255)
   tggtgatctcgttcatgctggtgcaaggcttcaagctgctggagcgctggttgcgcgtggatgcctgccaaggccgctgagacctcactgactgCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTac                                                                                                                                             >  2:760974‑M1/95‑150 (MQ=255)
   tCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTac                                                                                                                                             <  2:1746191/150‑1 (MQ=255)
    cAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACa                                                                                                                                            >  2:422154‑M2/1‑149 (MQ=255)
    cAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACa                                                                                                                                            >  1:1769486‑M2/1‑149 (MQ=255)
     gtgatcttgttcatgttggtgaaacgtttcaagctgcttgagcgctggttgcgcgtggatgcctgccaaggccgctgagacctcactgactgCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTacac                                                                                                                                           >  2:1852716‑M1/93‑150 (MQ=255)
     gtgatctcgttcatgctggtgcaaggcttcaagctgctggagcgctggttgcgcgtggatgcctgccaaggccgctgagacctcactgactgCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTacac                                                                                                                                           >  1:255428‑M1/93‑150 (MQ=255)
     aGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGgcgc                                                                                                                                                >  2:530291/1‑145 (MQ=255)
     aGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGATAGCGCCACATTCAAGGTTGgcgc                                                                                                                                                <  1:530291/145‑1 (MQ=255)
       gatctcgttcatgctggtgcaaggcttcaagctgctggagcgctggttgcgcgtggatgcctgccaaggccgctgagacctcactgactgCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCt                                                                                                                                         >  2:798020‑M1/91‑150 (MQ=255)
         aCTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACacctgt                                                                                                                                       >  1:1057667‑M2/1‑144 (MQ=255)
          cTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACacctgtg                                                                                                                                      >  2:53075‑M2/1‑143 (MQ=255)
             gTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACacctgtggga                                                                                                                                   >  1:1208361‑M2/1‑140 (MQ=255)
               cGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACTCCCGCTACCAACGGTATAGCGCCACATTCAAGGTTGGCGCTACacctgtgggagc                                                                                                                                 >  1:805684‑M2/1‑138 (MQ=255)
                 aGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACacctgtgggagcgg                                                                                                                               >  1:1855747‑M2/1‑136 (MQ=255)
                 aGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTAAacctgtgggcgcgg                                                                                                                               >  2:1763022‑M2/1‑136 (MQ=255)
                   ggCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGTATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACacctgtgggagcgggc                                                                                                                             >  1:253845‑M2/1‑134 (MQ=255)
                         gAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACacctgtgggagcgggcgtgccc                                                                                                                       >  1:1818222‑M2/1‑128 (MQ=255)
                                    ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATacctgtgg                                                                                                                                     <  2:127110‑M2/125‑9 (MQ=255)
                                    ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATacctgtgg                                                                                                                                     >  1:127110‑M2/1‑117 (MQ=255)
                                                             gggCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACac                                                                                                                                           >  2:687367‑M2/1‑92 (MQ=255)
                                                             gggCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACac                                                                                                                                           <  1:687367‑M2/94‑3 (MQ=255)
                                                                   ggCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTa                                                                                                                                              <  1:994754/85‑1 (MQ=255)
                                                                   ggCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTa                                                                                                                                              >  2:994754/1‑85 (MQ=255)
                                                                   ctgccaaggccgctgagacctcactgactgCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCt                                                                                                                                         <  1:563331‑M1/60‑1 (MQ=255)
                                                                   ctgccaaggccgctgagacctcactgactgCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCt                                                                                                                                         >  2:563331‑M1/31‑90 (MQ=255)
                                                                    tgccaaggccgctgagacctcactgactgCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTacac                                                                                                                                           <  1:1815972‑M1/58‑1 (MQ=255)
                                                                    tgccaaggccgctgagacctcactgactgCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTacac                                                                                                                                           >  2:1815972‑M1/30‑87 (MQ=255)
                                                                    tgccaaggccgctgagacctcactgactgCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGt                                                                                                                           <  1:981977‑M1/74‑1 (MQ=255)
                                                                    tgccaaggccgctgagacctcactgactgCTTCGCGGGCACGCCCGCTCCCACAGCTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGt                                                                                                                           >  2:981977‑M1/30‑103 (MQ=255)
                                                                                       tcactgactgCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTacac                                                                                                                                           <  2:1066769‑M1/58‑1 (MQ=255)
                                                                                       tcactgactgCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTacac                                                                                                                                           >  1:1066769‑M1/11‑68 (MQ=255)
                                                                                             gCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGgc                                                                                                                                                  <  1:1098488/55‑1 (MQ=255)
                                                                                             gCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGgc                                                                                                                                                  >  2:1098488/1‑55 (MQ=255)
                                                                                                      gGGGCACGCCCGCTCGCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTg                                          >  1:192815‑M1/2‑150 (MQ=255)
                                                                                                           ACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCa                                     <  2:1067795‑M1/150‑1 (MQ=255)
                                                                                                                   TCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                             <  2:1226273‑M1/150‑1 (MQ=255)
                                                                                                                     CCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCa                                     >  1:807741‑M1/1‑140 (MQ=255)
                                                                                                                     CCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCa                                     <  2:807741‑M1/140‑1 (MQ=255)
                                                                                                                     CCACAGGTATAGCGCCACATTCAAGGTTGGCGTTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                           <  2:1070244‑M1/150‑1 (MQ=255)
                                                                                                                      cacaGGTATAGCGCCACATTCCAGGTTGGCCCTACacctgtgcgagcgggcgtgcccgcgaagaggccctcaaaaccaccacacatccttgccatgcccaaccacctccccagcccccacctaagcccccccttcccggccctcgaccac                          <  2:501345‑M2/150‑116 (MQ=255)
                                                                                                                      cacaGGTATAGCGCCACATTCAAGGTTGGCGCTACacctgtgggagcgggcgtgcccgcgaagaggccctcaaaaaaacggcacatccttgccatgcccaaccacataaacagccaccacctaagcacccgcttccgggccctcgaccag                          <  1:329154‑M2/150‑116 (MQ=255)
                                                                                                                      CACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgag                          <  2:1592856‑M1/150‑1 (MQ=255)
                                                                                                                       acaGGTATAGCGCCACATTCAAGGTTGGCGCTACacctgtgggagcgggcgtgcccgcgaagaggccctcaaaaccaccacacatccttgccatgcccaaccacctccacagcccccacctaa                                                    <  2:1760404‑M2/123‑90 (MQ=255)
                                                                                                                       acaGGTATAGCGCCACATTCAAGGTTGGCGCTACacctgtgggagcgggcgtgcccgcgaagaggccctcaaaaccaccacacatccttgccatgcccaaccacctccacagcccccacctaa                                                    >  1:1760404‑M2/1‑34 (MQ=255)
                                                                                                                       ACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgaga                         <  1:829652‑M1/150‑1 (MQ=255)
                                                                                                                          GGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTTGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTATCAGCGGCCTTGGCAGGCATCCACTCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCa                      <  2:1654368‑M1/150‑1 (MQ=255)
                                                                                                                           gTATAGCGCCACATTCAAGGTTGGCGCTACacctgtgggagcgggcgtgcccgcgaagaggccctcaaaaccaccacacatccttgccatgcccaaccacctccacagcccccac                                                        <  2:1481080‑M2/115‑86 (MQ=255)
                                                                                                                           gTATAGCGCCACATTCAAGGTTGGCGCTACacctgtgggagcgggcgtgcccgcgaagaggccctcaaaaccaccacacatccttgccatgcccaaccacctccacagcccccac                                                        >  1:1481080‑M2/1‑30 (MQ=255)
                                                                                                                           gTATAGCGCCACATTCAAGGTTGGCGCTACacctgtgggagcgggcgtgcccgcgaagaggccctcaaaaccaccacacatccttg                                                                                     <  2:1195692‑M2/86‑57 (MQ=255)
                                                                                                                           gTATAGCGCCACATTCAAGGTTGGCGCTACacctgtgggagcgggcgtgcccgcgaagaggccctcaaaaccaccacacatccttg                                                                                     >  1:1195692‑M2/1‑30 (MQ=255)
                                                                                                                            TATAGCGCCACCTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTCAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCCGCTTGAAGCCTTGCACCAGCATGAACGAGATGAcc                    <  2:721155‑M1/150‑1 (MQ=255)
                                                                                                                             ataGCGCCACATTCAAGGTTGGCGCTACacctgtgggagcgggcgtgcccgcgaagaggccctcaaaaccacca                                                                                               <  2:380131‑M2/74‑47 (MQ=255)
                                                                                                                             ataGCGCCACATTCAAGGTTGGCGCTACacctgtgggagcgggcgtgcccgcgaagaggccctcaaaaccacca                                                                                               >  1:380131‑M2/1‑28 (MQ=255)
                                                                                                                             agagcGCCACATTCAAGGGGGGCGCTAGacctgggggagcgggcgtgcccgcgaagaggccctcaaaaccaccacacatcctggccatgcccaaccacctccacagcccccacctaagcacccgcttccaggccctcgaccagttcctga                   >  1:501345‑M2/3‑28 (MQ=255)
                                                                                                                              taGCGCCACATTCAAGGTTGGCGCTATacctgtgggagcgggcgtgcccgcgaagaggccctcaaaaccaccacacatccttgccatgcccaaccacctc                                                                    >  1:285004‑M2/1‑27 (MQ=255)
                                                                                                                              taGCGCCACATTCAAGGTTGGCGCTATacctgtgggagcgggcgtgcccgcgaagaggccctcaaaaccaccacacatccttgccatgcccaaccacctc                                                                    <  2:285004‑M2/100‑74 (MQ=255)
                                                                                                                                  gcCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggcag                                                                                   >  1:44561/1‑81 (MQ=255)
                                                                                                                                  gcCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggcag                                                                                   <  2:44561/81‑1 (MQ=255)
                                                                                                                                  gcCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAg              >  1:364407/1‑150 (MQ=255)
                                                                                                                                        aTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGAGAGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAacaa         >  1:602957/1‑149 (MQ=255)
                                                                                                                                         ttCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggcaggc                                                                                 <  2:793306/77‑1 (MQ=255)
                                                                                                                                         ttCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggcaggc                                                                                 >  1:793306/1‑77 (MQ=255)
                                                                                                                                         ttCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACg        >  1:1742211/1‑150 (MQ=255)
                                                                                                                                            aaGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGt                <  2:939609/139‑1 (MQ=255)
                                                                                                                                            aaGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGt                >  1:939609/1‑139 (MQ=255)
                                                                                                                                            aaGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCt     >  1:1181856/1‑150 (MQ=255)
                                                                                                                                               gTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTgcg  >  1:481358/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCG  >  NC_002947/337687‑337978

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: