Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,244,402 G→C 7.1% S221R (AGC→AGG aruH ← arginine‑‑pyruvate transaminase AruH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,244,4020GC 7.1% 353.4 / 14.5 114S221R (AGC→AGGaruHarginine‑‑pyruvate transaminase AruH
Reads supporting (aligned to +/- strand):  ref base G (54/51);  new base C (3/5);  total (58/56)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

Only 100 of 115 total aligned reads displayed.

GCCGTAGAGCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGT  >  NC_002947/4244258‑4244550
                                                                                                                                                |                                                                                                                                                      
gCCGTAGAGCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGcc                                                                                                                                                   <  2:1570935/150‑1 (MQ=255)
    tAGAGCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCCCCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGACCATGCCAGGCGGGCTGCgggag                                                                                                                                               <  2:1133633/150‑3 (MQ=255)
    tAGAGCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGCGCGCTCGGCCATGCCCGGCAGGCTGCCgggg                                                                                                                                               <  2:1480911/150‑1 (MQ=255)
      gagCATGCGAAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGTCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGCCTGCCgggcat                                                                                                                                             <  1:1613557/150‑4 (MQ=255)
      gagCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCt                                                                                                                                             <  2:1014394/150‑1 (MQ=255)
      gagCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGCCTGCCgggcat                                                                                                                                             <  2:1841394/150‑4 (MQ=255)
      gagCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGCCTGCCgggcatg                                                                                                                                            <  2:2050764/150‑5 (MQ=255)
              acaGGGCGGGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCt                                                                                                                                    <  2:6321/150‑1 (MQ=255)
              acaGGGCGAGGTTTTCCAGATGCGTGGGCAGTGCTGTCGAGCCGGCGACCCACCCCATGCACCAGCGGGTCATGGCATGCGATTTGGACAGGCTGTTGAAGGTTGCAGTGCGCTCGGCCATGGCCGGCAGGGTGCGGGGGCTGACATggt                                                                                                                                    >  2:1175455/1‑148 (MQ=255)
                  ggCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGcc                                                                                                                                >  1:1338314/1‑150 (MQ=255)
                   gCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCt                                                                                                                               <  2:422677/150‑1 (MQ=255)
                             ccAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATAcagca                                                                                                                    >  2:1685009/1‑150 (MQ=255)
                                  tGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCg                                                                                                               <  1:1682311/150‑1 (MQ=255)
                                   gCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCgggg                                                                                                                                               <  2:2285275/119‑1 (MQ=255)
                                   gCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCgggg                                                                                                                                               >  1:2285275/1‑119 (MQ=255)
                                    cGTGGCCAGTGCTGTCGAGCCGACCACACACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCt                                                                                                             <  2:2071833/150‑1 (MQ=255)
                                       ggCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTa                                                                                                          <  2:1851347/150‑1 (MQ=255)
                                       ggCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTa                                                                                                          <  2:1054561/150‑1 (MQ=255)
                                        gCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAg                                                                                                         >  1:912536/1‑150 (MQ=255)
                                         ccAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCGAGCTGCTGTatg                                                                                                        <  2:605768/150‑3 (MQ=255)
                                             tGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTc                                                                                                    >  1:60791/1‑150 (MQ=255)
                                              gCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGt                                                                                                           >  1:2220261/1‑142 (MQ=255)
                                              gCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGt                                                                                                           <  2:2220261/142‑1 (MQ=255)
                                               cTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGt                                                                                                  >  2:658533/1‑150 (MQ=255)
                                               cTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGt                                                                                                  <  1:1685009/150‑1 (MQ=255)
                                                tGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTc                                                                                                 <  1:1107310/150‑1 (MQ=255)
                                                  tCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCag                                                                                               >  2:586528/1‑150 (MQ=255)
                                                     aGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAgaga                                                                                            >  1:435849/1‑150 (MQ=255)
                                                        cGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCa                                                                                         >  1:831273/1‑150 (MQ=255)
                                                         gACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACAt                                                                                                                                       >  2:485828/1‑104 (MQ=255)
                                                         gACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACAt                                                                                                                                       <  1:485828/104‑1 (MQ=255)
                                                         gACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAg                                                                                        >  1:66801/1‑150 (MQ=255)
                                                          accaccCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGc                                                                                       <  2:639588/150‑1 (MQ=255)
                                                          accaccCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACAGTCGCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGc                                                                                       <  1:658533/150‑1 (MQ=255)
                                                            cacccaccCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCca                                                                                     <  2:912536/150‑1 (MQ=255)
                                                            cacccaccCCATGCGCCAGCCGGTCATGGCATGCGATATGTACAGGCTGTTCAACGTTGCAATGCGCTCGGCAACAAGCCCACGCTGGTACGTATCGGTATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCATAAGAGATGAGCca                                                                                     >  1:1151676/1‑150 (MQ=255)
                                                             acccacccCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCa                                                                                                >  2:2124004/1‑138 (MQ=255)
                                                             acccacccCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCa                                                                                                <  1:2124004/138‑1 (MQ=255)
                                                              cccacccCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTc                                                                                                                                   >  2:460753/1‑103 (MQ=255)
                                                              cccacccCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTc                                                                                                                                   <  1:460753/103‑1 (MQ=255)
                                                               ccacccCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAg                                                                                  >  1:1571918/1‑150 (MQ=255)
                                                                 acccCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGAt                                                                                >  1:2183386/1‑150 (MQ=255)
                                                                   cccATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTa                                                                                                          >  2:2047917/1‑122 (MQ=255)
                                                                   cccATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTa                                                                                                          <  1:2047917/122‑1 (MQ=255)
                                                                   cccATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGAGATCGGCCATCAACGCCACGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCa                                                                              >  2:2078899/1‑150 (MQ=255)
                                                                    ccATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTa                                                                                                          <  1:2408212/121‑1 (MQ=255)
                                                                    ccATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTa                                                                                                          >  2:2408212/1‑121 (MQ=255)
                                                                     cATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCa                                                                                         >  2:8917/1‑137 (MQ=255)
                                                                      aTGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTc                                                                                                    >  2:900511/1‑125 (MQ=255)
                                                                      aTGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCAGGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTc                                                                                                    <  1:900511/125‑1 (MQ=255)
                                                                             aGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGGGCGCTCGGCCATGACCGGCAGGCAGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCAAAGATCATCGGCAACGc                                                                    >  1:965228/1‑150 (MQ=255)
                                                                              gCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGca                                                                   <  1:509208/150‑1 (MQ=255)
                                                                                cGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCaga                                                                                              <  2:1298633/121‑1 (MQ=255)
                                                                                cGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCaga                                                                                              >  1:1298633/1‑121 (MQ=255)
                                                                                 ggTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAg                                                                >  1:1091214/1‑150 (MQ=255)
                                                                                      tGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACCGCAGCTCGCTGTAGACCTCGTCAGACATTAGCCACAGATCATGGGCAACGCACAGTTCgg                                                           >  2:1141473/1‑150 (MQ=255)
                                                                                       ggCATGCGAGTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCTACGCACAGTTgggg                                                          >  2:1651594/1‑149 (MQ=255)
                                                                                        gCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGcc                                                         >  2:1827436/1‑150 (MQ=255)
                                                                                          aTGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAg                                                       >  2:333296/1‑150 (MQ=255)
                                                                                           tGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAgc                                                      >  1:2027320/1‑150 (MQ=255)
                                                                                               aTTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATc                                                                                          <  2:1928303/110‑1 (MQ=255)
                                                                                               aTTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATc                                                                                          >  1:1928303/1‑110 (MQ=255)
                                                                                                       aGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTg                                          <  2:1861309/150‑1 (MQ=255)
                                                                                                        ggCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCa                                                                                         >  1:66754/1‑102 (MQ=255)
                                                                                                         gCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGgc                                        >  2:1157099/1‑150 (MQ=255)
                                                                                                           tGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGgcgc                                      >  2:197863/1‑150 (MQ=255)
                                                                                                           tGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGgcgc                                      >  2:2022162/1‑150 (MQ=255)
                                                                                                           tGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGCGGCTCACTTTCTCGCTTTCATACAGCAGCTCGCTGTAGACCTCGTCAGCGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCACCTGGGAGgcgc                                      <  2:677958/150‑1 (MQ=255)
                                                                                                             ttCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTTgcgcgc                                    <  1:1240344/150‑1 (MQ=255)
                                                                                                             ttCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGgcgcgc                                    <  2:435849/150‑1 (MQ=255)
                                                                                                             ttCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTCCCGCGGCTGACCTGATGGCCTTCGTACAGCAGCTCGCTGTCGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGgcgcgc                                    >  2:214407/1‑150 (MQ=255)
                                                                                                                aaCGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCcaca                                                                                   >  1:1915205/1‑100 (MQ=255)
                                                                                                                aaCGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGgcgc                                      >  2:570493/1‑145 (MQ=255)
                                                                                                                aaCGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGgcgc                                      <  1:570493/145‑1 (MQ=255)
                                                                                                                aaCGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGc                                 >  1:1936311/1‑150 (MQ=255)
                                                                                                                 aCGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCa                                <  1:197863/150‑1 (MQ=255)
                                                                                                                      gCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCaga                                                                                              <  2:1870581/83‑1 (MQ=255)
                                                                                                                      gCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCaga                                                                                              >  1:1870581/1‑83 (MQ=255)
                                                                                                                           gcgcTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGctcc                      <  2:831273/150‑1 (MQ=255)
                                                                                                                             gcTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTg                    >  1:964422/1‑150 (MQ=255)
                                                                                                                              cTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGa                   <  1:1819042/150‑1 (MQ=255)
                                                                                                                                cGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAg                 <  2:1589615/150‑1 (MQ=255)
                                                                                                                                cGGCCAGGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCTCTGAAg                 <  1:466264/150‑1 (MQ=255)
                                                                                                                                 ggCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAgg                <  1:1651594/150‑1 (MQ=255)
                                                                                                                                   ccATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGt              >  2:242046/1‑150 (MQ=255)
                                                                                                                                    cATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGtt             >  2:284875/1‑150 (MQ=255)
                                                                                                                                    cATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGtt             <  2:1338314/150‑1 (MQ=255)
                                                                                                                                      gccccGGCAGCCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAg                 <  2:1861308/142‑1 (MQ=255)
                                                                                                                                      gccccGGCAGCCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAg                 >  1:1861308/3‑144 (MQ=255)
                                                                                                                                       ccccGGCAGCCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGtt             <  1:768279/146‑1 (MQ=255)
                                                                                                                                       ccccGGCAGCCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGtt             >  2:768279/2‑147 (MQ=255)
                                                                                                                                       ccccGGCAGCCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTtgtg          >  1:799577/2‑150 (MQ=255)
                                                                                                                                         ccGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACagcagc                                                                                                                  >  1:1030195/1‑44 (MQ=255)
                                                                                                                                         ccGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACagcagc                                                                                                                  <  2:1030195/44‑1 (MQ=255)
                                                                                                                                         ccGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATgg                                                                           <  1:1935155/83‑1 (MQ=255)
                                                                                                                                         ccGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATgg                                                                           >  2:1935155/1‑83 (MQ=255)
                                                                                                                                         ccGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCGCTGAAGGGTTGTGCg        >  2:1683066/1‑150 (MQ=255)
                                                                                                                                          cGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCgg       <  2:1333833/150‑1 (MQ=255)
                                                                                                                                           ggCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCggg      <  2:2027320/150‑1 (MQ=255)
                                                                                                                                               ggCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGt  >  2:1865752/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
GCCGTAGAGCATGCACAGGGCGAGGTTTTCCAGATGCGTGGCCAGTGCTGTCGAGCCGACCACCCACCCCATGCGCCAGCCGGTCATGGCATGCGATTTGGACAGGCTGTTCAACGTTGCAGTGCGCTCGGCCATGCCCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGT  >  NC_002947/4244258‑4244550

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