Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 993,254 A→G 7.4% intergenic (+20/+23) pepN → / ← ssuB aminopeptidase N/aliphatic sulfonate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,2540AG 7.4% 307.9 / 12.5 108intergenic (+20/+23)pepN/ssuBaminopeptidase N/aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (53/47);  new base G (4/4);  total (57/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 111 total aligned reads displayed.

ACCGCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCAA  >  NC_000913/993105‑993402
                                                                                                                                                     |                                                                                                                                                    
aCCGCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCa                                                                                                                                                      >  2:1554349/1‑150 (MQ=255)
 ccGCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCaa                                                                                                                                                     >  1:146839/1‑150 (MQ=255)
     tGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAgc                                                                                                                                                 <  1:890484/150‑1 (MQ=255)
       aTTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAgcgc                                                                                                                                               >  1:298764/1‑150 (MQ=255)
       aTTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAgcgc                                                                                                                                               <  2:547259/150‑1 (MQ=255)
            ccTGAAACGTTACGATGCCAAACGTCAGGAGAAAAAGCGCGCGGCGCTGGAACAGTAGGAAGGGCTGGAAAATCTCTCTGGCTATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCa                                                                                                                                          >  1:752598/1‑150 (MQ=255)
             cTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCAt                                                                                                                                         <  1:1059002/150‑1 (MQ=255)
               gAAACGTTACGATGCCAAACGTCAGGAGAAAAAGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTc                                                                                                                                       <  2:428944/150‑1 (MQ=255)
                aaaCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTCg                                                                                                                                      <  1:1153988/150‑1 (MQ=255)
                aaaCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                                                                      >  1:1594109/1‑150 (MQ=255)
                aaaCGATACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTCg                                                                                                                                      <  1:1172267/150‑1 (MQ=255)
                  aCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTCggt                                                                                                                                    <  2:133967/150‑2 (MQ=255)
                  aCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                                                                    <  2:1687096/150‑1 (MQ=255)
                   cGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCgggg                                                                                                                                   <  1:361363/150‑1 (MQ=255)
                   cGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCgggg                                                                                                                                   >  2:455942/1‑150 (MQ=255)
                    gTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACAGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGa                                                                                                                                  >  2:1646281/1‑150 (MQ=255)
                          aTGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGTGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAc                                                                                                                            >  1:189953/1‑150 (MQ=255)
                             ccAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAaccc                                                                                                                         >  1:260130/1‑150 (MQ=255)
                                  cGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGttt                                                                                                                    <  1:1586832/150‑1 (MQ=255)
                                  cGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGttt                                                                                                                    >  2:1053471/1‑150 (MQ=255)
                                  cGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCCGGAAACTCTCTCTGGCGATCTGTACGAGAAGCTAACAAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGttt                                                                                                                    <  2:58287/150‑1 (MQ=255)
                                    tCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTtc                                                                                                                  >  2:695995/1‑150 (MQ=255)
                                       ggAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGTAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCa                                                                                                               <  2:1458635/150‑1 (MQ=255)
                                       ggAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTtctcc                                                                                                               >  2:1086592/1‑149 (MQ=255)
                                       ggAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCa                                                                                                               <  1:286233/150‑1 (MQ=255)
                                       ggAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCa                                                                                                               >  1:1482533/1‑150 (MQ=255)
                                        gagaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAg                                                                                                              <  1:1145845/150‑1 (MQ=255)
                                         agaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAgg                                                                                                             <  1:426859/150‑1 (MQ=255)
                                         agaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAgg                                                                                                             >  2:1281668/1‑150 (MQ=255)
                                         agaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAgg                                                                                                             <  2:1394637/150‑1 (MQ=255)
                                         agaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAgg                                                                                                             <  2:189953/150‑1 (MQ=255)
                                         agaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAgg                                                                                                             >  2:101843/1‑150 (MQ=255)
                                            aaaTGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGa                                                                                                          >  2:72405/1‑150 (MQ=255)
                                                        cgcTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTca                                                                                              >  2:555619/1‑150 (MQ=255)
                                                         gcTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcac                                                                                             <  1:72405/150‑1 (MQ=255)
                                                         gcTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcac                                                                                             >  2:793/1‑150 (MQ=255)
                                                         gcTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcac                                                                                             <  2:260130/150‑1 (MQ=255)
                                                          cTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcacc                                                                                            >  1:1386728/1‑150 (MQ=255)
                                                             gAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTGGATAAATTACCGAATGGCGGCAATAGAGCCGCCATTCTGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACg                                                                                         >  2:1221338/1‑150 (MQ=255)
                                                              aaCAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGt                                                                                        >  2:1685025/1‑150 (MQ=255)
                                                                cAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTg                                                                                      <  1:455942/150‑1 (MQ=255)
                                                                cAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTg                                                                                      >  1:393241/1‑150 (MQ=255)
                                                                  gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCTCTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACACCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCa                                                                                    >  1:67909/1‑150 (MQ=255)
                                                                    tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATc                                                                                  <  1:101843/150‑1 (MQ=255)
                                                                     gAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCa                                                                                 >  2:744132/1‑150 (MQ=255)
                                                                      aaaGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTATCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAc                                                                                <  1:255502/150‑1 (MQ=255)
                                                                       aaGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAcc                                                                               <  2:683683/150‑1 (MQ=255)
                                                                          ggCTGTAAAATCTCTCTGGCGATCTGTACGTGAAGATAAGTAAAGGACTGGCTTGATAAATACGCGAATCGCGGCAATAGCGCCGCCATTCGGGGTATTTTCCCCTGTTTTCTCTGGCGAGTTTCAGATTCACCACGTTGCATCACCCGt                                                                            <  1:1635844/150‑1 (MQ=255)
                                                                          ggCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCTGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACACGt                                                                            <  2:752598/150‑1 (MQ=255)
                                                                              ggAAAATCTCTCTTGCTATCTGTACGAGAAGATAACTATAGCACTGGCTTGATAACTAACCGAATGGCGGCAATAGCGCCGCCATTCTGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGATGCa                                                                        <  2:35055/150‑1 (MQ=255)
                                                                              ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGACATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCa                                                                        >  2:1103/1‑150 (MQ=255)
                                                                               gAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCaa                                                                       >  1:576293/1‑150 (MQ=255)
                                                                                 aaaTCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAaca                                                                     <  1:1002038/150‑1 (MQ=255)
                                                                                 aaaTCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAaca                                                                     <  2:1614464/150‑1 (MQ=255)
                                                                                  aaTCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAacac                                                                    <  2:373155/150‑1 (MQ=255)
                                                                                     ctctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                                                                 >  1:744656/1‑150 (MQ=255)
                                                                                        tctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGc                                                              >  1:923095/1‑150 (MQ=255)
                                                                                         ctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCt                                                             <  1:472388/150‑1 (MQ=255)
                                                                                          tGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCtt                                                            >  2:903728/1‑150 (MQ=255)
                                                                                           ggCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTc                                                           >  1:1265253/1‑150 (MQ=255)
                                                                                            gCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTcc                                                          <  2:424260/150‑1 (MQ=255)
                                                                                               aTCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGc                                                       <  1:695995/150‑1 (MQ=255)
                                                                                                 cTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTc                                                     >  1:382686/1‑150 (MQ=255)
                                                                                                    tACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGc                                                  <  2:1594109/150‑1 (MQ=255)
                                                                                                     aCGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGcg                                                 >  1:1109834/1‑149 (MQ=255)
                                                                                                     aCGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCt                                                 <  1:531252/150‑1 (MQ=255)
                                                                                                      cGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTa                                                >  2:184144/1‑150 (MQ=255)
                                                                                                        agaagaTAACTAAAGCACTGGCTTGATAAATAACCGAATCGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCATGTTTTCTCAGGCGAGTTACAGATTCACACCGTATCATCACCCGTTCCAACACCACCGCTACCAGCTCCGCAAAt                                              >  2:866978/1‑150 (MQ=255)
                                                                                                         gaagaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATc                                             <  2:1386728/150‑1 (MQ=255)
                                                                                                              tAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCAcc                                        >  1:1718103/1‑150 (MQ=255)
                                                                                                              tAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCAcc                                        >  2:919219/1‑150 (MQ=255)
                                                                                                               aaCTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCg                                       <  2:1265253/150‑1 (MQ=255)
                                                                                                               aaCTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCg                                       >  2:798841/1‑150 (MQ=255)
                                                                                                                  tAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAcc                                    >  2:1661821/1‑150 (MQ=255)
                                                                                                                   aaaGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAccc                                   >  1:1400086/1‑150 (MQ=255)
                                                                                                                         cTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCCCGTTGCCTCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTACTCGCGCCGACCCAAATCg                             >  1:314379/1‑150 (MQ=255)
                                                                                                                            gCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCCGCTCCGCTAATCGCACCGACCCTAATCGAcg                          >  2:493332/1‑150 (MQ=255)
                                                                                                                             cTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAcgc                         >  1:305619/1‑150 (MQ=255)
                                                                                                                              ttGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAcgcg                        <  1:1192094/150‑1 (MQ=255)
                                                                                                                                 ataaataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGAc                     <  2:576293/150‑1 (MQ=255)
                                                                                                                                  taaataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcaccac                                                                                          <  2:164444/80‑1 (MQ=255)
                                                                                                                                  taaataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcaccac                                                                                          >  1:164444/1‑80 (MQ=255)
                                                                                                                                  taaataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACg                    <  2:923095/150‑1 (MQ=255)
                                                                                                                                    attccccGAATGGCGGCGCTATTGCCGCCATTCGGGGTATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGCCCCTAATCGACGCGGACggg                  >  1:249040/6‑150 (MQ=255)
                                                                                                                                    aataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAg                                                        <  1:1661467/112‑1 (MQ=255)
                                                                                                                                    aataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAg                                                        >  2:1661467/1‑112 (MQ=255)
                                                                                                                                    aataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACggg                  >  2:646492/1‑150 (MQ=255)
                                                                                                                                      tccccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTACAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACggggg                >  2:516299/4‑150 (MQ=255)
                                                                                                                                        aCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACggggggt              >  1:584480/1‑148 (MQ=255)
                                                                                                                                           gAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATc           >  2:483899/1‑150 (MQ=255)
                                                                                                                                           gAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTACAGATTCAACACGTTGAATCCCCCGTTGCAACACCTCCGCTACCAGCTCCGCTAATCGCACCGACCCTAATCGAGCCGGACGGGGGAAATc           >  2:906860/1‑150 (MQ=255)
                                                                                                                                               ggcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACt       >  1:725825/1‑150 (MQ=255)
                                                                                                                                                 cggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGt     >  2:494623/1‑150 (MQ=255)
                                                                                                                                                 cggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGt     <  1:744132/150‑1 (MQ=255)
                                                                                                                                                  ggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGAATCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTc    <  2:789255/150‑1 (MQ=255)
                                                                                                                                                   gcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCa   >  1:112277/1‑150 (MQ=255)
                                                                                                                                                   gcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCa   <  2:1400086/150‑1 (MQ=255)
                                                                                                                                                   gcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCa   <  1:494623/150‑1 (MQ=255)
                                                                                                                                                   gcAATAGCGCCACCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCa   >  1:756963/1‑150 (MQ=255)
                                                                                                                                                    cAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCaa  >  2:1143827/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCGCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCAA  >  NC_000913/993105‑993402

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: