Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,961,605 G→A R80R (CGC→CGT gppA ← guanosine pentaphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,961,6050GA99.5% 765.0 / ‑5.2 217R80R (CGC→CGTgppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (109/107);  total (109/108)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

Only 100 of 225 total aligned reads displayed.

CACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAGGGC  >  NC_000913/3961515‑3961694
                                                                                          |                                                                                         
cACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTACCGTCGCCGTAGCGACAAcacg                                                                                         >  2:3794915/1‑93 (MQ=255)
cACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcacg                                                                                         >  1:853145/1‑93 (MQ=255)
 cccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGa                                                                                        <  2:1540454/92‑1 (MQ=255)
 aCCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGa                                                                                        >  1:2278738/1‑93 (MQ=255)
   cTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAt                                                                                      <  2:2675327/93‑1 (MQ=255)
     gTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAttt                                                                                    >  1:757225/1‑93 (MQ=255)
      tACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTg                                                                                   >  1:4636533/1‑93 (MQ=255)
      tACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTg                                                                                   <  1:1369123/93‑1 (MQ=255)
        ccGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCg                                                                                 >  1:4515690/1‑93 (MQ=255)
         cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCga                                                                                <  2:1594255/93‑1 (MQ=255)
         cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCga                                                                                >  1:2209248/1‑93 (MQ=255)
         cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCCCAGTAGCGACAACACGAATTTGCga                                                                                <  2:543040/93‑1 (MQ=255)
          ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgag                                                                               <  2:118557/93‑1 (MQ=255)
          ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgag                                                                               <  2:409308/93‑1 (MQ=255)
          ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCATAACGTCGCCGTAGCGACAACACGAATTTGCgag                                                                               <  2:443360/93‑1 (MQ=255)
          cgACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgag                                                                               <  1:3034597/92‑1 (MQ=255)
           gACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagg                                                                              >  1:4339258/1‑93 (MQ=255)
            aCAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgaggg                                                                             <  1:1838614/93‑1 (MQ=255)
              aaCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggag                                                                           <  1:1219184/93‑1 (MQ=255)
              aaCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggag                                                                           <  1:2498736/93‑1 (MQ=255)
              aaCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggag                                                                           >  2:123988/1‑93 (MQ=255)
              aaCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggag                                                                           >  1:3685626/1‑93 (MQ=255)
                ccGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGggagga                                                                         <  1:2387083/93‑1 (MQ=255)
                 cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGaa                                                                        >  2:3081820/1‑93 (MQ=255)
                 cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACAAATTTGCGAGGGAGGaa                                                                        >  2:3234032/1‑93 (MQ=255)
                  gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAat                                                                       <  1:843975/93‑1 (MQ=255)
                   agaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAata                                                                      <  1:2808534/93‑1 (MQ=255)
                     aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATc                                                                    <  2:2119289/93‑1 (MQ=255)
                      aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCt                                                                   <  2:4148179/93‑1 (MQ=255)
                      aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCt                                                                   <  1:3403459/93‑1 (MQ=255)
                       tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCtt                                                                  >  1:64850/1‑93 (MQ=255)
                       tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCtt                                                                  >  2:167072/1‑93 (MQ=255)
                         tCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGc                                                                >  1:1583459/1‑93 (MQ=255)
                         tCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGc                                                                <  1:2972650/93‑1 (MQ=255)
                          ccTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCa                                                               <  1:2012903/93‑1 (MQ=255)
                           cTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAg                                                              <  1:2852773/93‑1 (MQ=255)
                           cTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAg                                                              >  2:659499/1‑93 (MQ=255)
                             gcgcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGAc                                                            <  2:4300756/93‑1 (MQ=255)
                              cgcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACg                                                           <  2:353120/93‑1 (MQ=255)
                               gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGt                                                          >  1:324239/1‑93 (MQ=255)
                                 tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                                                        >  1:4189196/1‑93 (MQ=255)
                                   tGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAg                                                      <  2:2022472/93‑1 (MQ=255)
                                    ggCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGc                                                     <  2:4494225/93‑1 (MQ=255)
                                     gCAATAAAATCCCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCa                                                    <  1:2218583/93‑1 (MQ=255)
                                      cAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaa                                                   >  1:1935696/1‑93 (MQ=255)
                                      cAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaa                                                   >  1:2743697/1‑93 (MQ=255)
                                       attAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaaa                                                  >  2:53781/3‑93 (MQ=255)
                                       aaTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaga                                                  >  1:3831115/1‑91 (MQ=255)
                                         tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                                                <  2:1197691/93‑1 (MQ=255)
                                         tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                                                <  1:307971/93‑1 (MQ=255)
                                         aaaaaACGCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                                                <  2:1305276/92‑1 (MQ=255)
                                          aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGCGGAATATCTTGCAGACGTTCAGCAAACAg                                               >  1:1309316/1‑93 (MQ=255)
                                          aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg                                               >  1:2266673/1‑93 (MQ=255)
                                          aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg                                               <  2:4015159/93‑1 (MQ=255)
                                          aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg                                               <  2:3674067/93‑1 (MQ=255)
                                          aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg                                               <  1:3272987/93‑1 (MQ=255)
                                           aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg                                              >  2:18410/1‑93 (MQ=255)
                                            aaTCCCCCGCATTGACCGCAAGGCGTAACGTCTCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                                             <  1:1477963/93‑1 (MQ=255)
                                            aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                                             <  1:3123636/93‑1 (MQ=255)
                                                aCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaa                                         >  2:423222/1‑93 (MQ=255)
                                                aCCCGCATTGACGGAAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaa                                         >  1:1704887/1‑93 (MQ=255)
                                                 cccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                                        >  2:798843/1‑93 (MQ=255)
                                                 cccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                                        <  2:299348/93‑1 (MQ=255)
                                                  ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAAc                                       >  1:4434559/1‑93 (MQ=255)
                                                  ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAAc                                       >  1:3616542/1‑93 (MQ=255)
                                                   cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                                      <  1:3464775/93‑1 (MQ=255)
                                                   cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                                      <  1:3800808/93‑1 (MQ=255)
                                                    gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                                     >  1:3502067/1‑93 (MQ=255)
                                                    gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                                     >  2:3937708/1‑93 (MQ=255)
                                                    gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                                     >  2:2860356/1‑93 (MQ=255)
                                                     cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAtt                                    <  1:625884/93‑1 (MQ=255)
                                                     cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAtt                                    >  2:3539768/1‑93 (MQ=255)
                                                       ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGc                                  >  1:3955021/1‑93 (MQ=255)
                                                          aCGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCaa                               <  2:3309208/93‑1 (MQ=255)
                                                           cGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAGACATTGCCAAc                              <  1:4005435/93‑1 (MQ=255)
                                                           cGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAc                              <  1:3049796/93‑1 (MQ=255)
                                                           cGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAc                              >  1:2935351/1‑93 (MQ=255)
                                                            ggCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTCGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAcc                             >  1:1769997/1‑93 (MQ=255)
                                                             gCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcg                            >  1:2408448/1‑93 (MQ=255)
                                                               aaGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcgcg                          >  2:3026962/1‑93 (MQ=255)
                                                                aGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc                         <  2:2536646/93‑1 (MQ=255)
                                                                   cGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTcc                      <  1:1823769/93‑1 (MQ=255)
                                                                   cGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTcc                      >  2:3765613/1‑93 (MQ=255)
                                                                     tAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt                    >  1:4099879/1‑93 (MQ=255)
                                                                       acgtCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTg                  >  1:24736/4‑93 (MQ=255)
                                                                         gtCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCt                <  2:1094495/92‑1 (MQ=255)
                                                                          tCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCtt               >  2:3582456/1‑93 (MQ=255)
                                                                            gCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAGACATTGCCAACCGCGCTCCATTGCTTCa             >  2:2286200/1‑93 (MQ=255)
                                                                            gCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCa             >  2:3173412/1‑93 (MQ=255)
                                                                              cGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCAtt           <  1:973364/93‑1 (MQ=255)
                                                                                tAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTgg         <  2:2638314/93‑1 (MQ=255)
                                                                                tAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTgg         >  1:4267686/1‑93 (MQ=255)
                                                                                 aGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGga        <  2:3347241/93‑1 (MQ=255)
                                                                                    gACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGagag     <  2:2072662/93‑1 (MQ=255)
                                                                                    gACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGagag     >  2:1679509/1‑93 (MQ=255)
                                                                                     acaacaCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAgg    >  2:3902262/1‑93 (MQ=255)
                                                                                     acaacaCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAgg    >  1:2785001/1‑93 (MQ=255)
                                                                                      cAACACGAATTTGCGAGGGAGGAATATCTTGCAGAAGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAggg   >  2:4664762/1‑93 (MQ=255)
                                                                                       aaCACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAggga  >  1:1319874/1‑92 (MQ=255)
                                                                                       aaCACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAGGGc  >  1:3726145/1‑93 (MQ=255)
                                                                                          |                                                                                         
CACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAGGGC  >  NC_000913/3961515‑3961694

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