Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,523,429 A→G 6.0% K41K (AAA→AAG)  yncE → ATP‑binding protein, periplasmic, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,523,4290AG 6.0% 213.0 / 4.2 68K41K (AAA→AAG) yncEATP‑binding protein, periplasmic, function unknown
Reads supporting (aligned to +/- strand):  ref base A (28/35);  new base G (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-01

GTTCATCAAAGGGAGTCGTCATGCATTTACGTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGTGACGCAGGCGATCCATAACGATCTCAA  >  NC_000913/1523287‑1523575
                                                                                                                                              |                                                                                                                                                  
gTTCATCAAAGGGAGTCGTCATGCATTTACGTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGcc                                                                                                                                              <  1:437527/149‑1 (MQ=255)
    atcaAAGGGAGTCGTCATGCATTTACGTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTAc                                                                                                                                           <  4:216862/148‑1 (MQ=255)
    atcaAAGGGAGTCGTCATGCATTTACGTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTAc                                                                                                                                           >  3:216862/1‑148 (MQ=255)
    atcaAAGGGAGTCGTCATGCATTTACGTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACg                                                                                                                                          >  4:269131/1‑149 (MQ=255)
       aaaGGGAGTCGTCATGCATTTACGTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGaaa                                                                                                                                       <  1:182830/149‑1 (MQ=255)
       aaaGGGAGTCGTCATGCATTTACGTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGaaa                                                                                                                                       <  2:264001/149‑1 (MQ=255)
            gAGTCGTCATGCATTTACGTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCt                                                                                                                                  <  1:419403/149‑1 (MQ=255)
                cgtcATGCGTTTACTTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTAAGCAGGTCGCAGAAGAAATACTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTtata                                                                                                                              <  3:474218/149‑1 (MQ=255)
                  tcATGCATTTACGTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGc                                                                                                                            >  1:434708/1‑149 (MQ=255)
                   cATGCATTTACGTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGcc                                                                                                                           >  3:188918/1‑149 (MQ=255)
                    aTGCATTTACGTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCa                                                                                                                          <  4:92903/149‑1 (MQ=255)
                              gTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAAc                                                                                                                >  3:177503/1‑149 (MQ=255)
                                 aTCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAGTGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAAcgcg                                                                                                             >  4:45622/1‑149 (MQ=255)
                                 aTCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAGGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAAcgcg                                                                                                             >  2:278240/1‑149 (MQ=255)
                                  tCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACgcgc                                                                                                            <  3:71955/149‑1 (MQ=255)
                                  tCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACgcgc                                                                                                            >  3:535239/1‑149 (MQ=255)
                                     gTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCtgt                                                                                                         >  2:461819/1‑149 (MQ=255)
                                       tttCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTgg                                                                                                       <  1:461819/149‑1 (MQ=255)
                                           aTCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCg                                                                                                   >  4:187173/1‑149 (MQ=255)
                                            tCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGc                                                                                                  >  4:77246/1‑149 (MQ=255)
                                             cgcgCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGcc                                                                                                 >  3:253509/1‑149 (MQ=255)
                                             cgcgCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGcc                                                                                                 >  1:300943/1‑149 (MQ=255)
                                              gcgcCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCa                                                                                                <  4:71100/149‑1 (MQ=255)
                                                 ccTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACtt                                                                                             <  2:435574/149‑1 (MQ=255)
                                                  cTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTc                                                                                            <  3:516747/149‑1 (MQ=255)
                                                  cTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTc                                                                                            <  1:53238/149‑1 (MQ=255)
                                                  cTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTc                                                                                            <  3:472216/149‑1 (MQ=255)
                                                          ttCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGcc                                                                                    >  3:124914/1‑149 (MQ=255)
                                                          ttCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGcc                                                                                    >  2:221225/1‑149 (MQ=255)
                                                          ttCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCACAAAGcc                                                                                    <  1:252885/149‑1 (MQ=255)
                                                          ttCATTACTGGTAGGTGCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAGTGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGcc                                                                                    >  4:285397/1‑149 (MQ=255)
                                                                    ttAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGa                                                                          <  2:541508/149‑1 (MQ=255)
                                                                        gTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATaaa                                                                      >  2:586618/1‑149 (MQ=255)
                                                                        gTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATaaa                                                                      <  2:223546/149‑1 (MQ=255)
                                                                         ttCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAg                                                                     >  2:67852/1‑149 (MQ=255)
                                                                         ttCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAg                                                                     >  1:397230/1‑149 (MQ=255)
                                                                         ttCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAg                                                                     >  1:62234/1‑149 (MQ=255)
                                                                         ttCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAg                                                                     >  1:399143/1‑149 (MQ=255)
                                                                          tCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCtg                                                                                                          <  4:430358/111‑1 (MQ=255)
                                                                          tCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCtg                                                                                                          >  3:89191/1‑111 (MQ=255)
                                                                          tCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCtg                                                                                                          >  3:430358/1‑111 (MQ=255)
                                                                          tCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCtg                                                                                                          <  4:89191/111‑1 (MQ=255)
                                                                           cATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGt                                                                   <  3:216253/149‑1 (MQ=255)
                                                                                 ttgttgTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTg                                                             >  1:201965/1‑149 (MQ=255)
                                                                                 ttgttgTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTg                                                             >  3:240954/1‑149 (MQ=255)
                                                                                  tgttgtTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTgg                                                            >  4:177134/1‑149 (MQ=255)
                                                                                   gttgttTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGt                                                           <  1:103740/149‑1 (MQ=255)
                                                                                   gttgttTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGt                                                           <  3:471609/149‑1 (MQ=255)
                                                                                   gttgttTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGt                                                           <  3:515060/149‑1 (MQ=255)
                                                                                      gttTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAGGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTa                                                        <  3:45622/149‑1 (MQ=255)
                                                                                      gttTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAGGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTa                                                        <  1:278240/149‑1 (MQ=255)
                                                                                        ttCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATc                                                      <  4:552647/149‑1 (MQ=255)
                                                                                          catcatTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGt                                                    >  3:471453/1‑149 (MQ=255)
                                                                                           atcatTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTc                                                   >  3:164398/1‑149 (MQ=255)
                                                                                                 cAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATc                                             <  4:146535/149‑1 (MQ=255)
                                                                                                 cAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATc                                             <  4:28665/149‑1 (MQ=255)
                                                                                                 cAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGCGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATc                                             <  3:327251/149‑1 (MQ=255)
                                                                                                          ggCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCActc                                    <  4:124914/149‑1 (MQ=255)
                                                                                                                 gaagaaATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGt                             >  2:202406/1‑149 (MQ=255)
                                                                                                                  aagaaATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGTg                            >  1:490563/1‑149 (MQ=255)
                                                                                                                   agaaATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGTGa                           >  2:435060/1‑149 (MQ=255)
                                                                                                                   agaaATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGTGa                           <  2:332142/149‑1 (MQ=255)
                                                                                                                    gaaATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGTGAc                          <  4:546477/149‑1 (MQ=255)
                                                                                                                       aTGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGTGACGCa                       <  4:164398/149‑1 (MQ=255)
                                                                                                                       aTGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGTGACGCa                       <  1:435060/149‑1 (MQ=255)
                                                                                                                          ctgcGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGTGACGCAGGc                    <  1:163801/149‑1 (MQ=255)
                                                                                                                                       taggtaAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGTGACGCAGGCGATCCATAACGAt       <  2:399143/149‑1 (MQ=255)
                                                                                                                                       taggtaAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGTGACGCAGGCGATCCATAACGAt       <  1:537195/149‑1 (MQ=255)
                                                                                                                                            aaaGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGTGACGCAGGCGATCCATAACGATCTCaa  <  2:397230/149‑1 (MQ=255)
                                                                                                                                            aaaGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGTGACGCAGGCGATCCATAACGATCTCaa  <  4:188918/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GTTCATCAAAGGGAGTCGTCATGCATTTACGTCATCTGTTTTCATCGCGCCTGCGTGGTTCATTACTGTTAGGTTCATTGCTTGTTGTTTCATCATTCAGTACGCAGGCCGCAGAAGAAATGCTGCGTAAAGCGGTAGGTAAAGGTGCCTACGAAATGGCTTATAGCCAGCAAGAAAACGCGCTGTGGCTCGCCACTTCGCAAAGCCGCAAACTGGATAAAGGTGGCGTGGTTTATCGTCTTGATCCGGTCACTCTGGAAGTGACGCAGGCGATCCATAACGATCTCAA  >  NC_000913/1523287‑1523575

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑