Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,950,330 G→A 44.8% M363I (ATG→ATA argA → amino acid N‑acetyltransferase and inactive acetylglutamate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,950,3300GA44.8% 29.8 / 105.0 89M363I (ATG→ATAargAamino acid N‑acetyltransferase and inactive acetylglutamate kinase
Reads supporting (aligned to +/- strand):  ref base G (18/31);  new base A (18/22);  total (36/53)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

GGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTTAAGCA  >  NC_000913/2950231‑2950427
                                                                                                   |                                                                                                 
ggAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATg                                                                                                   <  1:675170/100‑1 (MQ=255)
  agaTGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGc                                                                                                 <  1:1019082/100‑1 (MQ=255)
   gaTGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGcc                                                                                                <  1:354875/100‑1 (MQ=255)
   gaTGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGcc                                                                                                <  2:338435/100‑1 (MQ=255)
    ctGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCt                                                                                               <  2:1882304/99‑1 (MQ=255)
     tGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCtg                                                                                              >  2:413495/1‑100 (MQ=255)
      ggAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCtag                                                                                             >  1:1882294/1‑98 (MQ=255)
      ggAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCtgt                                                                                             >  2:1951713/1‑100 (MQ=255)
       naaaTCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCtgtg                                                                                            >  2:904834/2‑100 (MQ=255)
       gAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCtgtg                                                                                            <  1:1932313/100‑1 (MQ=255)
         aaTCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCtgtgtg                                                                                          <  1:1066321/100‑1 (MQ=255)
         aaTCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCtgtgtg                                                                                          <  1:920920/100‑1 (MQ=255)
           tCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGc                                                                                        <  2:1898926/100‑1 (MQ=255)
             gACAAATTCATCATTATTCAGCGCGATAACACGACTATTGCCTGTGGCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAg                                                                                      <  1:604061/100‑1 (MQ=255)
             gACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTGCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAg                                                                                      <  1:2208291/100‑1 (MQ=255)
                 aaTTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCa                                                                                  >  2:290987/1‑100 (MQ=255)
                     cACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCg                                                                              >  2:1087437/1‑100 (MQ=255)
                     cACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCg                                                                              >  2:1684345/1‑100 (MQ=255)
                       ccATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGa                                                                            >  1:1833382/1‑100 (MQ=255)
                       ccATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGa                                                                            <  2:490572/100‑1 (MQ=255)
                       ccATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGa                                                                            >  2:688370/1‑100 (MQ=255)
                        cATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGAt                                                                           <  2:1839000/100‑1 (MQ=255)
                        cATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCCTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCgg                                                                             >  1:2228227/1‑98 (MQ=255)
                          ttattCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTa                                                                         <  1:1647409/100‑1 (MQ=255)
                             ttCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCg                                                                      <  1:40156/100‑1 (MQ=255)
                               cAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGcn                                                                    <  1:1887049/100‑2 (MQ=255)
                               cAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCa                                                                    >  2:2068555/1‑100 (MQ=255)
                               cAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCa                                                                    <  2:358752/100‑1 (MQ=255)
                                 gcgcGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGt                                                                  >  1:199957/1‑100 (MQ=255)
                                      aTAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTtcatc                                                             <  2:2038171/100‑1 (MQ=255)
                                      aTAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTtcatc                                                             <  1:2010231/100‑1 (MQ=255)
                                       tAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTtcatca                                                            <  2:572774/100‑1 (MQ=255)
                                        aaCACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCaa                                                           <  2:595813/100‑1 (MQ=255)
                                        aaCACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCaa                                                           >  2:2376263/1‑100 (MQ=255)
                                        aaCACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCaa                                                           <  2:2026592/100‑1 (MQ=255)
                                         acacGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAg                                                          >  2:2270610/1‑100 (MQ=255)
                                         acacGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAg                                                          >  2:246165/1‑100 (MQ=255)
                                          cacGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAgg                                                         >  2:653078/1‑100 (MQ=255)
                                          cacGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAgg                                                         <  1:1324191/100‑1 (MQ=255)
                                          cacGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAgg                                                         >  1:1925373/1‑100 (MQ=255)
                                          cacGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAgg                                                         >  2:2300724/1‑100 (MQ=255)
                                          cacGACTATTGCCTGCGCAGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGAGTAGCAGTTCACCCGGATTACCGCAGTTCATCAAgg                                                         >  1:100599/1‑100 (MQ=255)
                                           acGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAggg                                                        <  1:2496500/100‑1 (MQ=255)
                                            cGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGACGAGAAGATTGGGGAAATGGCCTGTGTGGCACTTCAGCCGGATTATCGCAGTTCATCAAgggg                                                       >  1:74682/1‑100 (MQ=255)
                                            cGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAgggg                                                       <  1:353746/100‑1 (MQ=255)
                                               cTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGa                                                    >  1:2126837/1‑100 (MQ=255)
                                                 aTTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCTGGATTACCGCAGTTCATCAAGGGGTGAAg                                                  >  1:1560045/1‑100 (MQ=255)
                                                 aTTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAg                                                  <  1:1425111/100‑1 (MQ=255)
                                                   ngCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGtt                                                >  2:765350/2‑100 (MQ=255)
                                                    gCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTc                                               >  1:67323/1‑100 (MQ=255)
                                                      cTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTctg                                             >  1:738075/1‑100 (MQ=255)
                                                       tGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTctgc                                            <  2:1241387/100‑1 (MQ=255)
                                                       tGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTctgc                                            >  1:1306413/1‑100 (MQ=255)
                                                       tGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTctgc                                            <  2:1833392/100‑1 (MQ=255)
                                                        gcgcCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTctgct                                           >  2:1588963/1‑100 (MQ=255)
                                                         cgccgcGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTctgctg                                          <  2:1673658/100‑1 (MQ=255)
                                                          gccgcGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTgg                                         <  2:811593/100‑1 (MQ=255)
                                                          gccgcGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTgg                                         <  2:811032/100‑1 (MQ=255)
                                                             gcgcTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAAc                                      <  2:2510179/100‑1 (MQ=255)
                                                              cgcTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACg                                     <  1:2418334/100‑1 (MQ=255)
                                                               gcTGTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGGAGTTCATCAAGGGGTGAAGTTCTGCTGGAACg                                     <  1:2278818/99‑1 (MQ=255)
                                                                 tctATCCGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCAt                                  <  1:1300187/100‑1 (MQ=255)
                                                                    aTCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGc                               <  2:2104706/100‑1 (MQ=255)
                                                                      ccGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCg                             <  1:2343335/100‑1 (MQ=255)
                                                                      ccGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCg                             <  2:2112973/100‑1 (MQ=255)
                                                                      ccGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGc                               <  1:547464/98‑1 (MQ=255)
                                                                      ccGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCg                             >  2:893275/1‑100 (MQ=255)
                                                                       cGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGc                            >  2:2153214/1‑100 (MQ=255)
                                                                       cGTTCCCGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGc                            <  2:948881/100‑1 (MQ=255)
                                                                          tcgcGGAAGAGAAGGTTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCTTCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCa                         <  1:1053202/97‑1 (MQ=255)
                                                                           cccGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAg                        >  2:1468192/1‑100 (MQ=255)
                                                                           cccGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAg                        <  2:1202389/100‑1 (MQ=255)
                                                                            ccGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAgg                       >  1:869161/1‑100 (MQ=255)
                                                                            ccGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAgg                       <  1:1320775/100‑1 (MQ=255)
                                                                            ccGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAgg                       <  2:800148/100‑1 (MQ=255)
                                                                            ccGGAAGAGAAGATTGGGGAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAgg                       >  2:1088010/1‑100 (MQ=255)
                                                                                aagagaagaTTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCACGGGGCGAGGTTCTGCTGGAACGCATTGCCGCTCAGGCGaa                   <  1:912778/100‑1 (MQ=255)
                                                                                 agagaagaTTTGGGAAATGGCCTGTGTGGCAGTTCACCCTGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGcta                    <  1:419037/98‑3 (MQ=255)
                                                                                    gaagaTTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCag               >  2:1351360/1‑100 (MQ=255)
                                                                                         ttGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCgg          <  1:723350/100‑1 (MQ=255)
                                                                                          tGGGGAAATAGCCTGGGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGc         >  1:1777004/1‑100 (MQ=255)
                                                                                           tgggAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCt        <  2:1123655/99‑1 (MQ=255)
                                                                                           ggggAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCt        <  1:496946/100‑1 (MQ=255)
                                                                                           ggggAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCt        >  2:1088635/1‑100 (MQ=255)
                                                                                            gggAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCta       >  1:1123651/1‑99 (MQ=255)
                                                                                             ggAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTTa      <  2:1501933/100‑1 (MQ=255)
                                                                                             ggAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTTa      <  2:1560054/100‑1 (MQ=255)
                                                                                             ggAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTTa      <  2:893138/100‑1 (MQ=255)
                                                                                             ggAAATAGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTTa      <  2:199959/100‑1 (MQ=255)
                                                                                                 aTGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTTAAGCa  <  2:2091541/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTTAAGCA  >  NC_000913/2950231‑2950427

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