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.6% 1878.3 / ‑5.1 517R80R (CGC→CGTgppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (272/243);  total (274/243)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

Only 100 of 540 total aligned reads displayed.

CACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAGGGC  >  NC_000913/3961515‑3961694
                                                                                          |                                                                                         
 aCCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGa                                                                                        >  1:2884185/1‑93 (MQ=255)
  ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGaa                                                                                       <  2:6430807/93‑1 (MQ=255)
    tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAtt                                                                                     <  1:7247668/93‑1 (MQ=255)
      tACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTg                                                                                   >  2:3876669/1‑93 (MQ=255)
        ccGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCg                                                                                 <  1:3708672/93‑1 (MQ=255)
         cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCga                                                                                <  1:3883913/93‑1 (MQ=255)
         cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCga                                                                                <  1:4155553/93‑1 (MQ=255)
          ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgag                                                                               <  2:683761/93‑1 (MQ=255)
            aCAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgaggg                                                                             <  2:6236439/93‑1 (MQ=255)
            aCAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgaggg                                                                             <  1:4892801/93‑1 (MQ=255)
              aaCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggag                                                                           <  1:3240467/93‑1 (MQ=255)
               aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggagg                                                                          >  2:6679828/1‑93 (MQ=255)
               aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggagg                                                                          >  2:9589577/1‑93 (MQ=255)
                ccGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGggagga                                                                         >  2:197340/1‑93 (MQ=255)
                 cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGTAGGaa                                                                        >  1:6829539/1‑93 (MQ=255)
                 cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGaa                                                                        <  1:858308/93‑1 (MQ=255)
                  gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAat                                                                       >  2:4269206/1‑93 (MQ=255)
                   agaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAata                                                                      >  1:7495429/1‑93 (MQ=255)
                     aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATc                                                                    >  1:8325893/1‑93 (MQ=255)
                       tttCCTGCGCTTTGGCAATACAATCACCCGCACTGACGGCAAGGCGTAACGTCGCCGTCGCGACAACAAGAATTTGCGAGGGATGAATATCtt                                                                  <  2:8527497/93‑1 (MQ=255)
                       tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCtt                                                                  >  2:8832396/1‑93 (MQ=255)
                       tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCtt                                                                  <  1:5221044/93‑1 (MQ=255)
                       tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCtt                                                                  >  2:9059505/1‑93 (MQ=255)
                         tCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATCTCTTGc                                                                >  1:2716782/1‑93 (MQ=255)
                          ccTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCa                                                               <  1:4526894/93‑1 (MQ=255)
                          ccTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCa                                                               <  2:873617/93‑1 (MQ=255)
                          ccTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCa                                                               <  2:7760889/93‑1 (MQ=255)
                            tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGa                                                             >  1:2466057/1‑93 (MQ=255)
                            tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGa                                                             <  2:8049897/93‑1 (MQ=255)
                             gcgcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGAc                                                            <  1:5853315/93‑1 (MQ=255)
                             gcgcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGAc                                                            <  2:9417267/93‑1 (MQ=255)
                             gcgcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGAc                                                            <  2:2453420/93‑1 (MQ=255)
                             gcgcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGAc                                                            >  2:1425204/1‑93 (MQ=255)
                               gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGt                                                          >  1:2023840/1‑93 (MQ=255)
                                cTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACTAATTTGCGAGGGAGGAATATCTTGCAGACGtt                                                         <  1:6296635/93‑1 (MQ=255)
                                cTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGtt                                                         <  1:2207367/93‑1 (MQ=255)
                                 tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                                                        >  1:10008069/1‑93 (MQ=255)
                                    ggCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGc                                                     <  1:7785752/93‑1 (MQ=255)
                                     gCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCa                                                    <  2:1707070/93‑1 (MQ=255)
                                     gCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCa                                                    >  2:157729/1‑93 (MQ=255)
                                        aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAAc                                                 <  1:6938872/93‑1 (MQ=255)
                                        aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAAc                                                 >  2:6997950/1‑93 (MQ=255)
                                         tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                                                >  2:8157434/1‑93 (MQ=255)
                                         tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGACGGAGGAATATCTTGCAGACGCTCAGCAAACa                                                >  2:537841/1‑93 (MQ=255)
                                          aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg                                               <  2:4448586/93‑1 (MQ=255)
                                           aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg                                              >  1:7592050/1‑93 (MQ=255)
                                           aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg                                              <  1:941953/93‑1 (MQ=255)
                                           aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg                                              <  2:3481948/93‑1 (MQ=255)
                                            aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                                             >  1:928965/1‑93 (MQ=255)
                                            aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                                             >  2:1868360/1‑93 (MQ=255)
                                             aTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcg                                            >  2:3218694/1‑93 (MQ=255)
                                                 gccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                                        >  1:2537475/2‑93 (MQ=255)
                                                 cccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                                        >  2:5917175/1‑93 (MQ=255)
                                                 cccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                                        >  2:1768966/1‑93 (MQ=255)
                                                  ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCACACAGGCGCAAAc                                       >  2:3726332/1‑93 (MQ=255)
                                                  ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAAc                                       <  2:7299165/93‑1 (MQ=255)
                                                  ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAAc                                       <  1:4241197/93‑1 (MQ=255)
                                                  ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAAc                                       >  1:6094937/1‑93 (MQ=255)
                                                   cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                                      >  2:121375/1‑93 (MQ=255)
                                                   cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                                      <  2:2034265/93‑1 (MQ=255)
                                                   cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                                      <  1:9801577/93‑1 (MQ=255)
                                                   cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                                      >  2:266207/1‑93 (MQ=255)
                                                   cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa                                      <  2:9177315/93‑1 (MQ=255)
                                                    gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                                     >  1:4836802/1‑93 (MQ=255)
                                                     cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAtt                                    <  2:7959510/93‑1 (MQ=255)
                                                     cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAtt                                    <  2:4794697/93‑1 (MQ=255)
                                                      aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTg                                   >  2:2353323/1‑93 (MQ=255)
                                                       ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGc                                  <  1:9098786/93‑1 (MQ=255)
                                                        tGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGcc                                 >  2:1481888/1‑93 (MQ=255)
                                                           cGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAc                              <  1:4460920/93‑1 (MQ=255)
                                                           cGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAc                              <  2:9783662/93‑1 (MQ=255)
                                                            ggCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAcc                             <  1:5935987/93‑1 (MQ=255)
                                                             gCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcg                            >  2:8406159/1‑93 (MQ=255)
                                                              cAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcgc                           >  2:55733/1‑93 (MQ=255)
                                                                aGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc                         >  2:3107956/1‑93 (MQ=255)
                                                                aGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc                         <  1:5642467/93‑1 (MQ=255)
                                                                aGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc                         >  1:1848096/1‑93 (MQ=255)
                                                                     tAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt                    <  1:1712819/93‑1 (MQ=255)
                                                                     tAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt                    >  1:6611138/1‑93 (MQ=255)
                                                                     tAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt                    >  1:3513263/1‑93 (MQ=255)
                                                                       acgtCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTg                  >  1:3631731/4‑93 (MQ=255)
                                                                         gtCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCt                >  1:8361719/2‑93 (MQ=255)
                                                                            gCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCa             >  1:1177117/1‑93 (MQ=255)
                                                                            gCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCa             <  1:899887/93‑1 (MQ=255)
                                                                             ccGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCAt            >  1:3114095/1‑93 (MQ=255)
                                                                              cGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCAtt           <  1:1606535/93‑1 (MQ=255)
                                                                              cGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCAtt           <  1:549808/93‑1 (MQ=255)
                                                                              cGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCAtt           <  1:9593965/93‑1 (MQ=255)
                                                                               gTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTg          >  2:4974743/1‑93 (MQ=255)
                                                                                tAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTgg         >  2:10004907/1‑93 (MQ=255)
                                                                                tAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCACAAACATTGCCAACCGCGCTCCATTGCTTCATTgg         >  2:2239626/1‑93 (MQ=255)
                                                                                 aGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGga        >  1:1468372/1‑93 (MQ=255)
                                                                                  gCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTCgag       >  2:6068061/1‑93 (MQ=255)
                                                                                   cGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGgaga      <  1:8147845/93‑1 (MQ=255)
                                                                                   cGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGgaga      <  1:8037129/93‑1 (MQ=255)
                                                                                   cGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGgaga      <  2:5160588/93‑1 (MQ=255)
                                                                                   cGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGgaga      >  2:5753504/1‑93 (MQ=255)
                                                                                     acaacaCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAgg    >  2:7547643/1‑93 (MQ=255)
                                                                                     acaacaCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGACAACCGCGCTCCATTGCTTCATTGGAGAgg    >  1:3028786/1‑93 (MQ=255)
                                                                                       aaCACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTCCCAACCGCGCTCCATTGCTTCATTGGAGAggga  >  2:2874354/1‑92 (MQ=255)
                                                                                          |                                                                                         
CACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATTGCTTCATTGGAGAGGGC  >  NC_000913/3961515‑3961694

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