Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,503,627 A→G 31.5% intergenic (‑196/‑431) insG ← / → yjhB IS4 transposase/putative MFS transporter, membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,503,6270AG31.5% 156.2 / 101.4 133intergenic (‑196/‑431)insG/yjhBIS4 transposase/putative MFS transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base A (43/47);  new base G (21/20);  total (66/67)
Fisher's exact test for biased strand distribution p-value = 8.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

GTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTG  >  NC_000913/4503528‑4503726
                                                                                                   |                                                                                                   
gtgtTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGta                                                                                                     <  2:1465584/100‑1 (MQ=255)
 tgtTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGtat                                                                                                    <  2:2932195/100‑1 (MQ=255)
  gtTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGtata                                                                                                   <  2:3043872/100‑1 (MQ=255)
      gTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAg                                                                                               <  2:666154/100‑1 (MQ=255)
      gTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAg                                                                                               <  1:2547071/100‑1 (MQ=255)
      gTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAg                                                                                               >  1:1680925/1‑100 (MQ=255)
       tAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAgg                                                                                              >  2:2353116/1‑100 (MQ=255)
       tAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAgg                                                                                              <  2:1266362/100‑1 (MQ=255)
       tAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAgg                                                                                              <  1:491634/100‑1 (MQ=255)
       tAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAgg                                                                                              <  1:638163/100‑1 (MQ=255)
       tAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGAATTGCTTGTTTTTTGATCTGGTATAACAgg                                                                                              >  1:1959963/1‑100 (MQ=255)
       tAACAAACTATCAGATACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACa                                                                                                >  2:917227/1‑98 (MQ=255)
        aaCAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGt                                                                                             <  2:954175/100‑1 (MQ=255)
        aaCAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGt                                                                                             >  2:1532469/1‑100 (MQ=255)
        aaCAAACTATCAGATAAAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGt                                                                                             <  2:2507609/100‑1 (MQ=255)
         aCAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGta                                                                                            <  1:492086/100‑1 (MQ=255)
         aCAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGta                                                                                            <  2:2153784/100‑1 (MQ=255)
         aCAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGta                                                                                            <  1:1295481/100‑1 (MQ=255)
               tATCAGCTTCAAAAATATGCTCAATTTTTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAgg                                                                                      <  1:1330050/100‑1 (MQ=255)
               tATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGTGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAgg                                                                                      <  2:2093493/100‑1 (MQ=255)
               tATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAgg                                                                                      <  1:938886/100‑1 (MQ=255)
                aTCAGCTACAAAAATATGCACAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTAGTATAACAGGTATAAAg                                                                                       >  2:283443/1‑98 (MQ=255)
                  cAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGTTGTAACAGTTATAAAGGtat                                                                                   >  2:467630/1‑100 (MQ=255)
                   aGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGtata                                                                                  >  1:2344689/1‑100 (MQ=255)
                    gCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATAc                                                                                 <  2:1374028/100‑1 (MQ=255)
                    gCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATAc                                                                                 >  1:2719241/1‑100 (MQ=255)
                     cTACAAAAATATGCTCAATTTGTTACATCAGTAACAAAACGCGTTTTGTTATGTGGATTTCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATAc                                                                                 <  2:97253/99‑1 (MQ=255)
                      tACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAg                                                                               >  1:1728251/1‑100 (MQ=255)
                      tACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAg                                                                               >  1:1341631/1‑100 (MQ=255)
                           aaaTTTGTTTAATTTGTGAGATCTGTAACAATACGCGTTTTGTTATGTGTATTGCTTGTTTTTTGATCTTTTATATCAGGTATAAAGGTATACAGAAAAg                                                                          <  2:1484027/100‑1 (MQ=255)
                           aaaTATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAg                                                                          <  2:1728231/100‑1 (MQ=255)
                                tGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAg                                                                     <  2:552335/100‑1 (MQ=255)
                                 gCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGa                                                                    >  2:1657931/1‑100 (MQ=255)
                                  cTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGaa                                                                   >  1:2633128/1‑100 (MQ=255)
                                  cTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGaa                                                                   >  1:1635103/1‑100 (MQ=255)
                                      aTTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATAc                                                               >  2:2533635/1‑100 (MQ=255)
                                      aTTTGTGAAATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATAc                                                               <  1:1259630/100‑1 (MQ=255)
                                        ttGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTg                                                             <  1:2852739/100‑1 (MQ=255)
                                           tgACATCAGTAAGAAAACGCGTTTTGTTATGTAGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGc                                                            <  2:2700285/98‑1 (MQ=255)
                                            aatccacatAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCaaa                                                         >  1:553315/9‑100 (MQ=38)
                                              cATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAgg                                                       <  2:2752479/100‑1 (MQ=255)
                                               aTCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGa                                                      <  2:696318/100‑1 (MQ=255)
                                                 cAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGaaa                                                    >  1:106473/1‑100 (MQ=255)
                                                   gTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAAc                                                  >  1:2578738/1‑100 (MQ=255)
                                                    tAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACa                                                 <  1:1662519/100‑1 (MQ=255)
                                                    tAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAAAAAGTATAAAAGTATACAGAAAAGCAAAGAAATAATGCAAAGGAAAACa                                                 >  1:421820/1‑100 (MQ=38)
                                                     aaCAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACa                                                 <  2:925347/99‑1 (MQ=255)
                                                     aaCAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAg                                                >  2:2430407/1‑100 (MQ=255)
                                                       cAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCt                                              >  1:2081139/1‑100 (MQ=255)
                                                        aaaaCGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCta                                             >  2:258125/1‑100 (MQ=255)
                                                          aaCGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCtata                                           <  1:1543862/100‑1 (MQ=255)
                                                            cgcgTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAAc                                         >  2:1542702/1‑100 (MQ=255)
                                                            cgcgTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAAc                                         <  2:1341619/100‑1 (MQ=255)
                                                               gTTTTGTTATGTGGATTGCTTTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAGGAAATACTGCAAAGGAAAACAGCTATAACg                                        <  2:1918036/98‑1 (MQ=255)
                                                               gTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTa                                      <  2:1298008/100‑1 (MQ=255)
                                                                 tttGTTATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAg                                    <  2:295516/100‑1 (MQ=255)
                                                                  ttGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGc                                   <  1:2918788/100‑1 (MQ=255)
                                                                     ttATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTaa                                >  1:1060302/1‑100 (MQ=255)
                                                                      tATGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTaaa                               <  1:1198880/100‑1 (MQ=255)
                                                                       aTGTGGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAg                              <  2:2362147/100‑1 (MQ=255)
                                                                        tgtgGATTGCTTGTTTTTTTATTTTTTATAAAAGGGAAAAAGGTATACAGAAAAGCAAAGAAATAGTGCAAACGAAAACAGCTATAACCTAAGCTAAAGt                             >  2:2400283/1‑100 (MQ=255)
                                                                        tgtgGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGt                             >  2:1284224/1‑100 (MQ=255)
                                                                        tgtgGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGt                             >  2:1810441/1‑100 (MQ=255)
                                                                        tgtgGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGt                             >  1:1973262/1‑100 (MQ=255)
                                                                        tgtgGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGt                             <  1:983655/100‑1 (MQ=255)
                                                                         gtgGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTACAAAGGAAAACAGCTATAACGTAAGCTAAAGta                            >  1:353360/1‑100 (MQ=255)
                                                                         gtgGATTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAACACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGta                            <  2:1607471/100‑1 (MQ=255)
                                                                           ggATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGtaat                          <  2:1476266/100‑1 (MQ=255)
                                                                            gATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGtaata                         <  1:2406272/100‑1 (MQ=255)
                                                                             aTTGCTTGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGtaataa                        >  1:2459780/1‑100 (MQ=255)
                                                                                  ttGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACctct                   >  1:2645074/1‑100 (MQ=255)
                                                                                  ttGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACctct                   >  1:2732150/1‑100 (MQ=255)
                                                                                   tGTTTTTTTTTTTTGTTTAACAGGGATAAAGGTATACAGAAAAGCAAAGAAATACTGGAAAGGAAAACAGCTATAACGTAAGCTAAAGTTATAACctct                   >  2:1507627/1‑99 (MQ=255)
                                                                                   tGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCtctc                  <  2:1716313/100‑1 (MQ=255)
                                                                                   tGTTTTTTGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCtctc                  <  1:2779594/100‑1 (MQ=255)
                                                                                   tGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCtctc                  <  1:2279600/100‑1 (MQ=255)
                                                                                   tGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCtctc                  <  2:2685949/100‑1 (MQ=255)
                                                                                   tGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCtctc                  >  2:1429269/1‑100 (MQ=255)
                                                                                   tGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCtctc                  <  2:1827784/100‑1 (MQ=255)
                                                                                     ttttttGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAg                <  2:1264047/100‑1 (MQ=255)
                                                                                     ttttttGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAg                <  2:106472/100‑1 (MQ=255)
                                                                                     ttttttGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAg                >  1:321156/1‑100 (MQ=255)
                                                                                     ttttttGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATAATGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTATCAg                >  2:2527790/1‑100 (MQ=255)
                                                                                      tttttGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGt               <  1:379456/100‑1 (MQ=255)
                                                                                       ttttGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTc              >  2:259015/1‑100 (MQ=255)
                                                                                       ttttGATCTAGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTc              <  2:1031275/100‑1 (MQ=255)
                                                                                         ttGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCtt            <  1:523419/100‑1 (MQ=255)
                                                                                         ttGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCtt            <  2:1924589/100‑1 (MQ=255)
                                                                                          tGATCTGGTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCttt           >  2:1411313/1‑100 (MQ=255)
                                                                                          tGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCtt            >  1:349140/1‑99 (MQ=255)
                                                                                           gATCTGGTGTAACAGGCATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTc          >  2:1835863/1‑100 (MQ=255)
                                                                                           gATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTc          <  2:524736/100‑1 (MQ=255)
                                                                                            aTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTcc         <  1:1601958/100‑1 (MQ=255)
                                                                                            aTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTcc         <  2:1027481/100‑1 (MQ=255)
                                                                                            aTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTcc         <  1:1095686/100‑1 (MQ=255)
                                                                                                ggTGTAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCAt     <  2:900976/100‑1 (MQ=255)
                                                                                                ggTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCAt     <  2:1930431/100‑1 (MQ=255)
                                                                                                ggTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCAt     >  1:2056037/1‑100 (MQ=255)
                                                                                                   ataACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTg  <  1:2408473/100‑1 (MQ=255)
                                                                                                   ataACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTg  <  1:587662/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTG  >  NC_000913/4503528‑4503726

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: