Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,147,319 A→T 5.3% C74* (TGT→TGA)  alkA ← 3‑methyl‑adenine DNA glycosylase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,147,3190AT 5.3% 235.1 / 3.5 76C74* (TGT→TGA) alkA3‑methyl‑adenine DNA glycosylase II
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base T (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

GCCTAAAATCGCCCGCACGCCCTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTG  >  NC_000913/2147172‑2147466
                                                                                                                                                   |                                                                                                                                                   
gCCTAAAATCGCCCGCACGCCCTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAACCAGGCGGCTCATTTTCGCCAGa                                                                                                                                                     >  3:64642/1‑148 (MQ=255)
   tAAAATCGCCCGCACGCCCTGCTCCAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACAct                                                                                                                                                 >  6:365518/1‑149 (MQ=255)
    aaaaTCGCCCGCACGCCCTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACActc                                                                                                                                                >  4:83433/1‑149 (MQ=255)
       aTCGCCCGCACGCCCTGCTCAAAAGCATCAACACAGCCGGGTACACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGGGGGGTACATTGCAGATCAAACAGGCGGCGCATTTTCGCCAGACACTCTGc                                                                                                                                             <  8:433263/149‑1 (MQ=255)
        tCGCCCGCACGCCCTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCg                                                                                                                                            >  2:256974/1‑149 (MQ=255)
        tCGCCCGCACGCCCTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCg                                                                                                                                            >  2:368238/1‑149 (MQ=255)
        tCGCCCGCAAGCCCTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCg                                                                                                                                            <  8:389441/149‑1 (MQ=255)
          gCCCGCACGCCCTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGc                                                                                                                                          <  2:230334/149‑1 (MQ=255)
                aCGCCCTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAgg                                                                                                                                    <  1:99114/149‑1 (MQ=255)
                 cGCCCTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGt                                                                                                                                   <  1:368238/149‑1 (MQ=255)
                 cGCCCTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGt                                                                                                                                   >  1:389225/1‑149 (MQ=255)
                  gCCCTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGtt                                                                                                                                  <  1:256974/149‑1 (MQ=255)
                  gCCCTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGtt                                                                                                                                  <  4:83430/149‑1 (MQ=255)
                    ccTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTc                                                                                                                                 >  1:296137/1‑148 (MQ=255)
                      tGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTa                                                                                                                               >  4:123659/1‑148 (MQ=255)
                           aaaaGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCt                                                                                                                          <  1:76157/148‑1 (MQ=255)
                           aaaaGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTg                                                                                                                         >  6:388020/1‑149 (MQ=255)
                                  tCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTa                                                                                                                   >  4:295574/1‑148 (MQ=255)
                                    aaCACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAAt                                                                                                                <  6:367229/149‑1 (MQ=255)
                                     acacAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAAtt                                                                                                               >  5:330468/1‑149 (MQ=255)
                                        caGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTAt                                                                                                            <  8:180191/149‑1 (MQ=255)
                                           ccGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTAtgtg                                                                                                         <  8:106597/149‑1 (MQ=255)
                                               gTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCaga                                                                                                     <  3:123659/149‑1 (MQ=255)
                                                      cAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATgg                                                                                               <  7:43708/148‑1 (MQ=255)
                                                      cAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATgg                                                                                               >  8:43708/1‑148 (MQ=255)
                                                      cAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGc                                                                                              >  8:245517/1‑149 (MQ=255)
                                                      cAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGc                                                                                              <  2:389225/149‑1 (MQ=255)
                                                        aTCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCgg                                                                                            <  3:383696/149‑1 (MQ=255)
                                                                 gcgccgcgccTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTCCAGAGTATGGCGGGCTATATcc                                                                                   >  4:374426/1‑149 (MQ=255)
                                                                   gccgcgccTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACAGTGCAGATAAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCgg                                                                                 <  8:413947/149‑1 (MQ=255)
                                                                     cgcgccTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGaa                                                                               <  5:365518/149‑1 (MQ=255)
                                                                      gcgccTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAAt                                                                              <  4:365370/149‑1 (MQ=255)
                                                                         ccTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGCCACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGc                                                                           >  1:423058/1‑149 (MQ=255)
                                                                         ccTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAg                                                                            >  1:277436/1‑148 (MQ=255)
                                                                         ccTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAg                                                                            >  7:421363/1‑148 (MQ=255)
                                                                          cTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCaga                                                                                                     <  6:109359/122‑1 (MQ=255)
                                                                          cTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCaga                                                                                                     >  5:109359/1‑122 (MQ=255)
                                                                          cTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAg                                                                            >  8:32830/1‑147 (MQ=255)
                                                                           tAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAg                                                                         >  5:139403/1‑149 (MQ=255)
                                                                             aCCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAAAAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTc                                                                       >  4:142417/1‑149 (MQ=255)
                                                                              ccTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTca                                                                      >  2:371792/1‑149 (MQ=255)
                                                                                  cccAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTcaccac                                                                  <  6:330468/149‑1 (MQ=255)
                                                                                         cACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGa                                                           <  2:296137/149‑1 (MQ=255)
                                                                                            cGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGa                                                                                <  7:28744/125‑1 (MQ=255)
                                                                                            cGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGa                                                                                >  8:28744/1‑125 (MQ=255)
                                                                                             gTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATAtt                                                       <  7:245517/149‑1 (MQ=255)
                                                                                              ttAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTc                                                      >  3:313625/1‑149 (MQ=255)
                                                                                               tAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCAcgc                                                                >  3:151092/1‑138 (MQ=255)
                                                                                               tAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCAcgc                                                                <  4:151092/138‑1 (MQ=255)
                                                                                               tAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCg                                                     >  1:263406/1‑149 (MQ=255)
                                                                                                aaCAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTACACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCg                                                     <  1:371792/148‑1 (MQ=255)
                                                                                                 aCAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGa                                                           <  5:57351/141‑1 (MQ=255)
                                                                                                 aCAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGa                                                           >  6:57351/1‑141 (MQ=255)
                                                                                                 aCAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGcc                                                   <  2:263406/149‑1 (MQ=255)
                                                                                                 aCAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACAACGCCGCGATATTCGcc                                                   >  1:427662/1‑149 (MQ=255)
                                                                                                  cAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGcc                                                   >  1:195301/1‑148 (MQ=255)
                                                                                                            gTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTc                                                                       >  1:266713/1‑118 (MQ=255)
                                                                                                            gTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTc                                                                       <  2:266713/118‑1 (MQ=255)
                                                                                                                  ttGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACggg                                  <  7:32830/149‑1 (MQ=255)
                                                                                                                    gCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCa                                >  8:403213/1‑149 (MQ=255)
                                                                                                                     cAGATCAAACAGGCGGCTCATTTTCGCCAGTCACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCa                                           <  7:125109/137‑1 (MQ=255)
                                                                                                                     cAGATCAAACAGGCGGCTCATTTTCGCCAGTCACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCa                                           <  7:131191/137‑1 (MQ=255)
                                                                                                                     cAGATCAAACAGGCGGCTCATTTTCGCCAGTCACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCa                                           >  8:125109/1‑137 (MQ=255)
                                                                                                                     cAGATCAAACAGGCGGCTCATTTTCGCCAGTCACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCa                                           >  8:131191/1‑137 (MQ=255)
                                                                                                                     cAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCa                                <  3:295574/148‑1 (MQ=255)
                                                                                                                       gATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCataa                             >  5:384154/1‑149 (MQ=255)
                                                                                                                           acaCAGGCGGCTAATTTTCTCCAGACACTCTGCGGCAACAGGTTCTAAACCTTCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCAtaataa                          <  4:64642/146‑1 (MQ=255)
                                                                                                                           aaaCAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCAtaataa                          <  5:366440/148‑1 (MQ=255)
                                                                                                                              cAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGt                      <  5:146194/149‑1 (MQ=255)
                                                                                                                                    gcTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGa                <  3:83433/149‑1 (MQ=255)
                                                                                                                                      tCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAAATTTCCGCGAcc              >  6:248978/1‑149 (MQ=255)
                                                                                                                                       cATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCg             >  2:426398/1‑149 (MQ=255)
                                                                                                                                         ttttCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCataa                             <  2:129778/131‑1 (MQ=255)
                                                                                                                                         ttttCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCataa                             >  1:129778/1‑131 (MQ=255)
                                                                                                                                           ttCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCCCGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTc         >  6:156928/1‑149 (MQ=255)
                                                                                                                                                cAGACACTCTGCGGAAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCacac    >  3:338286/1‑149 (MQ=255)
                                                                                                                                                  gACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTg  <  6:384154/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCCTAAAATCGCCCGCACGCCCTGCTCAAAAGCATCAACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTG  >  NC_000913/2147172‑2147466

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