Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,148,972 C→T Q72* (CAG→TAG)  fabH → 3‑oxoacyl‑[acyl carrier protein] synthase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,148,9720CT97.6% 129.7 / ‑4.0 41Q72* (CAG→TAG) fabH3‑oxoacyl‑[acyl carrier protein] synthase 3
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (17/23);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

ATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACCAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTGGCGCGCACCTGCGATCCAACCGATCGTG  >  NC_000913/1148734‑1149212
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
ataACCGAAAAGTTACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCTGACAAACTCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCTTTTCAACCATGGGCTTTGAAGCGTCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGa                                                                                                                                                                                                                                       <  1:204094/251‑1 (MQ=255)
 taACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTAAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATGTGGAAAGACTAGATTGGccggtt                                                                                                                                                                                                                                      >  2:165740/1‑247 (MQ=255)
       aaaaGGGACTTAGCGTACATTTATACGAAGATTAGTGGTACTGGCAGCTATCTGACCGAACAAGTGTGGACAAACGCCGATTTGGAAAAACTGGGGGACACCTCTGACGAGTGGATTGTCACTCGTACCGTTATACGCAAACGCCACATCGCCGCGACAAACGAATCCGTGTCAACCATGGGCTTTGAAGCGGCGACACTCGTAAGTGAGGTGGTGGGCATTGAGAAAGCCTAGATTGGCTTGATAGTtgt                                                                                                                                                                                                                                <  2:160192/251‑1 (MQ=255)
              aCTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACCCCTCTGACGAGTGGATTGTCACTCGTACCGGTAGCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACg                                                                                                                                                                                                                         <  2:337176/251‑1 (MQ=255)
              aCTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACg                                                                                                                                                                                                                         >  2:83943/1‑251 (MQ=255)
                tGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGTTTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCCCATTGCCGCGCCAAACGAAACCGTTTCAACCATGGTCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGGGAAGGACTATATTGGCCTGATCGTTGTGGTAACGAc                                                                                                                                                                                                                       >  2:2179/1‑251 (MQ=255)
                   gCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGAGGGCATTGAGAAAGCCTAAATTGGCCTGATCGTTTTGGCAACGACttt                                                                                                                                                                                                                    >  2:281365/1‑250 (MQ=255)
                       aCATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGCCGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCAGCGCCAAACGAAACCGTTACACCCATGGGCTTTGAAGCGGCGACACGCGCAATTGTGCTAGCGGGCATTGAGAAAGACTTGATTGGCCTTATCGTTGTGGCAACGACTTCTGCt                                                                                                                                                                                                                >  2:215458/1‑251 (MQ=255)
                         aTGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTGTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTac                                                                                                                                                                                                              <  2:133831/251‑1 (MQ=255)
                                       tattGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGGAAGACTAGATTGGCCTGATCGTTGTGGCAAAGACTTCTGCTACGCACGCTTTCCCgg                                                                                                                                                                                                >  1:334024/1‑250 (MQ=255)
                                                 ggCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCtt                                                                                                                                                                                       <  1:112494/250‑1 (MQ=255)
                                                            gCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAATCGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTGCTGCTACGCACGCTGTCCCGAGCGCAGCTTGTCAGATTCa                                                                                                                                                                            <  1:112677/249‑1 (MQ=255)
                                                                     aGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCCCGCTTTCCCGAGCGCAGCTGTGCAGATTCAAAGCATGTTgg                                                                                                                                                                 >  2:157172/1‑251 (MQ=255)
                                                                       ggCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGGTGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGc                                                                                                                                                               <  2:227221/250‑1 (MQ=255)
                                                                                    cGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGcc                                                                                                                                                  >  2:27071/1‑251 (MQ=255)
                                                                                      aTTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCg                                                                                                                                                >  1:7547/1‑251 (MQ=255)
                                                                                           gAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCAtt                                                                                                                                           <  1:83943/251‑1 (MQ=255)
                                                                                              aaaaTGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGa                                                                                                                                        <  1:167748/251‑1 (MQ=255)
                                                                                                        acacCTCTGACGAGTGGATTGTCACGCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTgcagca                                                                                                                              <  2:7547/251‑1 (MQ=255)
                                                                                                         cacCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGTCTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGCTTCAAAGAATTTTGGGCATTAAAGGTTGCCCGGCATTTGCCGTTgcagccc                                                                                                                             >  2:295944/1‑249 (MQ=255)
                                                                                                         cacCTATGACGAGTGGATTGTCACTCGTACCGGTATCCGCGACCGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGAATAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGATTTCCCGAGCGCAGCTTGTCAGATTCCATGCATGTTGGGCATTAAAGGTTGCACGGCATTTGAAGTTgcagatg                                                                                                                             >  2:267379/1‑248 (MQ=255)
                                                                                                           ccTCTGACGAGGGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTGGACGTTGCAGCAGcc                                                                                                                           <  1:86036/251‑1 (MQ=255)
                                                                                                                        gATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCa                                                                                                              <  1:135272/251‑1 (MQ=255)
                                                                                                                                       cGGTAGCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCg                                                                                               <  1:101006/251‑1 (MQ=255)
                                                                                                                                        ggTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGt                                                                                              >  1:347544/1‑251 (MQ=255)
                                                                                                                                             ccGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCa                                                                                                                                                                                  <  2:100764/162‑1 (MQ=255)
                                                                                                                                             ccGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGcc                                                                                          <  1:157172/250‑1 (MQ=255)
                                                                                                                                                   aCGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACCAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAAt                                                                                   >  2:163707/1‑251 (MQ=255)
                                                                                                                                                            tGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACCCGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATCCGTGatat                                                                          >  2:280883/1‑248 (MQ=255)
                                                                                                                                                                  gcCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTgggg                                                                    >  1:87070/1‑251 (MQ=255)
                                                                                                                                                                  gcCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTgggg                                                                    <  1:168074/251‑1 (MQ=255)
                                                                                                                                                                        cGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGGTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGa                                                              <  1:215458/251‑1 (MQ=255)
                                                                                                                                                                              cggttCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAAGACGTGAAATCTGGGGCGGTGAAGTATg                                                        <  1:295944/248‑1 (MQ=255)
                                                                                                                                                                                      ccATGGGCATTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGGGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGGTTCGCCTATGCATTAAGCGTAGCCGATCAATACGGGAAATCTGGGGGGGGGAAGAAGGCTCTgggc                                                >  2:50878/1‑249 (MQ=255)
                                                                                                                                                                                                 gAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCg                                      <  2:87070/250‑1 (MQ=255)
                                                                                                                                                                                                      ggcgTCACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATGTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTAc                                <  1:218953/251‑1 (MQ=255)
                                                                                                                                                                                                        cgACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTg                              <  1:27071/251‑1 (MQ=255)
                                                                                                                                                                                                             cgcgCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATTTTGGGCATTAAAGGTTGCCCGGTATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTGgcgcg                         <  2:292644/251‑1 (MQ=255)
                                                                                                                                                                                                                              gggCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTGGCGCGCACCTGCGATCCAACCg        <  1:21525/251‑1 (MQ=255)
                                                                                                                                                                                                                               ggCATTGAGACAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCAAGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCACTACGTGAAATCTGGGGCGGTGGAGGATGGTCTGGTCCTCGGTTCCGATGTACTGgccggcaccctcgaacacaacg        >  2:338202/1‑230 (MQ=255)
                                                                                                                                                                                                                                     tGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTGGCGCGCACCTGCGATCCAACCGATCGTg  >  1:67078/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
ATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACCAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTGGCGCGCACCTGCGATCCAACCGATCGTG  >  NC_000913/1148734‑1149212

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