Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 772,731 A→G Y156C (TAC→TGC)  cydB → cytochrome d terminal oxidase polypeptide subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913772,7310AG99.6% 865.9 / NA 270Y156C (TAC→TGC) cydBcytochrome d terminal oxidase polypeptide subunit II
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (147/121);  total (147/123)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 283 total aligned reads displayed.

GTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAA  >  NC_000913/772641‑772820
                                                                                          |                                                                                         
gTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGtgct                                                                                         >  2:4053404/1‑93 (MQ=255)
   gCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTaca                                                                                      >  2:4987283/1‑93 (MQ=255)
   gCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTaca                                                                                      >  1:6177575/1‑93 (MQ=255)
    cTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTacac                                                                                     >  1:5910426/1‑93 (MQ=255)
     ttCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTacacg                                                                                    >  1:806973/1‑92 (MQ=255)
     ttCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACAcc                                                                                    >  1:991593/1‑93 (MQ=255)
     ttCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACAcc                                                                                    >  1:2422887/1‑93 (MQ=255)
         ttCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTa                                                                                >  1:2200160/1‑93 (MQ=255)
          tCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTaa                                                                               <  2:431418/93‑1 (MQ=255)
            cgccgcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAAct                                                                             <  1:2331366/93‑1 (MQ=255)
            cgccgcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAAct                                                                             <  2:5632535/93‑1 (MQ=255)
             gccgcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAActt                                                                            <  1:131929/93‑1 (MQ=255)
              cggcTGGTAAGTGGTGTAGCGTTCGGTAACCGGTTTCAGGGCGTACCGGTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAActtg                                                                           <  1:348796/91‑2 (MQ=255)
              ccgcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAActtc                                                                           <  1:1765146/93‑1 (MQ=255)
                gcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAActtctt                                                                         <  2:5484592/93‑1 (MQ=255)
                 cTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACttcttc                                                                        >  1:2924600/1‑93 (MQ=255)
                  tGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTcc                                                                       >  2:469547/1‑93 (MQ=255)
                     tAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGt                                                                    <  1:6708600/93‑1 (MQ=255)
                      taTTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGtt                                                                   <  1:643345/92‑1 (MQ=255)
                       aTTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCATTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTg                                                                  >  2:1455189/1‑93 (MQ=255)
                        ttGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGc                                                                 <  1:833309/93‑1 (MQ=255)
                        ttGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGc                                                                 <  1:493825/93‑1 (MQ=255)
                         tGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCt                                                                >  2:7789733/1‑93 (MQ=255)
                           gtgtAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTa                                                              >  1:3139479/1‑93 (MQ=255)
                            tgtAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTaa                                                             >  2:3649639/1‑93 (MQ=255)
                             gtAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAAc                                                            <  2:565387/93‑1 (MQ=255)
                             gtAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAAc                                                            <  2:3990786/93‑1 (MQ=255)
                                gCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCg                                                         <  2:7302197/93‑1 (MQ=255)
                                  gTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGtt                                                       <  2:1702230/93‑1 (MQ=255)
                                  gTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGtt                                                       >  2:3283168/1‑93 (MQ=255)
                                  gTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGtt                                                       <  1:5986855/93‑1 (MQ=255)
                                  gTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGtt                                                       >  1:170156/1‑93 (MQ=255)
                                  gTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGtt                                                       >  2:6931806/1‑93 (MQ=255)
                                  gTTCGGTAACCTGGTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTCACTTCTTCCAGTTGCTTAACCCGtc                                                       >  2:1312383/1‑92 (MQ=255)
                                   ttCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTc                                                      <  2:4183483/93‑1 (MQ=255)
                                   ttCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTc                                                      >  1:3547890/1‑93 (MQ=255)
                                    tCGGTAACCTGTTGCAGTGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCg                                                     <  1:4168102/93‑1 (MQ=255)
                                    tCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCg                                                     >  2:6057909/1‑93 (MQ=255)
                                      ggTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGc                                                   >  1:5592948/1‑93 (MQ=255)
                                      ggTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGc                                                   >  2:4604389/1‑93 (MQ=255)
                                       gTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGcc                                                  >  2:3317464/1‑93 (MQ=255)
                                       gTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGcc                                                  <  2:7871536/93‑1 (MQ=255)
                                        tAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCct                                                 >  2:6726028/1‑93 (MQ=255)
                                        gaaCCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCct                                                 <  2:3029405/92‑1 (MQ=255)
                                          aCCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCctgc                                               <  2:4894218/93‑1 (MQ=255)
                                           ccTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCctgct                                              <  2:1765607/93‑1 (MQ=255)
                                            cTGTTGCAGGGCGTGCCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCctgctg                                             <  2:185822/93‑1 (MQ=255)
                                            cTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCctgctg                                             >  1:4434831/1‑93 (MQ=255)
                                              gTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGc                                           >  2:2189282/1‑93 (MQ=255)
                                               ttGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCa                                          <  1:2512509/93‑1 (MQ=255)
                                               ttGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCa                                          <  2:5015642/93‑1 (MQ=255)
                                               gtGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCCACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCa                                          >  2:7034855/2‑93 (MQ=255)
                                                 gCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAgg                                        <  1:122152/93‑1 (MQ=255)
                                                 gCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAgg                                        >  1:1946799/1‑93 (MQ=255)
                                                    gtgCGTTCCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgt                                     <  1:4292538/91‑1 (MQ=255)
                                                    gggCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgt                                     <  1:5204909/93‑1 (MQ=255)
                                                     ggCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtg                                    >  1:2608221/1‑93 (MQ=255)
                                                      gCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtgg                                   <  2:2091665/93‑1 (MQ=255)
                                                      gCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtgg                                   >  2:7495183/1‑93 (MQ=255)
                                                       cGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggt                                  <  2:1587789/93‑1 (MQ=255)
                                                       cGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggt                                  <  1:3332682/93‑1 (MQ=255)
                                                       cGTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggt                                  <  2:2939992/93‑1 (MQ=255)
                                                        gTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggtg                                 >  2:2812287/1‑93 (MQ=255)
                                                        gTACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggtg                                 >  1:5161066/1‑93 (MQ=255)
                                                         tACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggtgg                                >  2:5865463/1‑92 (MQ=255)
                                                         tACCGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGa                                <  1:326386/93‑1 (MQ=255)
                                                           ccGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGc                              >  2:5718426/1‑93 (MQ=255)
                                                           ccGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGc                              >  1:3885957/1‑93 (MQ=255)
                                                            cGTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCg                             >  2:7415077/1‑93 (MQ=255)
                                                             gTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGt                            >  1:1031000/1‑93 (MQ=255)
                                                             gTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGt                            >  1:7357023/1‑93 (MQ=255)
                                                             gTTCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGt                            >  2:5413835/1‑93 (MQ=255)
                                                             gTTCAACGTTGATGAATATCTGCGTCTGGGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTTGCAGGCGTGGTGTGCGt                            >  2:1654763/1‑93 (MQ=255)
                                                              ttCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTa                           <  1:6110874/93‑1 (MQ=255)
                                                              ttCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTa                           <  2:6442799/93‑1 (MQ=255)
                                                              ttCAACGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTa                           <  2:628379/93‑1 (MQ=255)
                                                                 aaCGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGGTCGGCCAGCTGGCAGGAGAGGGGAGCGTAggg                        >  2:6556913/1‑93 (MQ=255)
                                                                   cGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgat                      <  1:2242684/93‑1 (MQ=255)
                                                                   cGTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgat                      <  2:3434741/93‑1 (MQ=255)
                                                                    gTTGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgatg                     <  1:5445702/93‑1 (MQ=255)
                                                                     ttGATGAATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgatga                    >  1:2727332/1‑93 (MQ=255)
                                                                       gatgaATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATc                  >  2:5888492/1‑93 (MQ=255)
                                                                       gatgaATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATc                  >  1:6979856/1‑93 (MQ=255)
                                                                        atgaATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCa                 <  1:6531264/93‑1 (MQ=255)
                                                                         tgaATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCAt                >  2:2140953/1‑93 (MQ=255)
                                                                          gaATATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTGGGGATGATCAtt               >  1:5068403/1‑93 (MQ=255)
                                                                           aaTATCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTa              <  1:1860072/93‑1 (MQ=255)
                                                                            atatCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTAc             <  1:4024573/93‑1 (MQ=255)
                                                                             tatCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGGGCGTAGGGATGATCATTACt            >  1:1922719/1‑93 (MQ=255)
                                                                             tatCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACt            <  1:4715357/93‑1 (MQ=255)
                                                                             tatCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGGAGGGATGATCATTACt            >  2:2510240/1‑93 (MQ=255)
                                                                              atCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTc           <  2:93819/93‑1 (MQ=255)
                                                                              atCTGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTc           <  1:4987283/93‑1 (MQ=255)
                                                                               tctcCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACACGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCa          <  2:7226429/93‑1 (MQ=255)
                                                                                 tGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATGATTACTCAgg        >  1:4020215/1‑93 (MQ=255)
                                                                                 tGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCCTTACTCAgg        >  1:366467/1‑93 (MQ=255)
                                                                                 tGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAgg        <  2:5247173/93‑1 (MQ=255)
                                                                                 tGCGTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAgg        >  1:7094455/1‑93 (MQ=255)
                                                                                    gTCTGTGCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGCTGATCCTTACTCATGgcg     >  2:4149429/1‑93 (MQ=255)
                                                                                       tgtgCTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCaa  >  1:7588313/1‑93 (MQ=255)
                                                                                          |                                                                                         
GTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAA  >  NC_000913/772641‑772820

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