Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,418,628 A→G A286A (GCA→GCG b3269 → putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,418,6280AG99.6% 893.0 / ‑4.7 277A286A (GCA→GCGb3269putative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (153/123);  total (153/124)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

Only 100 of 280 total aligned reads displayed.

TTTAATTCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCAGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCATTGCG  >  NC_000913/3418538‑3418717
                                                                                          |                                                                                         
tttAATTCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGt                                                                                         <  2:1751647/93‑1 (MQ=255)
 ttAATTCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTg                                                                                        <  1:3439907/93‑1 (MQ=255)
   aaTTCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGcc                                                                                      <  1:1021045/93‑1 (MQ=255)
    aTTCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCt                                                                                     >  1:2136397/1‑93 (MQ=255)
      tCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTa                                                                                   >  1:1068817/1‑93 (MQ=255)
      tCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTa                                                                                   <  1:2995609/93‑1 (MQ=255)
           aaCTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTc                                                                              >  2:4830376/1‑93 (MQ=255)
           aaCTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTc                                                                              <  2:5725710/93‑1 (MQ=255)
            aCTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCa                                                                             >  1:728275/1‑93 (MQ=255)
             cTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCaa                                                                            >  2:2507959/1‑93 (MQ=255)
             cTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCaa                                                                            <  2:3188801/93‑1 (MQ=255)
               ggCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACa                                                                          <  1:2640157/93‑1 (MQ=255)
               ggCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACa                                                                          >  1:5188852/1‑93 (MQ=255)
               ggCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACa                                                                          >  2:4954155/1‑93 (MQ=255)
                 cAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATg                                                                        <  2:1869179/93‑1 (MQ=255)
                  aGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGa                                                                       >  2:2963974/1‑93 (MQ=255)
                  aGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGa                                                                       <  1:4906405/93‑1 (MQ=255)
                  aGCCTTAACGCTGGCACTTGCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGa                                                                       <  1:4989259/93‑1 (MQ=255)
                   gCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGaa                                                                      >  2:3430403/1‑93 (MQ=255)
                    ccTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAg                                                                     >  1:2441820/1‑93 (MQ=255)
                     cTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGc                                                                    <  1:4702953/93‑1 (MQ=255)
                     cTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGc                                                                    <  1:3706707/93‑1 (MQ=255)
                      ttAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCg                                                                   >  1:831485/1‑93 (MQ=255)
                       tAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCgg                                                                  >  1:6210404/1‑93 (MQ=255)
                        aaCGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGc                                                                 >  1:6317039/1‑93 (MQ=255)
                        aaCGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGc                                                                 >  2:2236600/1‑93 (MQ=255)
                        aaCGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGc                                                                 <  1:3650098/93‑1 (MQ=255)
                        aaCGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCAAGGCGGTGCCTTATGGTCAACATGAAGCGGc                                                                 <  1:1136094/93‑1 (MQ=255)
                         aCGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTGTCCGCGCTGGGGTCCAGGCGGTGCCTTATGGTCAATATGAAGCGGCt                                                                >  1:3313704/1‑93 (MQ=255)
                         aCGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCt                                                                <  2:2611820/93‑1 (MQ=255)
                               gCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCa                                                          <  1:3237117/93‑1 (MQ=255)
                                cACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCAc                                                         <  1:2550626/93‑1 (MQ=255)
                                 aCTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGCCACt                                                        <  2:5894107/93‑1 (MQ=255)
                                  cTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTg                                                       >  1:394417/1‑93 (MQ=255)
                                   tttCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTgg                                                      <  2:6304187/93‑1 (MQ=255)
                                     tCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGa                                                    <  2:3161845/93‑1 (MQ=255)
                                        gTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTa                                                 >  2:5141244/1‑93 (MQ=255)
                                          ttATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTAcc                                               >  2:3697790/1‑93 (MQ=255)
                                          ttATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTAcc                                               >  2:4815628/1‑93 (MQ=255)
                                           tataCATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTTCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTAccc                                              <  2:3576024/93‑1 (MQ=255)
                                           tataCATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTAccc                                              <  2:1801178/93‑1 (MQ=255)
                                             taCATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGGTTACCCaa                                            >  1:2241136/1‑93 (MQ=255)
                                             taCATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCaa                                            <  2:1158480/93‑1 (MQ=255)
                                              aCATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAAc                                           <  2:1004503/93‑1 (MQ=255)
                                               cATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTTCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAAcc                                          <  1:2349788/93‑1 (MQ=255)
                                                aTCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAAccc                                         >  1:561473/1‑93 (MQ=255)
                                                aTCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAAccc                                         >  2:5834714/1‑93 (MQ=255)
                                                aTCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAAccc                                         >  2:4653593/1‑93 (MQ=255)
                                                   tGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGt                                      <  2:323932/93‑1 (MQ=255)
                                                    gCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGtt                                     >  2:6021325/1‑93 (MQ=255)
                                                    gCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGtt                                     <  2:2378586/93‑1 (MQ=255)
                                                    gCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGtt                                     <  2:5060340/93‑1 (MQ=255)
                                                     cATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTta                                    >  2:5590013/1‑93 (MQ=255)
                                                       tttATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacg                                  >  2:1650231/1‑93 (MQ=255)
                                                        ttATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtccgc                                 >  1:1991022/1‑93 (MQ=255)
                                                        ttATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgc                                 <  1:2958067/93‑1 (MQ=255)
                                                         tATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgct                                >  2:759132/1‑93 (MQ=255)
                                                           tCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgctac                              >  1:1681530/1‑93 (MQ=255)
                                                           tCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgctac                              <  1:6103703/93‑1 (MQ=255)
                                                             gCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgctacgc                            >  1:3561784/1‑93 (MQ=255)
                                                              ccGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGcc                           <  2:4352430/93‑1 (MQ=255)
                                                              ccGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGcc                           <  2:3778508/93‑1 (MQ=255)
                                                              ccGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGCTACGCTACGcc                           >  1:809159/1‑93 (MQ=255)
                                                               cGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCa                          >  1:4481862/1‑93 (MQ=255)
                                                                gagaTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGTATTACCCAACCCGGTTACGCTACGCCAg                         >  1:1316623/1‑93 (MQ=255)
                                                                gagaTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAg                         >  2:461897/1‑93 (MQ=255)
                                                                 agaTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAgg                        >  2:3602887/1‑93 (MQ=255)
                                                                  gaTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTAGCCAACCCGGTTACGCTACGGCAGGt                       >  2:4291350/1‑93 (MQ=255)
                                                                  gaTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGt                       >  2:105235/1‑93 (MQ=255)
                                                                   aTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAgggc                      >  2:3062883/1‑91 (MQ=255)
                                                                   aTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTc                      >  1:5758556/1‑93 (MQ=255)
                                                                     tATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCat                    <  1:6705369/93‑1 (MQ=255)
                                                                     tATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCat                    <  2:6430563/93‑1 (MQ=255)
                                                                     tATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCat                    <  1:4020574/93‑1 (MQ=255)
                                                                       tCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCTGTCACTGGGATTACCCAACCCTGTTACGCTACGGCAGGTCatta                  >  1:3359532/1‑93 (MQ=255)
                                                                       tCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCatta                  <  2:3861065/93‑1 (MQ=255)
                                                                       tCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCatta                  <  2:6183715/93‑1 (MQ=255)
                                                                          gcgcTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTc               <  1:1826399/93‑1 (MQ=255)
                                                                          gcgcTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGGTACGCTACGCCAGGGCATTATTc               >  2:5426007/1‑93 (MQ=255)
                                                                           cgcTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCattattgc              <  2:3140912/93‑3 (MQ=255)
                                                                             cTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGCTAAGCTACGCCAGGTCATTATTcccc            >  2:4075611/1‑93 (MQ=255)
                                                                                ggATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAgg         >  2:656502/1‑93 (MQ=255)
                                                                                ggATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAgg         <  1:235642/93‑1 (MQ=255)
                                                                                 gATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGc        >  1:961822/1‑93 (MQ=255)
                                                                                  aTCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCa       >  1:1515225/1‑93 (MQ=255)
                                                                                   tCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCAt      <  1:4464073/93‑1 (MQ=255)
                                                                                   tCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCAt      >  2:2903602/1‑93 (MQ=255)
                                                                                   tCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCAt      <  2:3980571/93‑1 (MQ=255)
                                                                                    ccAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCAtt     <  1:3874779/93‑1 (MQ=255)
                                                                                    ccAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCAtt     >  1:378262/1‑93 (MQ=255)
                                                                                     cAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCATTg    >  2:4754841/1‑93 (MQ=255)
                                                                                     cAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGCTACGCTACGCCAGGTCATTATTCCCCAGGCATTg    >  1:5899037/1‑93 (MQ=255)
                                                                                      aGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCATTGc   <  2:6138939/93‑1 (MQ=255)
                                                                                      aGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCATTGc   <  1:686267/93‑1 (MQ=255)
                                                                                      aGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCATTGc   >  1:6009394/1‑93 (MQ=255)
                                                                                       ggCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCATTGCg  >  1:3160523/1‑93 (MQ=255)
                                                                                       ggCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCATTGCg  >  2:414877/1‑93 (MQ=255)
                                                                                       ggCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCATTGCg  <  2:1263034/93‑1 (MQ=255)
                                                                                       ggCGGTGCCTTATGGTCAACATGAAGCGGCTCCGGCACTGGGATTACCCAACCCGGTTAAGCTACGCCAGGTCATTATTCCCCAGGCATTGCg  >  1:4125679/1‑93 (MQ=255)
                                                                                       ggCGGTGCCTTATGGGCAACATGAAGCGGCTCGGTCACTAGGATTACCCAACCCCGTTACGCTACGCCAGGTCATTATTCCCCAGGCATTGCg  >  2:1969535/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTTAATTCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCAGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCATTGCG  >  NC_000913/3418538‑3418717

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