Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,258,563:1 +G 8.3% coding (184/1005 nt) asd → aspartate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,258,5631.G 8.3% 225.6 / 15.2 60coding (184/1005 nt)asdaspartate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (32/23);  new base G (3/2);  total (35/25)
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.27e-01

CCACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAAC  >  NC_002947/2258414‑2258712
                                                                                                                                                      |                                                                                                                                                      
ccACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCt                                                                                                                                                         >  2:298545/1‑150 (MQ=255)
ccACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCt                                                                                                                                                         >  1:1363556/1‑150 (MQ=255)
 cACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCtt                                                                                                                                                       <  1:245753/150‑1 (MQ=255)
  aCCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTc                                                                                                                                                      <  1:1884596/150‑1 (MQ=255)
      gCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAATTCGACAGCTTCGACt                                                                                                                                                 >  1:856456/1‑150 (MQ=255)
          cGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGAAGCTGTCGACTTCGc                                                                                                                                             <  2:29932/150‑1 (MQ=255)
              ggTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAg                                                                                                                                         <  2:1640408/150‑1 (MQ=255)
                     ccTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGc                                                                                                                                  >  2:232470/1‑150 (MQ=255)
                          gTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGcc                                                                                                                             <  1:1712436/150‑1 (MQ=255)
                               gATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCcttctt                                                                                                                        <  2:856456/150‑1 (MQ=255)
                                aTCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCttcttc                                                                                                                       >  2:172100/1‑150 (MQ=255)
                                     cgaggcaCTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTACGCGAAGTCGACAGCTTCGCCTTCGCCCAGGTCCAGCTGGCCTTCTTCGCCAc                                                                                                                  >  2:1765673/7‑150 (MQ=255)
                                      gaagaaCTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTg                                                                                                                                <  2:1392042/135‑1 (MQ=255)
                                      gaagaaCTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTg                                                                                                                                >  1:1392042/1‑135 (MQ=255)
                                       aagaaCTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCg                                                                                                                >  1:1863188/1‑150 (MQ=255)
                                            cTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGcgcc                                                                                                           >  2:523974/1‑150 (MQ=255)
                                                 aTTGCCGGTTGCCACTTTGCACCTGCTGGCCAGCATGGAATGGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAGACTGAGGGGGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCGTCTTCGCCACCGGCGCCGCCGt                                                                                                      <  2:1745972/150‑1 (MQ=255)
                                                  ttCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGAc                                                                                                                                                  >  2:399396/1‑105 (MQ=255)
                                                  ttCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGAc                                                                                                                                                  <  1:399396/105‑1 (MQ=255)
                                                     ccGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGc                                                                                                  <  2:1015701/150‑1 (MQ=255)
                                                             cACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGAACTGTTAGGACGGATAAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCtt                                                                                          >  2:1911045/1‑150 (MQ=255)
                                                              aCTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCttt                                                                                         <  1:298545/150‑1 (MQ=255)
                                                                  tGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCg                                                                                     <  1:1294044/150‑1 (MQ=255)
                                                                       ctgctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCg                                                                                     >  2:2591/1‑145 (MQ=255)
                                                                       ctgctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCg                                                                                     <  1:2591/145‑1 (MQ=255)
                                                                            ggCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGGCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCt                                                                           >  1:385809/1‑150 (MQ=255)
                                                                              ccAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGc                                                                         >  1:858914/1‑150 (MQ=255)
                                                                                 gCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGTCAAGGCGCTGCAgg                                                                      >  1:1626398/1‑150 (MQ=255)
                                                                                   ctggatccggcggGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGCCTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCg                                                                    <  1:1461657/143‑1 (MQ=255)
                                                                                    tGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCgg                                                                   <  1:1126557/150‑1 (MQ=255)
                                                                                    tGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCACCGCCGTCAGACGCTGCTGTGTCGGCGAGGCGCTGCAGGCgg                                                                   >  2:575742/1‑150 (MQ=255)
                                                                                         tCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCGAGGGCAAGCTTGCCGTCTTCGCCACCGGCGCCTCCGGCATCCGCAGCTTTGCCGTCAATGCGCTGCAGGCGGGCTgg                                                              >  2:1841630/1‑149 (MQ=255)
                                                                                            gcggGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACg                                                           <  2:1363556/150‑1 (MQ=255)
                                                                                              gggCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGt                                                         <  1:523974/150‑1 (MQ=255)
                                                                                                  agcagcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATc                                                     <  2:1863188/150‑1 (MQ=255)
                                                                                                     agcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGAc                                                  >  2:1130037/1‑150 (MQ=255)
                                                                                                      gcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAACTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGAcc                                                 >  1:1724947/1‑150 (MQ=255)
                                                                                                          gATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTc                                             >  2:755655/1‑150 (MQ=255)
                                                                                                              ttCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCAACGTCACCTCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCcggc                                         >  1:1554002/1‑150 (MQ=255)
                                                                                                                   tGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCt                                    <  2:858914/150‑1 (MQ=255)
                                                                                                                   tGGCAAGAAACAGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGATCATCGACCTGTCCGGCGGccc                                    >  2:715876/1‑149 (MQ=255)
                                                                                                                     gCAAGAAACTGAAGGTGCGCGAAGTCGACAGATTCAGGATTATGCTAGGCAACATAAGTTCCTTCTCCACCGACACCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTgg                                  >  1:514203/1‑150 (MQ=255)
                                                                                                                       aaGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGgac                                >  1:1520562/1‑150 (MQ=255)
                                                                                                                       aaGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGgac                                >  1:1645874/1‑150 (MQ=255)
                                                                                                                         gAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGgacga                              >  2:359512/1‑150 (MQ=255)
                                                                                                                             cTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgcc                          >  2:185586/1‑150 (MQ=255)
                                                                                                                              tGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccc                         >  2:1804638/1‑150 (MQ=255)
                                                                                                                                    tgggcgAAGTCGAAGCTGTCGACTTCGCGCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGGCATCGACCGTTCCGGCGGCCTGGACGACgccctggcac                   >  2:1260300/4‑148 (MQ=255)
                                                                                                                                    tGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggccc                   >  2:806264/1‑150 (MQ=255)
                                                                                                                                      cgcgAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggccctg                 <  1:172100/150‑1 (MQ=255)
                                                                                                                                       gcgAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGccctggccctgg                >  1:159997/1‑150 (MQ=255)
                                                                                                                                         gAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTg              >  2:614233/1‑150 (MQ=255)
                                                                                                                                         gAAGTCGAAGCTGTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccct                        <  2:953364/140‑1 (MQ=255)
                                                                                                                                         gAAGTCGAAGCTGTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccct                        >  1:953364/1‑140 (MQ=255)
                                                                                                                                         gAAGTCGAAGCTGTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTg              >  1:1759830/1‑150 (MQ=255)
                                                                                                                                            gTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGATGTACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCg           <  1:1573590/150‑1 (MQ=255)
                                                                                                                                               gACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa        <  1:614233/150‑1 (MQ=255)
                                                                                                                                               gACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa        <  1:1911045/150‑1 (MQ=255)
                                                                                                                                               gACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa        >  1:1687627/1‑150 (MQ=255)
                                                                                                                                                aCAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAg       >  1:5742/1‑150 (MQ=255)
                                                                                                                                                    cTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCaa                                                                                 <  2:1344176/72‑1 (MQ=255)
                                                                                                                                                    cTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCaa                                                                                 >  1:1344176/1‑72 (MQ=255)
                                                                                                                                                    cTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCATACGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCaa   <  1:1841630/150‑1 (MQ=255)
                                                                                                                                                     ttCGACTTCGCCCAGGTCGTGGTTGCCCACAGCGACACGCTGGCCGTCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAAc  <  2:1645874/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CCACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCTTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAAC  >  NC_002947/2258414‑2258712

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: