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.3% 1014.7 / ‑5.8 284R80R (CGC→CGT) gppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (150/130);  total (150/134)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 293 total aligned reads displayed.

CACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAGGGC  >  NC_000913/3961515‑3961694
                                                                                          |                                                                                         
cACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcacg                                                                                         >  2:4631677/1‑93 (MQ=255)
cACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcacg                                                                                         >  2:1858308/1‑93 (MQ=255)
 aCCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGa                                                                                        >  2:2725053/1‑93 (MQ=255)
  ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGaa                                                                                       <  2:4482535/93‑1 (MQ=255)
  ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGaa                                                                                       <  2:4310698/93‑1 (MQ=255)
  ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGaa                                                                                       <  2:3851943/93‑1 (MQ=255)
  ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGaa                                                                                       >  2:4148505/1‑93 (MQ=255)
  ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATAACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGaa                                                                                       >  1:3240684/1‑93 (MQ=255)
   cTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAt                                                                                      >  2:2944025/1‑93 (MQ=255)
   cTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAt                                                                                      <  2:946674/93‑1 (MQ=255)
    tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAtt                                                                                     <  1:4285967/93‑1 (MQ=255)
    tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAtt                                                                                     <  2:2365705/93‑1 (MQ=255)
    tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAtt                                                                                     <  2:1676372/93‑1 (MQ=255)
    tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAtt                                                                                     <  1:3383359/93‑1 (MQ=255)
     gTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCCCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAttt                                                                                    <  2:3084011/93‑1 (MQ=255)
     gTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAttt                                                                                    <  2:4266244/93‑1 (MQ=255)
        ccGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCg                                                                                 >  1:2878307/1‑93 (MQ=255)
         cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCga                                                                                <  2:1513018/93‑1 (MQ=255)
         cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCga                                                                                <  2:523688/93‑1 (MQ=255)
          ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgag                                                                               >  1:1220525/1‑93 (MQ=255)
          ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgag                                                                               >  2:3180104/1‑93 (MQ=255)
          ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgag                                                                               <  1:129615/93‑1 (MQ=255)
          ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgag                                                                               <  2:2165732/93‑1 (MQ=255)
             cAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgaggga                                                                            >  2:3777800/1‑93 (MQ=255)
               aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggagg                                                                          <  2:1975730/93‑1 (MQ=255)
                 cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGaa                                                                        >  1:1669117/1‑93 (MQ=255)
                 cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGaa                                                                        >  1:1130298/1‑93 (MQ=255)
                  gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAat                                                                       >  1:3481546/1‑93 (MQ=255)
                  gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAat                                                                       >  2:2644161/1‑93 (MQ=255)
                  gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAat                                                                       >  2:92215/1‑93 (MQ=255)
                  gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAat                                                                       >  2:4546883/1‑93 (MQ=255)
                  gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAat                                                                       >  1:3829499/1‑93 (MQ=255)
                    gaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAatat                                                                     <  1:2372010/93‑1 (MQ=255)
                    gaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAatat                                                                     <  1:1263596/93‑1 (MQ=255)
                     aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATc                                                                    <  2:3733573/93‑1 (MQ=255)
                     aaTTCCCTGCGCTTTGACAATAAAATCGCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATc                                                                    <  2:3508891/93‑1 (MQ=255)
                      aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCt                                                                   <  2:477205/93‑1 (MQ=255)
                      aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCt                                                                   <  1:461728/93‑1 (MQ=255)
                      aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCt                                                                   >  2:180407/1‑93 (MQ=255)
                      aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCt                                                                   >  1:1538491/1‑93 (MQ=255)
                       tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCtt                                                                  >  1:4402318/1‑93 (MQ=255)
                        ttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTg                                                                 <  1:2149255/93‑1 (MQ=255)
                         tCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCATAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGc                                                                <  2:3952693/93‑1 (MQ=255)
                          ccTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCa                                                               <  2:715693/93‑1 (MQ=255)
                           cTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAg                                                              <  1:2203708/93‑1 (MQ=255)
                               gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGt                                                          <  1:3431526/93‑1 (MQ=255)
                               gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCATAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGt                                                          >  1:4652938/1‑93 (MQ=255)
                                cTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGtt                                                         >  1:4455790/1‑93 (MQ=255)
                                 tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                                                        <  1:4677740/93‑1 (MQ=255)
                                 tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                                                        >  1:4467214/1‑93 (MQ=255)
                                 tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGGCGCCGTAGCGACAACCCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                                                        <  2:2605363/93‑1 (MQ=255)
                                  ttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCa                                                       <  2:4681730/93‑1 (MQ=255)
                                       aaTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaaa                                                  >  1:3566048/1‑93 (MQ=255)
                                        aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGGCAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAAc                                                 >  2:3712592/1‑93 (MQ=255)
                                        aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAAc                                                 >  1:1070166/1‑93 (MQ=255)
                                        aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAAc                                                 >  2:2873706/1‑93 (MQ=255)
                                         tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                                                <  2:3036008/93‑1 (MQ=255)
                                         tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                                                >  1:4042426/1‑93 (MQ=255)
                                         tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                                                <  1:3662734/93‑1 (MQ=255)
                                         tAAAATCACCCGCATTGACGGCAAGGCGTAACGCCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                                                >  1:650027/1‑93 (MQ=255)
                                          aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg                                               >  1:4513640/1‑93 (MQ=255)
                                           aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg                                              >  1:481282/1‑93 (MQ=255)
                                             aTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcg                                            >  1:1003650/1‑93 (MQ=255)
                                             aTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcg                                            >  1:3726420/1‑93 (MQ=255)
                                                aCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaa                                         >  1:3549545/1‑93 (MQ=255)
                                                 cccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                                        >  1:1417301/1‑93 (MQ=255)
                                                 cccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                                        >  2:3116911/1‑93 (MQ=255)
                                                     cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAtt                                    <  1:674273/93‑1 (MQ=255)
                                                       ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGc                                  <  1:3022515/93‑1 (MQ=255)
                                                       ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGc                                  >  2:3438619/1‑93 (MQ=255)
                                                       ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGc                                  >  1:2276383/1‑93 (MQ=255)
                                                        tGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGcc                                 >  1:891304/1‑93 (MQ=255)
                                                           cGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAc                              >  2:1136091/1‑93 (MQ=255)
                                                            ggCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAcc                             <  2:869949/93‑1 (MQ=255)
                                                                aGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc                         >  1:866657/1‑93 (MQ=255)
                                                                aGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc                         >  2:4109188/1‑93 (MQ=255)
                                                                  gCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTc                       <  1:1772247/93‑1 (MQ=255)
                                                                    gTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCa                     <  1:2755287/93‑1 (MQ=255)
                                                                    gTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCa                     >  1:542479/1‑93 (MQ=255)
                                                                     tAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt                    >  2:676030/1‑93 (MQ=255)
                                                                     tAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt                    >  2:1660001/1‑93 (MQ=255)
                                                                       acgtCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTg                  >  1:1042367/4‑93 (MQ=255)
                                                                        cgtCGCCGGAGCGACAACACGAATTTGCGAGGCAGGAATATCTTCCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGc                 <  2:1246556/91‑1 (MQ=255)
                                                                          tCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCtt               >  1:1571596/1‑93 (MQ=255)
                                                                          tCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCtt               >  2:2495452/1‑93 (MQ=255)
                                                                            gCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCa             >  1:3093066/1‑93 (MQ=255)
                                                                            gCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCa             <  1:4251112/93‑1 (MQ=255)
                                                                             ccGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCAt            <  2:2800192/93‑1 (MQ=255)
                                                                              cGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCAtt           <  2:1225016/93‑1 (MQ=255)
                                                                              cGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCAtt           >  1:710620/1‑93 (MQ=255)
                                                                               gTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCGGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTg          >  1:2606218/1‑93 (MQ=255)
                                                                                tAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTgg         >  1:2589759/1‑93 (MQ=255)
                                                                                 aGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGga        >  2:2247196/1‑93 (MQ=255)
                                                                                 aGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGga        <  1:171364/93‑1 (MQ=255)
                                                                                  gCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGgag       >  2:1607863/1‑93 (MQ=255)
                                                                                   cGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGgaga      <  1:2852208/93‑1 (MQ=255)
                                                                                    gACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGagag     >  2:4038555/1‑93 (MQ=255)
                                                                                    gACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGagag     <  2:3328/93‑1 (MQ=255)
                                                                                      cAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAggg   <  2:1010019/93‑1 (MQ=255)
                                                                                       aaCACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAGGGc  >  2:1333149/1‑93 (MQ=255)
                                                                                          |                                                                                         
CACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAGGGC  >  NC_000913/3961515‑3961694

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