Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 309,281 (G)6→7 intergenic (‑31/+27) yagY ← / ← yagZ hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913309,2751.G92.9% 682.3 / 43.6 225intergenic (‑25/+33)yagY/yagZhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (5/11);  new base G (126/82);  total (132/93)
Fisher's exact test for biased strand distribution p-value = 3.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 266 total aligned reads displayed.

ACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCT  >  NC_000913/309183‑309362
                                                                                             |                                                                                         
aCATCCGGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCa                                                                                            <  2:2728745/93‑1 (MQ=255)
 cATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAg                                                                                          <  2:5301763/93‑1 (MQ=255)
 cATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAg                                                                                          >  2:2639193/1‑93 (MQ=255)
  aTCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgg                                                                                        <  2:4305664/93‑1 (MQ=255)
   tCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAggg                                                                                       <  2:3006583/93‑1 (MQ=255)
   tCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAggg                                                                                       <  1:5217758/93‑1 (MQ=255)
    ccAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgggg                                                                                      <  2:663990/93‑1 (MQ=255)
    ccAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgggg                                                                                      >  1:5745394/1‑93 (MQ=255)
     cAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAggggg                                                                                    <  1:118067/93‑1 (MQ=255)
     cAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAggggg                                                                                    <  1:6059173/93‑1 (MQ=255)
       gcgcCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAggggggg                                                                                  >  2:3804614/1‑92 (MQ=255)
       gcgcCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGc                                                                                  >  1:961682/1‑93 (MQ=255)
        cgcCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgggggggc                                                                                 <  1:4971991/93‑3 (MQ=255)
         gcCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgggggggct                                                                                >  1:2378582/1‑90 (MQ=255)
         gcCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgggggggct                                                                                >  2:93953/1‑90 (MQ=255)
             gggCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTAccc                                                                             <  1:5461322/93‑1 (MQ=255)
             gggCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTAccc                                                                             >  1:3879573/1‑93 (MQ=255)
             gggCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTAccc                                                                             >  1:1505346/1‑93 (MQ=255)
             gggCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTAcccc                                                                            <  2:350771/93‑1 (MQ=255)
              ggCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTAcccc                                                                            <  2:2365604/93‑1 (MQ=255)
              ggCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTAcccc                                                                            <  1:3735049/93‑1 (MQ=255)
              ggCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTAcccc                                                                            <  1:4431267/93‑1 (MQ=255)
                 cGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCctgc                                                                         >  1:4889481/1‑93 (MQ=255)
                  ggAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCctgct                                                                        >  1:1455595/1‑93 (MQ=255)
                  ggAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCctgct                                                                        >  2:5517674/1‑93 (MQ=255)
                     gaTATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGt                                                                     >  2:4526998/1‑93 (MQ=255)
                       tatTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTAc                                                                   <  2:867499/93‑1 (MQ=255)
                       tatTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTAc                                                                   >  2:759073/1‑93 (MQ=255)
                        atTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACa                                                                  >  1:2774944/1‑93 (MQ=255)
                        atTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACa                                                                  <  2:6338972/93‑1 (MQ=255)
                        atTCCGGAAAACAGCAGAGAGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACa                                                                  <  1:5584474/93‑1 (MQ=255)
                          tCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATc                                                                <  1:1311977/93‑1 (MQ=255)
                            cGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCaga                                                             <  1:4531719/93‑1 (MQ=255)
                             ggAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCaga                                                             <  1:722012/93‑1 (MQ=255)
                                aaaCAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCagagag                                                          <  2:361255/93‑1 (MQ=255)
                                aaaCAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCagagag                                                          >  1:6665884/1‑93 (MQ=255)
                                aaaCAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAgagaga                                                         <  2:341203/93‑1 (MQ=255)
                                  aCAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGAt                                                        >  1:189419/1‑93 (MQ=255)
                                  aCAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGAt                                                        <  2:4735630/93‑1 (MQ=255)
                                  aCAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGAt                                                        <  1:778710/93‑1 (MQ=255)
                                  aCAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGAt                                                        >  2:811509/1‑93 (MQ=255)
                                    agcagAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTa                                                      >  1:515563/1‑93 (MQ=255)
                                    agcagAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTa                                                      >  2:5006048/1‑93 (MQ=255)
                                      cagAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAAc                                                    <  1:3509787/93‑1 (MQ=255)
                                       agagCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACt                                                   <  1:1255888/93‑1 (MQ=255)
                                       agagCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACt                                                   >  2:4866454/1‑93 (MQ=255)
                                        gagCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTg                                                  >  2:2984261/1‑93 (MQ=255)
                                        gagCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTg                                                  >  1:1675192/1‑93 (MQ=255)
                                        gagCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTg                                                  >  2:5412352/1‑93 (MQ=255)
                                         agCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTgg                                                 <  1:3392797/93‑1 (MQ=255)
                                           cGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTACCTGGTc                                               <  1:3315502/93‑1 (MQ=255)
                                            gagaGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTcc                                              <  2:3922401/93‑1 (MQ=255)
                                            gagaGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTcc                                              >  2:1679556/1‑93 (MQ=255)
                                             agagGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCa                                             >  2:6184627/1‑93 (MQ=255)
                                             agagGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCa                                             >  1:2668722/1‑93 (MQ=255)
                                               agGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAgg                                           >  1:4092146/1‑93 (MQ=255)
                                                ggCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGt                                          <  2:511928/93‑1 (MQ=255)
                                                ggCAGAAGATGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGt                                          >  1:6547199/1‑93 (MQ=255)
                                                   aGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTcgc                                       <  1:4120032/93‑1 (MQ=255)
                                                   aGAAGGCGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTcgc                                       >  1:538121/1‑93 (MQ=255)
                                                    gAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTcgcg                                      >  2:1971744/1‑93 (MQ=255)
                                                    gAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTcgcg                                      >  1:2306764/1‑93 (MQ=255)
                                                    gAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTcgcg                                      >  1:3243671/1‑93 (MQ=255)
                                                     aaGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGt                                     >  1:3199632/1‑93 (MQ=255)
                                                     aaGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGt                                     >  1:4789223/1‑93 (MQ=255)
                                                     aaGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTc                                    >  2:2284563/1‑93 (MQ=255)
                                                      aGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTc                                    >  1:6190654/1‑93 (MQ=255)
                                                      aGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTc                                    >  1:2262069/1‑93 (MQ=255)
                                                      aGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCg                                   >  1:3695762/1‑93 (MQ=255)
                                                         tGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGaa                                 >  2:2869386/1‑93 (MQ=255)
                                                          gCTTTTTCATAATATTTCCGGCCAGGAATAAAGCAGGGGATTCTCCCCCTCCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAAc                                <  1:1984787/93‑1 (MQ=255)
                                                             ttttCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGt                             >  2:1347651/1‑93 (MQ=255)
                                                             ttttCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGt                             >  2:4666866/1‑93 (MQ=255)
                                                             ttttCATAACAATTCCGTCCAGGAATAAAGCAGTGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAGCTGGTCCAGGTCGCGTCGAACTGt                             <  1:504598/93‑1 (MQ=255)
                                                               ttCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTAc                           >  1:1550221/1‑93 (MQ=255)
                                                                tCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACg                          >  2:3216247/1‑93 (MQ=255)
                                                                tCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGc                         <  2:3838995/93‑1 (MQ=255)
                                                                  aTAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCt                        >  1:5429601/1‑93 (MQ=255)
                                                                  aTAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCt                        >  2:5144648/1‑93 (MQ=255)
                                                                   tAACAATTGCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGGCGCGTCGAACTGTACGCTa                       >  1:6529274/1‑93 (MQ=255)
                                                                   tAACAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTa                       <  1:2944127/93‑1 (MQ=255)
                                                                     aCAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAAc                     <  1:6352550/93‑1 (MQ=255)
                                                                     aCAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAAc                     >  2:6051035/1‑93 (MQ=255)
                                                                      cAATTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTAAATCAGAGAGATTACTTGGCCAAGGGCGGGTCGATCTTTACGCTACCg                    >  2:1446200/1‑93 (MQ=255)
                                                                       aaTTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGt                   >  2:4659611/1‑93 (MQ=255)
                                                                       aaTTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGt                   >  2:1066860/1‑93 (MQ=255)
                                                                       aaTTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGt                   <  1:1038730/93‑1 (MQ=255)
                                                                        aTTCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTc                  >  2:5741236/1‑93 (MQ=255)
                                                                         ttCCGTCCAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTcg                 >  2:1562019/1‑93 (MQ=255)
                                                                                cAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTcc          >  2:1646598/1‑93 (MQ=255)
                                                                                cAGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTcc          >  1:5569204/1‑93 (MQ=255)
                                                                                 aGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCa         >  2:5655415/1‑93 (MQ=255)
                                                                                 aGGAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCa         >  2:6216481/1‑93 (MQ=255)
                                                                                  ggAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAg        >  2:6575850/1‑93 (MQ=255)
                                                                                  ggAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGACGCTCCAg        >  1:927195/1‑93 (MQ=255)
                                                                                   gAATAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGa       <  2:4598468/93‑1 (MQ=255)
                                                                                    aaTAAAGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGAt      >  2:249801/1‑93 (MQ=255)
                                                                                      tAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGcc   <  2:4813149/93‑1 (MQ=255)
                                                                                       aaaGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGcc   <  2:589355/93‑1 (MQ=255)
                                                                                       aaaGCAGGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGcc   <  1:1848571/93‑1 (MQ=255)
                                                                                             |                                                                                         
ACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCT  >  NC_000913/309183‑309362

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: