Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,298,438 T→A 5.1% N137Y (AAC→TAC) ‡ yraR ← putative nucleoside‑diphosphate‑sugar epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,298,4380TA 5.1% 316.9 / 4.9 97N137Y (AAC→TAC) ‡yraRputative nucleoside‑diphosphate‑sugar epimerase
Reads supporting (aligned to +/- strand):  ref base T (42/50);  new base A (2/3);  total (44/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CATGACCCGCGCCACGTCTCTCGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTCGCCGCCCGGTTAATGCGGTATCCACTACCAGCGTGT  >  NC_000913/3298289‑3298584
                                                                                                                                                     |                                                                                                                                                  
cATGACCCGCGCCACGTCTCTCGCATCAATGGATTTCAAGTTACATGGTAACAAAAGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGtt                                                                                                                                                    >  1:991165/1‑150 (MQ=255)
       cgcgCCACGTCTCTCGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCg                                                                                                                                             <  1:1408792/150‑1 (MQ=255)
       cgcgCCACGTCTCTCGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCg                                                                                                                                             <  1:1403092/150‑1 (MQ=255)
       cgcgCCACGTCTCTCGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGCCCAGAACTGGCCg                                                                                                                                             >  1:385698/1‑150 (MQ=255)
         cgcCACGTCTCTCGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGAt                                                                                                                                           <  2:1520918/150‑1 (MQ=255)
          gcCACGTCTCTCGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGAtt                                                                                                                                          <  2:52174/150‑1 (MQ=255)
          gcCACGTCTCTCGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGCCCAGAACTGGCCGaaa                                                                                                                                          <  1:539677/150‑3 (MQ=255)
          gcCACGTCTCTCGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGCCCAGAACTGGCCGaaa                                                                                                                                          <  2:385698/150‑3 (MQ=255)
            cACGTCTCTCGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTaa                                                                                                                                        >  2:1203884/1‑150 (MQ=255)
                tctctcGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGct                                                                                                                                    >  1:1261683/1‑150 (MQ=255)
                tctctcGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGct                                                                                                                                    <  1:1113893/150‑1 (MQ=255)
                 ctctcGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGctt                                                                                                                                   >  2:892863/1‑150 (MQ=255)
                   ctcGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGcttct                                                                                                                                 >  1:1312596/1‑150 (MQ=255)
                   ctcGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGcttct                                                                                                                                 >  2:871799/1‑150 (MQ=255)
                   ctcGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGAGATTAATGcttct                                                                                                                                 >  1:14000/1‑150 (MQ=255)
                      gCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTcc                                                                                                                              >  2:830129/1‑150 (MQ=255)
                       cATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACTATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGACAGTTCTGGGCGATTAATGCttcttcaa                                                                                                                             >  2:542089/1‑148 (MQ=255)
                       cATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCa                                                                                                                             <  2:1345165/150‑1 (MQ=255)
                         tCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCAtc                                                                                                                           <  2:1296283/150‑1 (MQ=255)
                          cAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCAtct                                                                                                                          <  2:1312596/150‑1 (MQ=255)
                          cAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCAtct                                                                                                                          >  2:1268961/1‑150 (MQ=255)
                          cAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCAtct                                                                                                                          <  2:924769/150‑1 (MQ=255)
                           aaTGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCAtctc                                                                                                                         >  2:194749/1‑150 (MQ=255)
                           aaTGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCAtctc                                                                                                                         >  1:238104/1‑150 (MQ=255)
                            aTGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTcc                                                                                                                        <  2:195548/150‑1 (MQ=255)
                            aTGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTcc                                                                                                                        <  2:1210939/150‑1 (MQ=255)
                            aTGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTcc                                                                                                                        <  1:1390352/150‑1 (MQ=255)
                                 tttCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCttt                                                                                                                   >  1:655789/1‑150 (MQ=255)
                                  ttCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTg                                                                                                                  >  2:675058/1‑150 (MQ=255)
                                   tCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGa                                                                                                                 <  2:655789/150‑1 (MQ=255)
                                     cAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGAcg                                                                                                               <  1:206486/150‑1 (MQ=255)
                                       gTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGAcgcg                                                                                                             <  1:194749/150‑1 (MQ=255)
                                        ttACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCgg                                                                                                            <  1:1104195/150‑1 (MQ=255)
                                        ttACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCgg                                                                                                            <  1:1109907/150‑1 (MQ=255)
                                        ttACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCgg                                                                                                            <  1:1129953/150‑1 (MQ=255)
                                         tACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGt                                                                                                           <  1:768802/150‑1 (MQ=255)
                                         tACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGt                                                                                                           <  1:1339536/150‑1 (MQ=255)
                                          aCCTGGTAACAAACGGAACAGCGTCGCAAAAAGCGTTTCCTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGAGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAAAGCTTCTACCATCTCCACTTTGACTCGGtt                                                                                                          >  1:685478/1‑150 (MQ=255)
                                          aCCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGtt                                                                                                          <  2:1454785/150‑1 (MQ=255)
                                           ccTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTta                                                                                                         <  1:520412/150‑1 (MQ=255)
                                              ggTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGGTTGCGACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAg                                                                                                      <  2:1220182/150‑1 (MQ=255)
                                                tAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGCCCAGAACTGGCCGaaat                                                                                                                                         >  1:14454/1‑110 (MQ=255)
                                                tAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGCCCAGAACTGGCCGaaat                                                                                                                                         <  2:14454/113‑4 (MQ=255)
                                                   cAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGaaaaa                                                                                                 >  1:1275872/1‑150 (MQ=255)
                                                      aCGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACgg                                                                                              >  1:968583/1‑150 (MQ=255)
                                                       ggtatcAGCGGCGCATTAAGCGTTTCGTTCTTGCGGTTTTTGGTACGATAGCCCAGTAACATCGACTGGGTTGCAATGGTAAATTTCGGACAGTTCTGGGCTATTAATGCTTCTTACATCTCGCCTTTGACGCGGTTATAGAAAAACgtc                                                                                             <  2:1437930/145‑3 (MQ=255)
                                                       cGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGc                                                                                             <  2:1267586/150‑1 (MQ=255)
                                                        ggAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCAGCTCCCCTTTGACGCGGTTATAGAAAAACGGCg                                                                                            >  2:396240/1‑150 (MQ=255)
                                                            cAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTg                                                                                        <  1:846543/150‑1 (MQ=255)
                                                             aGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTgg                                                                                       >  2:38803/1‑150 (MQ=255)
                                                                     aaaaaGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGTTTATAGAAAAACGGCGAGTGGGCATTGGc                                                                               <  1:1053041/150‑1 (MQ=255)
                                                                     aaaaaGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCAtt                                                                                  <  2:1632540/147‑1 (MQ=255)
                                                                     aaaaaGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCAtt                                                                                  >  1:1632540/1‑147 (MQ=255)
                                                                     aaaaaGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGc                                                                               <  2:968583/150‑1 (MQ=255)
                                                                        aaGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCAcc                                                                            <  1:830129/150‑1 (MQ=255)
                                                                          gCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCAAGTGGGCATTGGCACCCa                                                                          <  2:766536/150‑1 (MQ=255)
                                                                             tttCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCg                                                                       <  2:1087917/150‑1 (MQ=255)
                                                                                 gTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACt                                                                   <  2:709522/150‑1 (MQ=255)
                                                                                 gTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACt                                                                   >  1:709522/1‑150 (MQ=255)
                                                                                  ttCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTg                                                                  <  2:983778/150‑1 (MQ=255)
                                                                                     aTCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACt                                                               >  1:687394/1‑150 (MQ=255)
                                                                                     aTCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTCGCACCCATCGCACTGAca                                                               >  2:1306957/1‑149 (MQ=255)
                                                                                      tCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTa                                                              <  1:675058/150‑1 (MQ=255)
                                                                                        cGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTAcc                                                            <  1:1203884/150‑1 (MQ=255)
                                                                                         gCTGTTTGCTACGATCGCACAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCATTTCTTGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCa                                                           <  1:542089/150‑1 (MQ=255)
                                                                                            gTTTGCTACGATCGCCCACTAACATCGACGGGGGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTATTGCTTCTTCCAACTCCCCTTTGACGCGGTTATAGAAATACGGCGAGTGGGCATTGGAACCCATCCCACTGACTACCAACa                                                        >  2:71628/1‑150 (MQ=255)
                                                                                              ttGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACata                                                      <  1:931498/150‑1 (MQ=255)
                                                                                               tGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACatat                                                     >  1:426286/1‑150 (MQ=255)
                                                                                                gCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATAtg                                                    <  2:238104/150‑1 (MQ=255)
                                                                                                 cTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATAtgc                                                   >  2:773023/1‑150 (MQ=255)
                                                                                                  tACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATAtgct                                                  >  1:1500008/1‑150 (MQ=255)
                                                                                                    cGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATAtgctgc                                                >  1:503977/1‑150 (MQ=255)
                                                                                                    cGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATAtgctgc                                                >  1:922564/1‑150 (MQ=255)
                                                                                                    cGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATAtgctgc                                                <  2:155186/150‑1 (MQ=255)
                                                                                                      aTCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTgcgc                                              >  2:1085539/1‑150 (MQ=255)
                                                                                                        cGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTgcgcgc                                            <  1:773023/150‑1 (MQ=255)
                                                                                                         gCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGcc                                           >  1:1454596/1‑150 (MQ=255)
                                                                                                         gCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGcc                                           <  1:1400363/150‑1 (MQ=255)
                                                                                                         gCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGcc                                           >  1:463705/1‑150 (MQ=255)
                                                                                                          cccAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGccc                                          >  1:693301/1‑150 (MQ=255)
                                                                                                            cAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAg                                        >  1:1106586/1‑150 (MQ=255)
                                                                                                              gTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTc                                      >  2:1448736/1‑150 (MQ=255)
                                                                                                              gTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTc                                      <  1:339781/150‑1 (MQ=255)
                                                                                                                aaCATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTcgc                                    <  1:871799/150‑1 (MQ=255)
                                                                                                                  cATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTcgccg                                  >  1:613614/1‑150 (MQ=255)
                                                                                                                      gACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAg                                        >  2:922401/1‑140 (MQ=255)
                                                                                                                      gACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAg                                        <  1:922401/140‑1 (MQ=255)
                                                                                                                      gACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTCGCCGCCCg                              >  1:1021187/1‑150 (MQ=255)
                                                                                                                       aCGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTCGCCGCCCgg                             <  1:485210/150‑1 (MQ=255)
                                                                                                                          ggCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTCGCCGCCCGGTTa                          >  1:744498/1‑150 (MQ=255)
                                                                                                                                   tGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTCGCCGCCCGGTTAATGCGGTAt                 <  2:1167938/150‑1 (MQ=255)
                                                                                                                                     gTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTCGCCGCCCGGTTAATGCGGTATcc               >  2:144560/1‑150 (MQ=255)
                                                                                                                                     gTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTCGCCGCCCGGTTAATGCGGTATcc               >  2:1406982/1‑150 (MQ=255)
                                                                                                                                      tCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTCGCCGCCCGGTTAATGCGGTATCCa              >  1:492408/1‑150 (MQ=255)
                                                                                                                                        aaTTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTCGCCGCCCGGTTAATGCGGTATCCACt            <  1:930468/150‑1 (MQ=255)
                                                                                                                                         aTTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTCGCCGCCCGGTTAATGCGGTATCCACTa           <  2:627880/150‑1 (MQ=255)
                                                                                                                                                  aGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTCGCCGCCCGGTTAATGCGGTATCCACTACCAGCgtgt  <  1:1212942/150‑1 (MQ=255)
                                                                                                                                                  aGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATAGCACTGACTACCAACATATAATGCGAGCCCAGTCGACGCCCGGTTAAAGCGGTATCCACTAACAGAgtca  >  1:705742/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CATGACCCGCGCCACGTCTCTCGCATCAATGGATTTCCAGTTACCTGGTAACAAACGGAACAGCGGCGCAAAAAGCGTTTCGTTCATCCGCTGTTTGCTACGATCGCCCAGTAACATCGACGGGCGAGCAATGGTCAATTTCGGCCAGTTCTGGGCGATTAATGCTTCTTCCATCTCCCCTTTGACGCGGTTATAGAAAAACGGCGAGTGGGCATTGGCACCCATCGCACTGACTACCAACATATGCTGCGCGCCCAGTCGCCGCCCGGTTAATGCGGTATCCACTACCAGCGTGT  >  NC_000913/3298289‑3298584

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