Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,265,717 G→T 6.9% R245L (CGC→CTC)  folC → bifunctional folypolyglutamate synthase/dihydrofolate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,265,7170GT 6.9% 230.0 / 3.6 72R245L (CGC→CTC) folCbifunctional folypolyglutamate synthase/dihydrofolate synthase
Reads supporting (aligned to +/- strand):  ref base G (31/36);  new base T (3/2);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

TCTGGATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGACCTCCGTAATGT  >  NC_002947/2265571‑2265864
                                                                                                                                                  |                                                                                                                                                   
tCTGGATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGa                                                                                                                                                  <  2:1718409/150‑1 (MQ=255)
tCTGGATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGa                                                                                                                                                  >  2:1116482/1‑150 (MQ=255)
   ggATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCt                                                                                                                                               >  1:1054049/1‑150 (MQ=255)
         gCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCttt                                                                                                                                         >  1:89953/1‑150 (MQ=255)
           ccGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGc                                                                                                                                       <  2:1926946/150‑1 (MQ=255)
            cGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCt                                                                                                                                      <  2:89953/150‑1 (MQ=255)
              ggCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTgg                                                                                                                                    >  2:1160634/1‑150 (MQ=255)
               gCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGa                                                                                                                                   <  1:240309/150‑1 (MQ=255)
                ccGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGAc                                                                                                                                  <  2:1164385/150‑1 (MQ=255)
                ccGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGAc                                                                                                                                  <  1:1489810/150‑1 (MQ=255)
                      cTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGctc                                                                                                                            >  2:829525/1‑148 (MQ=255)
                        ggACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCAt                                                                                                                          <  1:196232/150‑1 (MQ=255)
                           cAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGCTTCTGGCGTGGTACGGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGa                                                                                                                       <  1:235505/150‑1 (MQ=255)
                              ggCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGTGTGGTACTGCCGCTGATGGGGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGaaaa                                                                                                                    <  1:615301/150‑1 (MQ=255)
                                  gccgAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGcc                                                                                                                >  1:213801/1‑150 (MQ=255)
                                     gAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGCTGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTGGCTGGACCTGCCCATGGGAAACGCCACg                                                                                                             >  2:815837/1‑150 (MQ=255)
                                      aaCTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGc                                                                                                            >  1:123848/1‑150 (MQ=255)
                                       aCTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCt                                                                                                           >  2:1083154/1‑150 (MQ=255)
                                                gcCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGc                                                                                                                              >  1:837921/1‑122 (MQ=255)
                                                gcCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGc                                                                                                                              <  2:837921/122‑1 (MQ=255)
                                                 ccACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAg                                                                                                 <  1:116702/150‑1 (MQ=255)
                                                          cTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGCGCAGCTCGACCTGc                                                                                                                                             <  1:1758326/97‑1 (MQ=255)
                                                          cTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGCGCAGCTCGACCTGc                                                                                                                                             >  2:1758326/1‑97 (MQ=255)
                                                          cTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCctg                                                                                        <  1:1083154/150‑1 (MQ=255)
                                                            gcgcGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCctgct                                                                                      <  2:438532/150‑1 (MQ=255)
                                                               cgGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGAt                                                                                   <  1:1140494/150‑1 (MQ=255)
                                                                  tCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGGCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATggg                                                                                <  1:815837/150‑1 (MQ=255)
                                                                    gggAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATgggg                                                                               >  1:1637550/1‑149 (MQ=255)
                                                                    gggAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATgggg                                                                               <  2:1637550/149‑1 (MQ=255)
                                                                            gACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTgg                                                                      <  1:1116482/150‑1 (MQ=255)
                                                                            gACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTgg                                                                      <  1:829525/150‑1 (MQ=255)
                                                                               cTGGCCAGTGCCGAAGGCTGCTGGCGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCATTTGCTGGACCTGCACATGGAAGTGCCCAGGGTGTCGTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAAc                                                                   <  1:1305610/150‑1 (MQ=255)
                                                                                tGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAAcg                                                                  <  2:1760070/150‑1 (MQ=255)
                                                                                 ggCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGCGCTGCGCTACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAAcgc                                                                 <  2:1515853/150‑1 (MQ=255)
                                                                                  gCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAAcgcg                                                                <  2:1178854/150‑1 (MQ=255)
                                                                                   ccAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACgcgg                                                               <  2:1054049/150‑1 (MQ=255)
                                                                                    cAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGaaaa                                                                                                                    <  1:1267749/96‑1 (MQ=255)
                                                                                    cAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGaaaa                                                                                                                    >  2:1267749/1‑96 (MQ=255)
                                                                                    cAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGATAATAGTGCGCAGGTCGAGGTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTACTGATGGGGCTCCCGTGGaagacagc                                                              >  2:781749/1‑144 (MQ=255)
                                                                                         ccGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAg                                                                                                 >  1:177261/1‑110 (MQ=255)
                                                                                         ccGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAg                                                                                                 <  2:177261/110‑1 (MQ=255)
                                                                                         ccGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAg                                                         >  1:693374/1‑150 (MQ=255)
                                                                                          cGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATggg                                                                                >  1:125956/1‑126 (MQ=255)
                                                                                          cGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATggg                                                                                <  2:125956/126‑1 (MQ=255)
                                                                                          cGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCaa                                                          <  1:1297997/148‑1 (MQ=255)
                                                                                          cGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCaa                                                          >  2:1297997/1‑148 (MQ=255)
                                                                                              ggCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGc                                                    <  2:1846390/150‑1 (MQ=255)
                                                                                                 tgctgGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGAATACCGTGGAACGCGGCGCAAGTCCGCCAg                                                 >  1:673/1‑150 (MQ=255)
                                                                                                   ctgGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCtt                                                                                             <  2:838726/104‑1 (MQ=255)
                                                                                                   ctgGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCtt                                                                                             >  1:838726/1‑104 (MQ=255)
                                                                                                          tGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCCGGCgctcttt                                        <  2:1449571/150‑5 (MQ=255)
                                                                                                          tGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGctgctg                                        >  2:1007135/1‑150 (MQ=255)
                                                                                                            gCGTGGTACTGCCGCTGATGGTGCGCAGGTCGCGCAGCTCGACCTg                                                                                                                                              >  2:911354/1‑46 (MQ=21)
                                                                                                            gCGTGGTACTGCCGCTGATGGTGCGCAGGTCGCGCAGCTCGACCTg                                                                                                                                              <  1:911354/46‑1 (MQ=21)
                                                                                                            gCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCttt                                                                                            <  2:962846/96‑1 (MQ=255)
                                                                                                            gCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCttt                                                                                            >  1:962846/1‑96 (MQ=255)
                                                                                                               tGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGACGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATAc                                   >  2:324941/1‑150 (MQ=255)
                                                                                                                    cTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTCCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTa                              >  1:339629/1‑150 (MQ=255)
                                                                                                                     tGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTAt                             >  2:937395/1‑150 (MQ=255)
                                                                                                                         gCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATgggg                                                                               >  2:812035/1‑96 (MQ=255)
                                                                                                                         gCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATgggg                                                                               <  1:812035/96‑1 (MQ=255)
                                                                                                                         gCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACt                         >  1:474187/1‑150 (MQ=255)
                                                                                                                           tGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTgg                       <  2:1337945/150‑1 (MQ=255)
                                                                                                                                    gcAGGTCGCGCAGCTCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGa              >  1:463629/1‑150 (MQ=255)
                                                                                                                                     cAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGAc             >  2:94524/1‑150 (MQ=255)
                                                                                                                                       ggTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGACCt           <  1:1272630/150‑1 (MQ=255)
                                                                                                                                       ggTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGACCt           >  2:1093248/1‑150 (MQ=255)
                                                                                                                                        gTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATggg                                                                                >  2:173060/1‑80 (MQ=255)
                                                                                                                                        gTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATggg                                                                                <  1:173060/80‑1 (MQ=255)
                                                                                                                                        gTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGACCTc          >  1:1837250/1‑150 (MQ=255)
                                                                                                                                              cTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTTTATACCCGTATCGCTGGGCGTCGCGACCGCCGaaaa    >  2:1859478/1‑149 (MQ=255)
                                                                                                                                               tGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGACCTCCGTAAtg   <  1:937395/150‑1 (MQ=255)
                                                                                                                                                gcgcGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGACCTCCGTAAtgt  <  2:1837250/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCTGGATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTTTCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGACCTCCGTAATGT  >  NC_002947/2265571‑2265864

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: