Predicted mutation
evidence seq id position mutation annotation gene description
RA AC_000091 2,684,480 (A)7→8 intergenic (‑317/‑11) glyA ← / → hmp serine hydroxymethyltransferase/fused nitric oxide dioxygenase and dihydropteridine reductase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AC_0000912,684,4731.A99.1% 472.4 / ‑3.0 115intergenic (‑310/‑18)glyA/hmpserine hydroxymethyltransferase/fused nitric oxide dioxygenase and dihydropteridine reductase 2
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (66/48);  total (66/49)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-01

Only 100 of 134 total aligned reads displayed.

ATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCA  >  AC_000091/2684384‑2684556
                                                                                          |                                                                                   
atcaCATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGc                                                                                      <  1:892835/90‑1 (MQ=255)
atcaCATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGc                                                                                      <  2:1597138/90‑1 (MQ=255)
  cacaTTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCaa                                                                                   <  2:3434800/90‑1 (MQ=255)
   acaTTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCaaa                                                                                  <  2:3062471/90‑1 (MQ=255)
    caTTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCaaaa                                                                                 >  2:13921/1‑90 (MQ=255)
     aTTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCaaaaa                                                                                >  1:2605038/1‑90 (MQ=255)
     aTTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCaaaaa                                                                                >  2:1846501/1‑90 (MQ=255)
     aTTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCaaaaa                                                                                >  2:107497/1‑90 (MQ=255)
       tctttctGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCaaaaaaa                                                                              >  1:2540517/1‑90 (MQ=255)
         tttctGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCaaaaaaaag                                                                            >  1:676659/1‑88 (MQ=255)
         tttctGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCaaaaaaaag                                                                            <  1:1867419/90‑3 (MQ=255)
           tctGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCaaaaaaaagga                                                                          <  2:563818/90‑5 (MQ=255)
             tGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAg                                                                         <  1:2764130/90‑1 (MQ=255)
             tGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAg                                                                         <  1:355627/90‑1 (MQ=255)
             tGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAg                                                                         >  2:1185565/1‑90 (MQ=255)
              gAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGa                                                                        <  2:1529937/90‑1 (MQ=255)
              gAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGa                                                                        >  2:1886941/1‑90 (MQ=255)
                aaaaCACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGGcc                                                                      >  2:1577244/1‑90 (MQ=255)
                aaaaCACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGAcc                                                                      <  1:1222374/90‑1 (MQ=255)
                 aaaCACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCa                                                                     >  2:351636/1‑90 (MQ=255)
                 aaaCACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCa                                                                     >  1:2216914/1‑90 (MQ=255)
                  aaCACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCat                                                                    <  2:303212/90‑1 (MQ=255)
                  aaCACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCat                                                                    >  2:2609248/1‑90 (MQ=255)
                  aaCACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCat                                                                    >  2:31307/1‑90 (MQ=255)
                   acatCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCata                                                                   <  2:3445269/90‑1 (MQ=255)
                    cacCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCatat                                                                  >  2:721636/1‑90 (MQ=255)
                    cacCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCatat                                                                  >  1:1850240/1‑90 (MQ=255)
                     acCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATg                                                                 <  1:2286467/90‑1 (MQ=255)
                     acCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATg                                                                 >  1:1599699/1‑90 (MQ=255)
                      ccAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGc                                                                >  1:128768/1‑90 (MQ=255)
                      ccAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGc                                                                <  1:94831/90‑1 (MQ=255)
                       cAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCt                                                               >  1:1172046/1‑90 (MQ=255)
                       cAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCt                                                               >  1:764760/1‑90 (MQ=255)
                        aaaGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCtt                                                              <  1:840401/90‑1 (MQ=255)
                        aaaGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACACCAATTAAGATGCAAAAAAAAGGAAGACCATATGCtt                                                              >  1:2708812/1‑90 (MQ=255)
                          aGAACCATTTACATTGCATGGCTATTTTTTATAAGATGCATTTGAGATACATCAACTTAGATGCAAAAAAAAGGAAGACCATATGCTTGa                                                            <  2:3230937/90‑1 (MQ=255)
                          aGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGa                                                            >  2:2503790/1‑90 (MQ=255)
                           gAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGAc                                                           <  1:631177/90‑1 (MQ=255)
                            aaCCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACg                                                          <  1:617718/90‑1 (MQ=255)
                            aaCCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACg                                                          >  2:1487964/1‑90 (MQ=255)
                            aaCCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGc                                                         <  2:1377702/90‑1 (MQ=255)
                              ccATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCt                                                        >  1:1162528/1‑90 (MQ=255)
                               cATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTc                                                       >  2:2683255/1‑90 (MQ=255)
                               cATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTc                                                       >  2:2824881/1‑90 (MQ=255)
                               cATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTc                                                       >  1:1147926/1‑90 (MQ=255)
                                aTTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCa                                                      >  2:3174549/1‑90 (MQ=255)
                                aTTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCa                                                      >  2:787977/1‑90 (MQ=255)
                                 tttACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCaa                                                     <  2:246375/90‑1 (MQ=255)
                                     cATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCa                                                 <  2:1637844/90‑1 (MQ=255)
                                      aTTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCAt                                                <  1:2131398/90‑1 (MQ=255)
                                      aTTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCAt                                                <  1:2122472/90‑1 (MQ=255)
                                       ttGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATc                                               <  1:2228654/90‑1 (MQ=255)
                                        tGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCg                                              <  2:2928055/90‑1 (MQ=255)
                                        tGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCg                                              <  1:1664745/90‑1 (MQ=255)
                                         gCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGc                                             >  2:958485/1‑90 (MQ=255)
                                         gCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGc                                             >  1:501651/1‑90 (MQ=255)
                                          cAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCt                                            >  2:2417019/1‑90 (MQ=255)
                                          cAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCt                                            <  1:615137/90‑1 (MQ=255)
                                          cAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCACt                                            <  1:305017/90‑1 (MQ=255)
                                            gggCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTAc                                          >  1:2139143/1‑90 (MQ=255)
                                              gCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAg                                        >  1:3455231/1‑90 (MQ=255)
                                              gCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAg                                        >  1:512483/1‑90 (MQ=255)
                                               cTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGt                                       >  1:1539738/1‑90 (MQ=255)
                                                 aTTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTaa                                     <  2:523938/90‑1 (MQ=255)
                                                  ttttttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTaaa                                    <  1:3117156/90‑1 (MQ=255)
                                                    ttttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAg                                  <  1:173372/90‑1 (MQ=255)
                                                    ttttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAg                                  <  1:2317837/90‑1 (MQ=255)
                                                    ttttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAg                                  >  2:1601956/1‑90 (MQ=255)
                                                     tttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGc                                 <  1:3119496/90‑1 (MQ=255)
                                                     tttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGc                                 >  2:1207720/1‑90 (MQ=255)
                                                      ttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGcc                                <  1:1143089/90‑1 (MQ=255)
                                                      ttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGcc                                <  2:826688/90‑1 (MQ=255)
                                                        ataAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGccac                              <  2:3368278/90‑1 (MQ=255)
                                                         taAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGccacc                             <  2:2157400/90‑1 (MQ=255)
                                                            gATGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATc                          >  2:2438646/1‑90 (MQ=255)
                                                             aTGCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATcc                         <  2:1725797/90‑1 (MQ=255)
                                                               gCATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCt                       >  2:1177440/1‑90 (MQ=255)
                                                                cATTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCtt                      <  1:502779/90‑1 (MQ=255)
                                                                 aTTTGAGATACATCANTTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCttt                     >  1:3208258/1‑90 (MQ=255)
                                                                 aTTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCGTCCCttt                     >  1:2066304/1‑90 (MQ=255)
                                                                 aTTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCttt                     <  2:1355632/90‑1 (MQ=255)
                                                                 aTTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCttt                     >  2:954008/1‑90 (MQ=255)
                                                                 aTTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCttt                     >  2:569133/1‑90 (MQ=255)
                                                                 aTTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCttt                     <  1:490414/90‑1 (MQ=255)
                                                                 aTTTGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCttt                     >  2:257807/1‑90 (MQ=255)
                                                                   ttGAGATACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTAc                   >  1:1319944/1‑90 (MQ=255)
                                                                     gagaTACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACtg                 <  2:2490145/90‑1 (MQ=255)
                                                                      agaTACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACtgg                <  1:1467774/90‑1 (MQ=255)
                                                                         tACATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACtggtgg             >  2:1658882/1‑90 (MQ=255)
                                                                          aCATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGa            >  1:2829142/1‑90 (MQ=255)
                                                                           cATCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGaa           >  1:2256826/1‑90 (MQ=255)
                                                                            aTCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGaaa          <  1:3377308/90‑1 (MQ=255)
                                                                             tCAATTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAAc         <  2:2713142/90‑1 (MQ=255)
                                                                               aaTTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACgg       <  2:3061130/90‑1 (MQ=255)
                                                                                aTTAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACggg      >  1:776220/1‑90 (MQ=255)
                                                                                 ttAAGATGCAAAAAAAAGGAAGACCATATGCTCGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTTCTGGTGGAAACgggg     <  2:14930/90‑1 (MQ=255)
                                                                                  tAAGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGc    >  2:1888896/1‑90 (MQ=255)
                                                                                   aaGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGcc   >  1:1849660/1‑90 (MQ=255)
                                                                                    aGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCa  >  2:458221/1‑90 (MQ=255)
                                                                                    aGATGCAAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCa  >  2:753206/1‑90 (MQ=255)
                                                                                          |                                                                                   
ATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCA  >  AC_000091/2684384‑2684556

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