Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,063,271 C→G 11.3% intergenic (+85/‑122) PP_1838 → / → PP_1839 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,063,2710CG11.3% 191.8 / 13.3 68intergenic (+85/‑122)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/27);  new base G (4/4);  total (37/31)
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.39e-01

TCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCAT  >  NC_002947/2063122‑2063412
                                                                                                                                                     |                                                                                                                                             
tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                                               <  2:2175488/150‑1 (MQ=255)
tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                                               >  2:383820/1‑150 (MQ=255)
  cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCga                                                                                                                                             >  2:1899820‑M2/1‑148 (MQ=255)
   aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGgag                                                                                                                                            <  2:252556‑M2/150‑4 (MQ=255)
      ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGgagcgg                                                                                                                                         >  1:1046565‑M2/1‑144 (MQ=255)
      ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGgagcgg                                                                                                                                         >  1:252556‑M2/1‑144 (MQ=255)
       cgcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCgagcggg                                                                                                                                        <  2:1798747‑M2/150‑8 (MQ=255)
          gACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGTGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCTCAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCgagcgggcaa                                                                                                                                     >  2:253323‑M2/1‑140 (MQ=255)
          gACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCgagcgggcaa                                                                                                                                     >  1:383820‑M2/1‑140 (MQ=255)
                        tgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCgagcgggcaagcccgcgaatgggc                                                                                                                       >  2:1296527‑M2/1‑126 (MQ=255)
                                                aCGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGGCCCGCACAATTTTCCCGGACTGaa                                                                                               <  2:1623445/150‑1 (MQ=255)
                                                     acacCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTa                                                                                          >  1:2074658/1‑150 (MQ=255)
                                                       acCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCgagcgggcaagcccgcgaatg                                                                                                                          >  2:1579753‑M2/1‑95 (MQ=255)
                                                       acCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCgagcgggcaagcccgcgaatg                                                                                                                          <  1:1579753‑M2/116‑22 (MQ=255)
                                                       acCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATc                                                                                        >  1:91764/1‑150 (MQ=255)
                                                       acCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATc                                                                                        <  1:1798747/150‑1 (MQ=255)
                                                                       gCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAg                                                                        <  2:91764/150‑1 (MQ=255)
                                                                                   cccATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTTTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTa                                                                                          >  1:80041/1‑120 (MQ=255)
                                                                                   cccATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTa                                                                                          <  2:80041/120‑1 (MQ=255)
                                                                                    ccATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATc                                                           >  2:2393680/1‑150 (MQ=255)
                                                                                    ccATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATc                                                           >  2:635390/1‑150 (MQ=255)
                                                                                     cATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGGGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATca                                                          <  2:388929/150‑1 (MQ=255)
                                                                                          gcgGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCacacac                                                     >  1:1203296/1‑150 (MQ=255)
                                                                                              gCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTgg                                                 >  1:98843/1‑150 (MQ=255)
                                                                                                ttGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCTGCACCATTTTCCAGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgtg                                               >  1:1475355/1‑150 (MQ=14)
                                                                                                   cccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                            >  1:967916/1‑150 (MQ=255)
                                                                                                    ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcggtgtgccgtgcaccacgccggcggtgttatccagcaggctcttg                                                  >  2:2260733‑M2/1‑50 (MQ=255)
                                                                                                    ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcggtgtgccgtgcaccacgccggcggtgttatccagcaggctcttg                                                  <  2:1277993‑M2/143‑94 (MQ=255)
                                                                                                    ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcggtgtgccgtgcaccacgccggcggtgttatccagcaggctcttg                                                  >  1:1277993‑M2/1‑50 (MQ=255)
                                                                                                    ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcggtgtgccgtgcaccacgccggcggtgttatccagcaggctcttg                                                  <  1:2260733‑M2/143‑94 (MQ=255)
                                                                                                    ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                                                 >  1:636291/1‑128 (MQ=255)
                                                                                                    ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                                                 <  2:636291/128‑1 (MQ=255)
                                                                                                    ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                                           >  1:1443389/1‑150 (MQ=255)
                                                                                                      gCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                                         <  2:1443389/150‑1 (MQ=255)
                                                                                                          ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                                     >  1:80100/1‑150 (MQ=255)
                                                                                                           aacagGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGGgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcggtgtgccgtgcaccacgccggcggtgttatccagcaggctcttggtcgcggtctggat                                    <  1:1047385‑M2/149‑108 (MQ=255)
                                                                                                            acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGc                                                              <  1:1353005/123‑1 (MQ=255)
                                                                                                            acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGc                                                              >  2:1353005/1‑123 (MQ=255)
                                                                                                            acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                                   >  2:441552/1‑150 (MQ=255)
                                                                                                              aGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga                                 <  2:98843/150‑1 (MQ=255)
                                                                                                              aGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga                                 >  1:2447013/1‑150 (MQ=255)
                                                                                                               ggCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGagag                                >  2:1060199/1‑150 (MQ=255)
                                                                                                                gCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcggtgtgccgtgcaccacgccggcggtgttatc                                                               >  1:2145223‑M2/1‑38 (MQ=255)
                                                                                                                gCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcggtgtgccgtgcaccacgccggcggtgttatc                                                               <  2:2145223‑M2/118‑81 (MQ=255)
                                                                                                                gCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg                               <  2:2447013/150‑1 (MQ=255)
                                                                                                                   ttGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTc                            <  2:967916/150‑1 (MQ=255)
                                                                                                                   ctGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAACCCCGCTAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCCAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTc                            >  1:1241055/2‑150 (MQ=255)
                                                                                                                       acaGGGTTCGGACCTTGTGCAGGTGCTGTGCgagcggacaagggcgcaaatgggccgcaacgcggccccggtttacttcggtgtgccgtgcaccacgccggcggtgttatcgagcaggctcttggtcgcggtctggatgtacggttacag                        <  2:587509‑M2/150‑120 (MQ=255)
                                                                                                                       acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt                        <  2:2074658/150‑1 (MQ=255)
                                                                                                                         aGGGTTCGGACCTTGTGCAGGTCCTGTGGgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcggtgtgccgtgcaccacgccggcggtgttatccagcaggctcttggtcgcggtctggatatacgcttccagct                      <  2:151598‑M2/150‑122 (MQ=255)
                                                                                                                         aGGGTTCGGACCTTGTGCAGGTCCTGTGGgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcggtgtgccgtgcaccacgccggcggtgttatccagcaggctcttggtcgcggtctggatatacgcttccagct                      <  1:2165650‑M2/150‑122 (MQ=255)
                                                                                                                         aGGGTTCGGACCTTGTGCAGGTCCTGTGCgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcg                                                                                              <  1:1327299‑M2/78‑50 (MQ=255)
                                                                                                                         aGGGTTCGGACCTTGTGCAGGTCCTGTCCgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcg                                                                                              >  2:1327299‑M2/1‑29 (MQ=255)
                                                                                                                               cGGACCTTGTGCAGGTCCTGTGGgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcggtgtgccgtgcaccacgccggcggtgttatccagcaggctcttggtcgcggtctggatatacgcttccagcttcttga                >  1:151598‑M2/1‑23 (MQ=255)
                                                                                                                                ggACCTTGTGCAGGTCCTGTGCgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcggtgtgccgtgcaccacgccggcggtgttatccagcaggctcttggtcgcggtctggatatacgcttccagcttcttga                <  1:36438‑M2/149‑128 (MQ=255)
                                                                                                                                ggACCTTGTGCAGGTCCTGTGCgagcgggcaagcccgcgaatgggccgcaacgcggccccggtttacttcggtgtgccgtgcaccacgccggcggtgttatccagcaggctcttggtcgcggtctggatatacgcttccagcttcttga                >  2:36438‑M2/1‑22 (MQ=255)
                                                                                                                                ggACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAag               >  1:827693/1‑150 (MQ=255)
                                                                                                                                 gACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTAtt                                                                   >  2:1001914/1‑97 (MQ=255)
                                                                                                                                 gACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTc                                                                  <  2:1086964/98‑1 (MQ=255)
                                                                                                                                 gACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTc                                                                  >  1:1086964/1‑98 (MQ=255)
                                                                                                                                 gACCTTGTGCAGGTCATGGTCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTAtt                                                                   <  1:1001914/97‑1 (MQ=39)
                                                                                                                                   ccTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGt            <  2:2234209/150‑1 (MQ=255)
                                                                                                                                   ccTTGTGCAGGCACTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGt            <  2:2234317/150‑1 (MQ=255)
                                                                                                                                    cTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCaa                     >  1:1579459/1‑140 (MQ=255)
                                                                                                                                    cTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCaa                     <  2:1579459/140‑1 (MQ=255)
                                                                                                                                      tgtgCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGt                                                                       <  2:264943/88‑1 (MQ=255)
                                                                                                                                      tgtgCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGt                                                                       >  1:264943/1‑88 (MQ=255)
                                                                                                                                             gTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCAt  >  2:2066331/1‑150 (MQ=255)
                                                                                                                                                cTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCaa                                                                              >  2:1621220/1‑71 (MQ=255)
                                                                                                                                                cTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCaa                                                                              <  1:1621220/71‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
TCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCAT  >  NC_002947/2063122‑2063412

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: