Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,081,371 C→T 6.1% T72T (ACG→ACA)  lysC ← aspartokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,081,3710CT 6.1% 192.1 / 2.3 65T72T (ACG→ACA) lysCaspartokinase
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base T (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CTGGTCGTCGATTTGCAGAATGCGCGCCTTGTTGTGCGAGCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAACCACCAGGTCGTCGCCCGCTTCACGGTGTTTCTTGACCTTTTCGG  >  NC_002947/5081230‑5081519
                                                                                                                                             |                                                                                                                                                    
ctGGTCGTCGATTTGCAGAATGCGCGCCTTGTTGTGCGAGCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGTTGACCGGCTCACCCGTCGacac                                                                                                                                              >  2:1859639/1‑150 (MQ=255)
     cgtcgATTTGCAGAATGCGCGCCTTGTTGTGCGAGCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCa                                                                                                                                         >  2:1744984/1‑150 (MQ=255)
      gtcgATTTGCAGAATGCGCGCCTTGTTGTGCGAGCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCAc                                                                                                                                        >  2:1363445/1‑150 (MQ=255)
       tcgATTTGCAGAATGCGCGCCTTGTTGTGCGAGCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCTGTCGACACGATCACg                                                                                                                                       >  2:1736636/1‑150 (MQ=255)
       tcgATTTGCAGAATGCGCGCCTTGTTGTGCGAGCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCTGTCGACACGATCACg                                                                                                                                       >  2:1735716/1‑150 (MQ=255)
                       gcgcCTTGTTGTGCGAGCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGaa                                                                                                                       >  2:836487/1‑150 (MQ=255)
                       gcgcCTTGTTGTGCGAGCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGaa                                                                                                                       >  1:584249/1‑150 (MQ=255)
                       gcgcCTTGTTGTGCGAGCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGaa                                                                                                                       >  2:1423075/1‑150 (MQ=255)
                                  tgCGAGCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATc                                                                                                                                          >  2:352402/1‑120 (MQ=255)
                                  tgCGAGCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATc                                                                                                                                          <  1:352402/120‑1 (MQ=255)
                                     gAGCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTAGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTg                                                                                                         <  1:977803/150‑1 (MQ=255)
                                       gCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGAt                                                                                                       <  1:1739282/150‑1 (MQ=255)
                                             ggTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAATCGGCTGATCGGTGATCTGCtt                                                                                                 >  1:324745/1‑150 (MQ=255)
                                              gTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTg                                                                                                <  1:212233/150‑1 (MQ=255)
                                                     tACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGTTCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGt                                                                                         <  1:541042/150‑1 (MQ=255)
                                                     tACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCa                                                                                                                              <  2:1435385/113‑1 (MQ=255)
                                                     tACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCa                                                                                                                              >  1:1435385/1‑113 (MQ=255)
                                                     tACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGt                                                                                         >  1:115775/1‑150 (MQ=255)
                                                      aCGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTc                                                                                        <  1:1744984/150‑1 (MQ=255)
                                                        gCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCAc                                                                                                                                        <  2:1533494/100‑1 (MQ=255)
                                                        gCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCAc                                                                                                                                        >  1:1533494/1‑100 (MQ=255)
                                                        gCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGa                                                                                      >  1:1237617/1‑150 (MQ=255)
                                                           ccTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCa                                                                                   <  1:1423075/150‑1 (MQ=255)
                                                            cTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAg                                                                                  <  1:836487/150‑1 (MQ=255)
                                                             tGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAgg                                                                                 <  2:115775/150‑1 (MQ=255)
                                                               gTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGAt                                                                                                       <  2:92260/126‑1 (MQ=255)
                                                               gTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGAt                                                                                                       >  1:92260/1‑126 (MQ=255)
                                                                  gCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGAtt                                                                            >  1:1803103/1‑150 (MQ=255)
                                                                  gCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGAtt                                                                            <  1:1363445/150‑1 (MQ=255)
                                                                     ggTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGTCACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGACCAGGTCGATCAGGCGATTGGt                                                                         >  1:678249/1‑150 (MQ=255)
                                                                     ggTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCgg                                                                                                                   <  1:1324379/108‑1 (MQ=255)
                                                                     ggTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCgg                                                                                                                   >  2:1324379/1‑108 (MQ=255)
                                                                     ggTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGt                                                                         >  1:1803315/1‑150 (MQ=255)
                                                                     ggTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGt                                                                         >  1:1695900/1‑150 (MQ=255)
                                                                      gtgtAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCtt                                                                                                 >  2:714778/1‑125 (MQ=255)
                                                                      gtgtAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCtt                                                                                                 <  1:714778/125‑1 (MQ=255)
                                                                           ggACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCAcc                                                                   <  2:584249/150‑1 (MQ=255)
                                                                            gACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCg                                                                  >  2:1238577/1‑150 (MQ=255)
                                                                            gACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTCGCTGCACCg                                                                  >  1:81996/1‑150 (MQ=255)
                                                                               acCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAg                                                                                           <  2:1460881/122‑1 (MQ=255)
                                                                               acCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAg                                                                                           >  1:1460881/1‑122 (MQ=255)
                                                                               acCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGt                                                                                         >  1:256368/1‑124 (MQ=255)
                                                                               acCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGt                                                                                         <  2:256368/124‑1 (MQ=255)
                                                                                      gCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTc                                                                                        >  2:1704365/1‑118 (MQ=255)
                                                                                      gCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTc                                                                                        <  1:1704365/118‑1 (MQ=255)
                                                                                        acacCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGAc                                                      >  1:1613189/1‑150 (MQ=255)
                                                                                          accacGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATc                                                                                                      <  2:252538/100‑1 (MQ=255)
                                                                                          accacGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATc                                                                                                      >  1:252538/1‑100 (MQ=255)
                                                                                          accacGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATc                                                                                    >  1:1595144/1‑118 (MQ=255)
                                                                                          accacGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATc                                                                                    <  2:1595144/118‑1 (MQ=255)
                                                                                          accacGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAg                                                    >  1:405599/1‑150 (MQ=255)
                                                                                            cacGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGca                                                  >  1:200455/1‑150 (MQ=255)
                                                                                                cTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGCCAGcacatc                                              <  2:1449868/150‑3 (MQ=255)
                                                                                                    aTCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAaccacc                                          >  2:1161593/1‑150 (MQ=255)
                                                                                                     tCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGc                                                                                              <  2:1581974/97‑1 (MQ=255)
                                                                                                     tCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGc                                                                                              >  1:1581974/1‑97 (MQ=255)
                                                                                                       aaGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAACCACCAgg                                       <  2:1237617/150‑1 (MQ=255)
                                                                                                       aaGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAACCACCAgg                                       <  2:1803103/150‑1 (MQ=255)
                                                                                                        aGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAACCACCAGgt                                      <  2:1803315/150‑1 (MQ=255)
                                                                                                           ccATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAACCACCAGgtcgt                                   >  1:1653073/1‑150 (MQ=255)
                                                                                                           ccATGGTCAGCAGGGCAATGGTGACCTCCTCACCCCTCCACACGATCACCTACAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGCCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAACCACCAGgtcgt                                   <  1:1859639/150‑1 (MQ=255)
                                                                                                            cATGGTCAGCAGGGCAATGGTGACCTGCTCACCTGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAACCACCAGgtcgtc                                  <  1:1735716/150‑1 (MQ=255)
                                                                                                            cATGGTCAGCAGGGCAATGGTGACCTGCTCACCTGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAACCACCAGgtcgtc                                  <  1:1736636/150‑1 (MQ=255)
                                                                                                                                  aCCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAACCACCAGGTCGTCGCCCGCTTCACGGTGTTTCTTg            >  2:786652/1‑150 (MQ=255)
                                                                                                                                   ccggcgcACCCGTCGACACGAGCACGTCCAGTTCACGAGGAACCGGCTGATTGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAACCACCAGGTCGTCGCCCGCTTCACGGTGTTTCTTGa           <  2:324745/144‑1 (MQ=255)
                                                                                                                                           cccGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAACCACCAGGTCGTCGCCCGCTTCACGGTGTTTATTGACCTTTTCg   <  1:1161593/150‑1 (MQ=255)
                                                                                                                                            ccGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGcc                                                                                             >  2:1111199/1‑59 (MQ=255)
                                                                                                                                            ccGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGcc                                                                                             <  1:1111199/59‑1 (MQ=255)
                                                                                                                                            ccGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAACCACCAGGTCGTCGCCCGCTTCACGGTGTTTCTTGACCTTTTCgg  >  2:560620/1‑150 (MQ=255)
                                                                                                                                             cGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCAcc                                                                   <  1:1470105/84‑1 (MQ=255)
                                                                                                                                             cGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCAcc                                                                   >  2:1470105/1‑84 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CTGGTCGTCGATTTGCAGAATGCGCGCCTTGTTGTGCGAGCTGTCGGTGAGAATACGCACCTGGTTGCCGGTGTAGGACACCGCTGGCACACCACGCTTGATCAAGGCCATGGTCAGCAGGGCAATGGTGACCTGCTCACCCGTCGACACGATCACGTCCAGTTCACGAGGAACCGGCTGATCGGTGATCTGCTTGGCCAGGTCGATCAGGCGATTGGTTTCACCGCTCATGGCCGACAGCACAACCACCAGGTCGTCGCCCGCTTCACGGTGTTTCTTGACCTTTTCGG  >  NC_002947/5081230‑5081519

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