Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 569,161 A→C 15.4% D87A (GAC→GCC) ‡ ybcK → DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913569,1610AC15.4% 359.6 / 61.1 149D87A (GAC→GCC) ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base A (70/55);  new base C (10/13);  total (80/69)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

Only 100 of 152 total aligned reads displayed.

AAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCT  >  NC_000913/569062‑569259
                                                                                                   |                                                                                                  
aaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAGTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGa                                                                                                    >  2:3564939/1‑100 (MQ=255)
aaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGa                                                                                                    >  2:3069991/1‑100 (MQ=255)
  ggAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACa                                                                                                  >  1:1466537/1‑100 (MQ=255)
  ggAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACa                                                                                                  >  2:3339307/1‑100 (MQ=255)
    aaaGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGa                                                                                                >  2:2977586/1‑100 (MQ=255)
     aaGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGAc                                                                                               >  2:3063645/1‑100 (MQ=255)
        cATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACttt                                                                                            >  2:1748201/1‑100 (MQ=255)
         aTGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTc                                                                                           <  2:991441/100‑1 (MQ=255)
          tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                          <  1:2584344/100‑1 (MQ=255)
          tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                          >  2:223574/1‑100 (MQ=255)
          tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                          >  1:3716209/1‑100 (MQ=255)
            cataaTGAGGAGCTATTTTCGAATTTGTAGATGCTATAGAGCATGGTTATATTATGACAGCAACTACATAGTTAGTTGAAAGTCTGTCCAGCCTTTCacc                                                                                        <  1:2704220/97‑3 (MQ=255)
            cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  2:2022490/100‑2 (MQ=255)
            cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  2:3673864/100‑2 (MQ=255)
            cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                        >  2:1209666/1‑100 (MQ=255)
            cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                        >  1:320525/1‑100 (MQ=255)
             acaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                        <  1:745082/99‑2 (MQ=255)
             acaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTAATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAaga                                                                                       >  1:2118845/1‑100 (MQ=255)
              caATCAGGAGCTTTTTCGGAATTTTTAGATGCTATATAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaacta                                                                                      <  2:2673154/100‑4 (MQ=38)
              caATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAagag                                                                                      >  2:46143/1‑100 (MQ=255)
               aaTCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaactaa                                                                                     <  1:1441529/100‑5 (MQ=39)
               aaTCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaacta                                                                                      <  1:3621441/99‑4 (MQ=39)
                    ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCCGGAACTACATTGTTCGTTGAAAGTCTGGACAGACTTTCAAGTGAAAAAg                                                                                >  2:3643395/1‑100 (MQ=255)
                    ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAg                                                                                >  1:3627391/1‑100 (MQ=255)
                    ggAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAg                                                                                >  2:2730343/1‑100 (MQ=255)
                      aGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTc                                                                              >  1:3147750/1‑100 (MQ=255)
                      aGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAAATACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGt                                                                               >  1:3147597/1‑99 (MQ=255)
                        cTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCgg                                                                            <  1:419132/100‑1 (MQ=255)
                        cTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCgg                                                                            >  2:113562/1‑100 (MQ=255)
                             tCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAg                                                                       <  2:2788889/100‑1 (MQ=255)
                              cGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGc                                                                      <  2:1606380/100‑1 (MQ=255)
                              cGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGc                                                                      <  2:346244/100‑1 (MQ=255)
                               ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCg                                                                     <  2:1574667/100‑1 (MQ=255)
                                   tttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTg                                                                 <  2:435804/100‑1 (MQ=255)
                                   tttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTg                                                                 <  1:2766400/100‑1 (MQ=255)
                                      ttAGATGCTATAGAGCATGGTTATATAGTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGTGATTGAAc                                                              <  2:1447647/100‑1 (MQ=255)
                                       tAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACg                                                             <  1:1745564/100‑1 (MQ=255)
                                        atattcTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGt                                                            <  1:2265592/95‑1 (MQ=255)
                                           tGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTcgg                                                         <  2:2586379/100‑3 (MQ=255)
                                           tGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTg                                                         >  1:3398262/1‑100 (MQ=255)
                                             cTATAGAGCAAGATTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACATTCTAGAAAAAAAGACGGTGAAGCGATTGAATGt                                                            >  1:1535843/1‑95 (MQ=255)
                                              tataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGaaa                                                      >  1:1820465/1‑100 (MQ=255)
                                              tataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGaaa                                                      >  2:1822560/1‑100 (MQ=255)
                                               ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAt                                                     <  1:1451386/100‑1 (MQ=255)
                                               ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAt                                                     >  2:869860/1‑100 (MQ=255)
                                                  gagCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGa                                                  >  1:668652/1‑100 (MQ=255)
                                                   agCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                 <  2:2469372/100‑1 (MQ=255)
                                                   agCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                 <  2:958992/100‑1 (MQ=255)
                                                    gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGaag                                                >  1:1349512/1‑98 (MQ=255)
                                                      aTGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtttt                                              >  2:1855627/1‑100 (MQ=255)
                                                      aTGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtttt                                              <  2:1065488/100‑1 (MQ=255)
                                                      aTGGTTATATATTGAAATGAACTACATTATTAGTTGAAAGTGTGTACAGACTTTTAAGAGAAATAGTTGGTTTAGCGATTGAACGTCAGAAATTGAtt                                                <  1:19224/98‑1 (MQ=255)
                                                       tGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCCGGACAGACTTTCAAGAAAAGAAGTAGGTGAATCGATTGAACGTCTGAAATTGATTTTg                                             >  1:2072254/1‑100 (MQ=255)
                                                          ttATATATTGCCAGGAACTACATTGTTAGTTGGAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAAt                                          >  2:3463816/1‑100 (MQ=255)
                                                           tatataTTGCGAGGAACTACATTGTTAGTTGATAGTAGGGGCAGACTGTCAAGAGATCGATTCGGTCCAGCTATTGGACGGCTGCAATTGATTTTGTATc                                         <  1:1719274/100‑1 (MQ=255)
                                                           tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATc                                         <  2:1373727/100‑1 (MQ=255)
                                                             tatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCAc                                       >  1:3462144/1‑100 (MQ=255)
                                                             tatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCGGAAATTGATTTTGAATCAc                                       <  1:1418828/100‑1 (MQ=255)
                                                              atatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACg                                      <  2:202319/100‑1 (MQ=255)
                                                              atatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACg                                      <  2:2118855/100‑1 (MQ=255)
                                                                   gCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAtt                                 >  1:926684/1‑100 (MQ=255)
                                                                    ccAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTg                                >  2:1747009/1‑100 (MQ=255)
                                                                      aGGAACTACATTGTTAGTTGAAAGTCTTGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGAt                              <  2:1264261/100‑1 (MQ=255)
                                                                       ggAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATg                             >  2:320220/1‑100 (MQ=255)
                                                                         aaCTAGATTGTTAGTTGAAAGTCTGGATAGACAAACAAAAAAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATg                             <  1:186002/98‑1 (MQ=255)
                                                                          aCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTta                          <  2:683996/100‑1 (MQ=255)
                                                                          aCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTta                          <  1:1032312/100‑1 (MQ=255)
                                                                           cTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtat                         >  2:1858789/1‑100 (MQ=255)
                                                                           cTACATTGGTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtat                         >  2:3497810/1‑100 (MQ=255)
                                                                            tACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtata                        <  1:1684417/100‑1 (MQ=255)
                                                                            tACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtata                        <  1:932/100‑1 (MQ=255)
                                                                             aCATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATaa                       >  2:3346835/1‑100 (MQ=255)
                                                                              cATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAAc                      <  1:1568899/100‑1 (MQ=255)
                                                                               tttttatcttGTAAGTCTCGACAGACTTTCAAGAGATAAAGTCAGCGAAGCGATTGAACGTATGAAATTGATTTTGAATCACGGTATTGATGTTATAAc                      <  1:2191185/91‑1 (MQ=255)
                                                                               aTTGTTAGTTTAAAGTCTGGACAGACTTTCAAGAGAAAAATTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTCAATCACGGTATTGATGTTATAAc                      <  2:3232985/99‑1 (MQ=255)
                                                                                ttGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActc                    >  1:2304681/1‑100 (MQ=255)
                                                                                 tGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActct                   >  2:2564621/1‑100 (MQ=255)
                                                                                     aGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc               >  2:712154/1‑100 (MQ=255)
                                                                                     aGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc               <  2:1493583/100‑1 (MQ=255)
                                                                                      gTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg              >  2:1178109/1‑100 (MQ=255)
                                                                                      cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg              >  1:280699/2‑100 (MQ=255)
                                                                                      cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg              >  2:2519359/2‑100 (MQ=255)
                                                                                       ttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGa             >  2:108982/1‑100 (MQ=255)
                                                                                       ttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGa             >  2:826786/1‑100 (MQ=255)
                                                                                        tGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc            >  1:3275172/1‑100 (MQ=255)
                                                                                        tGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc            >  2:2146791/1‑100 (MQ=255)
                                                                                        tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc            <  1:2745796/100‑1 (MQ=255)
                                                                                         gAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa           >  1:2192971/1‑100 (MQ=255)
                                                                                          aaaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGGGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACaa          <  1:6699/100‑1 (MQ=255)
                                                                                            aGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa        <  2:2177809/100‑1 (MQ=255)
                                                                                            aGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa        <  2:1355341/100‑1 (MQ=255)
                                                                                             gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATAc       >  2:959017/1‑100 (MQ=255)
                                                                                             gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATATTGAATCAAGGTATTGATGTTATAAATCTTTGCGACAATAc       >  1:3800564/1‑100 (MQ=255)
                                                                                                tGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGt    >  2:116560/1‑100 (MQ=255)
                                                                                                 ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc   >  2:3007563/1‑100 (MQ=255)
                                                                                                 ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc   <  1:552564/100‑1 (MQ=255)
                                                                                                 ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc   >  2:1774148/1‑100 (MQ=255)
                                                                                                  gACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTGGCGACAATACAGTCt  >  2:3761531/1‑100 (MQ=255)
                                                                                                  gACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGAATTGATGTTATAACTCTTTGCGAAAATACAGTCt  >  2:3748806/1‑100 (MQ=255)
                                                                                                   aCAGACTTTCAAGAGAAAAAGTCGGTGAAGCTATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCt  <  1:759726/99‑1 (MQ=255)
                                                                                                   |                                                                                                  
AAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCT  >  NC_000913/569062‑569259

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