Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,017,557 A→G E116G (GAG→GGG)  yigP → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,017,5570AG99.7% 1183.1 / ‑6.2 380E116G (GAG→GGG) yigPhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (196/183);  minor base C (1/0);  total (197/183)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 390 total aligned reads displayed.

AGCTTACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGAGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGTTCCTGC  >  NC_000913/4017467‑4017649
                                                                                          |                                                                                             
agcTTACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGt                                                                                             <  2:1979727/93‑1 (MQ=255)
agcTTACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGt                                                                                             <  2:286539/93‑1 (MQ=255)
 gcTTACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGtt                                                                                            >  2:4231880/1‑93 (MQ=255)
  cTTACCGCACTGATTCGCAGGGGTGAGCTGGAAGTGCAGTGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTc                                                                                           >  1:2022181/1‑93 (MQ=255)
    tACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGa                                                                                         <  1:3489931/93‑1 (MQ=255)
    tACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGa                                                                                         <  1:1896394/93‑1 (MQ=255)
      ccGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGAcc                                                                                       <  2:3362528/93‑1 (MQ=255)
       cGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGAccc                                                                                      >  2:4029382/1‑93 (MQ=255)
        gCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCt                                                                                     >  2:867875/1‑93 (MQ=255)
          aCTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGc                                                                                   <  2:6953283/93‑1 (MQ=255)
          aCTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGc                                                                                   <  2:7628554/93‑1 (MQ=255)
             gATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGa                                                                                >  1:624335/1‑93 (MQ=255)
              aTTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCGGGGTTCGACCCTGCGGaa                                                                               >  1:5322558/1‑93 (MQ=255)
              aTTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGaa                                                                               >  1:1334123/1‑93 (MQ=255)
              aTTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGaa                                                                               <  1:2984442/93‑1 (MQ=255)
               ttCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGCAc                                                                              >  2:6653970/1‑93 (MQ=255)
               ttCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAAc                                                                              <  2:6541695/93‑1 (MQ=255)
               ttCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAAc                                                                              >  1:1086791/1‑93 (MQ=255)
                   cAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAActgct                                                                          <  2:133579/93‑1 (MQ=255)
                   cAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAActgct                                                                          <  1:5801862/93‑1 (MQ=255)
                     gtggtgAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTgg                                                                        <  2:1048241/93‑1 (MQ=255)
                     gtggtgAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTgg                                                                        >  1:7292921/1‑93 (MQ=255)
                      tggtgAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGc                                                                       >  1:1213055/1‑93 (MQ=255)
                       ggtgAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGcc                                                                      >  1:4082509/1‑93 (MQ=255)
                           aGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCtt                                                                 <  1:465092/93‑1 (MQ=255)
                            gCTGGAGGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTta                                                                <  1:321098/93‑1 (MQ=255)
                            gCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTta                                                                >  1:7249323/1‑93 (MQ=255)
                            gCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTta                                                                <  1:1637271/93‑1 (MQ=255)
                             cTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTtat                                                               <  2:7482596/93‑1 (MQ=255)
                              tGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGGACTGCTGGCCCCTtata                                                              <  1:5700902/93‑1 (MQ=255)
                               ggAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATAc                                                             <  2:3774433/93‑1 (MQ=255)
                               ggAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATAc                                                             <  2:4718743/93‑1 (MQ=255)
                               ggAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATAc                                                             <  1:7882166/93‑1 (MQ=255)
                               ggAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATAc                                                             <  1:3576801/93‑1 (MQ=255)
                                gAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGTCCCCTTATAcc                                                            >  1:4582634/1‑93 (MQ=255)
                                  aGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCgg                                                          <  2:7302263/93‑1 (MQ=255)
                                   gTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCggg                                                         >  1:5645902/1‑92 (MQ=255)
                                   gTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGt                                                         >  1:6666390/1‑93 (MQ=255)
                                    tGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTg                                                        >  2:3319716/1‑93 (MQ=255)
                                      cAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGat                                                      >  2:260823/1‑93 (MQ=255)
                                      cAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGat                                                      <  2:5665079/93‑1 (MQ=255)
                                       aGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGata                                                     >  1:5447797/1‑93 (MQ=255)
                                       aGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGata                                                     <  2:593137/93‑1 (MQ=255)
                                       aGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGata                                                     <  1:3799986/93‑1 (MQ=255)
                                       aGGGCGATATTCAGGTGGTACAATACTTCGTTGCGCTCGGAGATCTGGCAGGGGTCGACCCTGCGGAACTGCTTGCCCCTTATACCGGTGata                                                     >  2:2917099/1‑93 (MQ=255)
                                        gggCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGatat                                                    <  1:744972/93‑1 (MQ=255)
                                        gggCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGatat                                                    >  2:5996588/1‑93 (MQ=255)
                                        gggCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGatat                                                    >  1:5400968/1‑93 (MQ=255)
                                        gggCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGatat                                                    <  1:5015430/93‑1 (MQ=255)
                                            gATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATcgcc                                                >  2:2319938/1‑93 (MQ=255)
                                            gATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATcgcc                                                >  1:4032684/1‑93 (MQ=255)
                                            gATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATcgcc                                                >  1:3936084/1‑93 (MQ=255)
                                              tatTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATcgccgc                                              <  1:3693855/93‑1 (MQ=255)
                                               atTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCt                                             >  1:2439715/1‑93 (MQ=255)
                                               atTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCt                                             <  1:2164330/93‑1 (MQ=255)
                                                ttCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTTATACCGGTGATATcgccgcg                                             >  2:7846654/1‑92 (MQ=255)
                                                 tCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGa                                           >  2:5358450/1‑93 (MQ=255)
                                                 tCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGa                                           <  2:1649984/93‑1 (MQ=255)
                                                 tCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGa                                           <  1:5835613/93‑1 (MQ=255)
                                                   aGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAg                                         >  2:1217128/1‑93 (MQ=255)
                                                   aGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAg                                         <  1:810882/93‑1 (MQ=255)
                                                       ggtgCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAAt                                     <  1:7627674/93‑1 (MQ=255)
                                                       ggtgCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAAt                                     >  2:4375162/1‑93 (MQ=255)
                                                         tgCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCa                                   <  1:1492511/93‑1 (MQ=255)
                                                         tgCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCa                                   >  1:2304731/1‑93 (MQ=255)
                                                         tgCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCa                                   <  1:6886964/93‑1 (MQ=255)
                                                          gCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAg                                  >  2:5831925/1‑93 (MQ=255)
                                                           cAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGc                                 >  1:6192545/1‑93 (MQ=255)
                                                            aaaaCTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCa                                >  1:4803793/1‑93 (MQ=255)
                                                            aaaaCTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCa                                <  2:5419892/93‑1 (MQ=255)
                                                               aCTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAg                             <  1:1012304/93‑1 (MQ=255)
                                                               aCTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAg                             >  1:4932586/1‑93 (MQ=255)
                                                                cTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGc                            >  2:1487189/1‑93 (MQ=255)
                                                                cTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGc                            >  2:1821731/1‑93 (MQ=255)
                                                                cTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGc                            <  1:5171804/93‑1 (MQ=255)
                                                                cTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGc                            >  2:503701/1‑93 (MQ=255)
                                                                cTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGc                            >  1:6408272/1‑93 (MQ=255)
                                                                   cGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCAt                         <  2:2206033/93‑1 (MQ=255)
                                                                    gTTGCGCTGGCAGATCTGGTAGGGTTCGACCCCGCGGAACTGCTGGCTCCTTATACCGGTGATATCGCTGCTGTATGAATCAGCAAAGCCATg                        >  1:7033936/1‑93 (MQ=255)
                                                                    gTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATg                        <  1:5971647/93‑1 (MQ=255)
                                                                       gcgcTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATgcgc                     <  1:3247180/93‑1 (MQ=255)
                                                                         gcTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCCAAGCCATGCGCgg                   >  2:5518588/1‑93 (MQ=255)
                                                                         gcTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCgg                   >  1:6822277/1‑93 (MQ=255)
                                                                         gcTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCgg                   >  1:7752019/1‑93 (MQ=255)
                                                                          cTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGa                  <  2:6716428/93‑1 (MQ=255)
                                                                            ggCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCCGCAAAGCCATGCGCGGAgg                <  1:3103932/93‑1 (MQ=255)
                                                                            ggCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAgg                >  2:754105/1‑93 (MQ=255)
                                                                            ggCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAgg                <  2:3609721/93‑1 (MQ=255)
                                                                             gCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGgc               <  2:4936693/93‑1 (MQ=255)
                                                                                 aTCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCaa           >  1:6503915/1‑93 (MQ=255)
                                                                                 aTCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCaa           >  2:4634168/1‑93 (MQ=255)
                                                                                 aTCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCaa           >  1:451412/1‑93 (MQ=255)
                                                                                   cTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAATGAATCAGCAAAGCCATGCGCGGAGGCGCAAAg         <  2:7725045/93‑1 (MQ=255)
                                                                                    tGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGt        >  2:3179925/1‑93 (MQ=255)
                                                                                    tGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGt        <  1:4238133/93‑1 (MQ=255)
                                                                                     ggCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGtt       <  2:2453041/93‑1 (MQ=255)
                                                                                     ggCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGtt       <  2:6237929/93‑1 (MQ=255)
                                                                                      gCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGGGGCGCAAAGTTc      >  1:2681759/1‑93 (MQ=255)
                                                                                       cAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGTTcc     <  1:286940/93‑1 (MQ=255)
                                                                                          agTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGTTCCTGc  <  2:6458866/93‑1 (MQ=255)
                                                                                          |                                                                                             
AGCTTACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGAGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGTTCCTGC  >  NC_000913/4017467‑4017649

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