Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 220,562 C→G 7.0% intergenic (+19/‑109) PP_0168 → / → PP_5431 putative surface adhesion protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947220,5620CG 7.0% 218.2 / 6.3 71intergenic (+19/‑109)PP_0168/PP_5431putative surface adhesion protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base G (2/3);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

CGAAGGTGGTATCGCCAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTACGTGCAACGCAA  >  NC_002947/220417‑220708
                                                                                                                                                 |                                                                                                                                                  
cGAAGGTGGTATCGCCAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGCCCCTCTGACAGACCACCGGGGCCGCACAGCg                                                                                                                                                <  1:398630/150‑1 (MQ=255)
        gtATCGCCAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGc                                                                                                                                        >  2:164214/1‑150 (MQ=255)
        gtATCGCCAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGACCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTCACAGACCACCGGGGCCGCACAGCGCCCCCGGc                                                                                                                                        >  2:1102559/1‑150 (MQ=255)
               cAACGCCGACGTGACCATCAAGCTGGAAGGCGTCCACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTcc                                                                                                                                 >  2:588110/1‑150 (MQ=255)
                                tCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAgc                                                                                                                >  2:191590/1‑150 (MQ=255)
                                     cTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATa                                                                                                           <  1:164214/150‑1 (MQ=255)
                                       ggAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTc                                                                                                                                  >  2:396704/1‑125 (MQ=255)
                                       ggAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTc                                                                                                                                  <  1:396704/125‑1 (MQ=255)
                                            gCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGgc                                                                                                    <  2:744334/150‑1 (MQ=255)
                                                 aaCTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTg                                                                                                      <  2:1376967/143‑1 (MQ=255)
                                                 aaCTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTg                                                                                                      >  1:1376967/1‑143 (MQ=255)
                                                 aaCTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCgg                                                                                               >  1:1374193/1‑150 (MQ=255)
                                                 aaCTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCgg                                                                                               <  1:395345/150‑1 (MQ=255)
                                                    tGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTc                                                                                                                                  >  2:1327367/1‑112 (MQ=255)
                                                    tGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTc                                                                                                                                  <  1:1327367/112‑1 (MQ=255)
                                                       tCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTgc                                                                                         >  1:423324/1‑150 (MQ=255)
                                                         tAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTgcgc                                                                                       >  2:1313133/1‑150 (MQ=255)
                                                               gctgATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGACCACCACGCATACTCTGGCGCCGGGCCTGCGCGCCCgg                                                                                 >  2:113946/4‑150 (MQ=255)
                                                                    tCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATc                                                                            >  2:1088213/1‑150 (MQ=255)
                                                                      aaCTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCTGTGCGGCCCCggt                                                                                                                                        <  1:1204575/88‑2 (MQ=38)
                                                                      aaCTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCTGTGCGGCCCCggt                                                                                                                                        >  2:1204575/1‑87 (MQ=38)
                                                                       aCTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCttt                                                                         >  1:1285080/1‑150 (MQ=255)
                                                                       aCTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCttt                                                                         >  1:732791/1‑150 (MQ=255)
                                                                       aCTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCttt                                                                         >  2:92364/1‑150 (MQ=255)
                                                                         tCGTTGATCAGCGGCGCCGATCCGACCCAACCCCTTCGGAACAGTTACAGCTGACAGACCACCGGGGCCGCACACCTGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGc                                                                       <  1:672719/150‑1 (MQ=255)
                                                                         tCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACCCACCACCGGGGCCGCCCAGCGGCCCCGGCTATTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGc                                                                       >  1:532924/1‑150 (MQ=255)
                                                                              gATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCgg                                                                                               <  2:909496/121‑1 (MQ=255)
                                                                              gATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCgg                                                                                               >  1:909496/1‑121 (MQ=255)
                                                                                tCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACa                                                                                                                  <  2:844648/100‑1 (MQ=255)
                                                                                tCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACa                                                                                                                  >  1:844648/1‑100 (MQ=255)
                                                                                tCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGc                                                                                             >  1:408942/1‑121 (MQ=255)
                                                                                tCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGc                                                                                             <  2:408942/121‑1 (MQ=255)
                                                                                     ggcgCCGATCCGACCATCATCATCCACACCCAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGTGATCTTTGCCTATGCTGCTGt                                                           <  1:1102559/150‑1 (MQ=255)
                                                                                     ggcgCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGt                                                           <  1:495596/150‑1 (MQ=255)
                                                                                      gcgcCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCGCCGGGGCCGCTGTGCGGCCCCggtg                                                                                                                                       <  2:411473/73‑3 (MQ=25)
                                                                                      gcgcCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCGCCGGGGCCGCTGTGCGGCCCCggtg                                                                                                                                       >  1:411473/1‑71 (MQ=25)
                                                                                        gcCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTa                                                        <  1:113946/150‑1 (MQ=255)
                                                                                         ccGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATATTTGCCTATGCTGCTGTCTAc                                                       <  2:423324/150‑1 (MQ=255)
                                                                                         ccGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCCCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTAc                                                       <  2:1353021/150‑1 (MQ=255)
                                                                                           gATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACca                                                     >  1:643975/1‑150 (MQ=255)
                                                                                             tCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTg                                                                                                                               >  2:1312273/1‑74 (MQ=255)
                                                                                             tCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTg                                                                                                                               <  1:1312273/74‑1 (MQ=255)
                                                                                                gACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAg                                                >  1:1351720/1‑150 (MQ=255)
                                                                                                caCCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAg                                                <  1:468010/149‑1 (MQ=255)
                                                                                                  ccATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGa                                              >  1:915524/1‑150 (MQ=255)
                                                                                                      catcatcCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAaaca                                          >  2:1345022/1‑150 (MQ=255)
                                                                                                       atcatcCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAaacaa                                         <  2:1351720/150‑1 (MQ=255)
                                                                                                            ccACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAaacaa                                         <  2:1212069/145‑1 (MQ=255)
                                                                                                            ccACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAaacaa                                         >  1:1212069/1‑145 (MQ=255)
                                                                                                             cacaACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATAATCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGAc                                   >  2:1294890/1‑150 (MQ=255)
                                                                                                               caacaaGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGa                                 <  2:134433/150‑1 (MQ=255)
                                                                                                                aacaaGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCg                                                                                  <  2:1339628/100‑1 (MQ=255)
                                                                                                                aacaaGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCg                                                                                  >  1:1339628/1‑100 (MQ=255)
                                                                                                                 acaaGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAgg                               <  1:588110/150‑1 (MQ=255)
                                                                                                                   aaGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGa                             <  2:532924/150‑1 (MQ=255)
                                                                                                                     ggACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACa                           >  1:923579/1‑150 (MQ=255)
                                                                                                                         aGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGc                                                                       >  1:1036080/1‑102 (MQ=255)
                                                                                                                         aGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGc                                                                       <  2:1036080/102‑1 (MQ=255)
                                                                                                                         aGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGc                       <  1:92364/150‑1 (MQ=255)
                                                                                                                           cTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTAccc                                                                                                                         <  2:799548/50‑1 (MQ=255)
                                                                                                                           cTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTAccc                                                                                                                         >  1:799548/1‑50 (MQ=255)
                                                                                                                                 gACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCt               <  1:1313133/150‑1 (MQ=255)
                                                                                                                                 gACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCGGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCt               >  2:378494/1‑150 (MQ=255)
                                                                                                                                  accaccGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTa              <  1:191590/150‑1 (MQ=255)
                                                                                                                                   ccaccGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTAc             <  2:643975/150‑1 (MQ=255)
                                                                                                                                   ccaccGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTAc             <  2:923579/150‑1 (MQ=255)
                                                                                                                                      ccGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTACGTg          >  1:2657/1‑150 (MQ=255)
                                                                                                                                        ggggCCGCTGTGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTgc                                                                                         <  2:596117/69‑1 (MQ=38)
                                                                                                                                            ccGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTACGTgcaacgc    >  2:607067/1‑150 (MQ=255)
                                                                                                                                              gCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTACGTgcaacgcaa  >  2:843794/1‑150 (MQ=255)
                                                                                                                                              gCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTACGTgcaacgcaa  <  2:1285080/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGAAGGTGGTATCGCCAACGCCGACGTGACCATCAAGCTGGAAGGCGTCAACTGGTCTAACCAGACGATCAACTCGTTGATCAGCGGCGCCGATCCGACCATCATCATCCACAACAAGGACAGCTGACAGACCACCGGGGCCGCACAGCGGCCCCGGCTCTTTCCTGCTACCCGCCAACAGCGCATACTCTGGCGCCGGGCCTGCGCGCCCGGCGATCTTTGCCTATGCTGCTGTCTACCAACAAGGAAACAACGTGACGAGGGACATCGCCATGTTCTACGTGCAACGCAA  >  NC_002947/220417‑220708

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