Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,988,115 C→T D318N (GAT→AAT)  hemX ← PF04375 family protein HemX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,988,1150CT86.8% 126.4 / 12.1 53D318N (GAT→AAT) hemXPF04375 family protein HemX
Reads supporting (aligned to +/- strand):  ref base C (5/2);  new base T (20/26);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

TTTAGCATGGTTATTCTCCTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCTGATGGAACTGGA  >  NC_000913/3987875‑3988361
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
tttAGCATGGTTATTCTCCTTGCGGCGCGGCTGCCGGAGTATAAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCCGCCTCGCCTGGCTTTGCAGGGATTCCGGAAGATCCATCGAGATCTTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGTGGGCATCATCAGTATTGTAGTAAGcc                                                                                                                                                                                                                                              >  1:319890/1‑250 (MQ=255)
      aTGGTTATTCTCCTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTAcc                                                                                                                                                                                                                                        >  1:414610/1‑251 (MQ=255)
             ttCTCCTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGCGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTgg                                                                                                                                                                                                                                 >  1:192481/1‑251 (MQ=255)
              tctcCTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTgg                                                                                                                                                                                                                                 >  1:175856/1‑250 (MQ=255)
                 ccTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGAc                                                                                                                                                                                                                             <  1:33043/251‑1 (MQ=255)
                 ccTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGAc                                                                                                                                                                                                                             >  1:24129/1‑251 (MQ=255)
                   ttGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGt                                                                                                                                                                                                                           >  1:345936/1‑251 (MQ=255)
                         cgcggcTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCa                                                                                                                                                                                                                     <  1:108217/251‑1 (MQ=255)
                           cggcTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAgc                                                                                                                                                                                                                   <  1:202229/251‑1 (MQ=255)
                            ggcTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGTCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCAGCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAgcg                                                                                                                                                                                                                  <  1:134250/251‑1 (MQ=255)
                              cTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGCACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGcc                                                                                                                                                                                                                >  2:361618/1‑251 (MQ=255)
                                gCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGcc                                                                                                                                                                                                                <  1:347274/249‑1 (MQ=255)
                                     aGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGAt                                                                                                                                                                                                         >  1:132597/1‑251 (MQ=255)
                                                  ggtgcaggtgcaGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTg                                                                                                                                                                                            <  1:96127/251‑1 (MQ=255)
                                                   gtgcaggtgcagGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTgg                                                                                                                                                                                           <  1:66011/251‑1 (MQ=255)
                                                      caggtgcaggtTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGCCTCTTCCTGGTGa                                                                                                                                                                                        <  2:224643/251‑1 (MQ=255)
                                                                     gttCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCAGCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtgt                                                                                                                                                                         <  1:239603/250‑1 (MQ=255)
                                                                             gTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAc                                                                                                                                                                 <  1:361618/251‑1 (MQ=255)
                                                                               ccccGCTGCCGGTTGAACAAGCGGGTTACGAACGCGATCCTCCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCCTCGGGATATTTTGCTGGCTTAACTGGTACACCTCGTCGAGGAACGCTTTGGTGTTGTCATCATCAGTATTTTAGTAAGCACGAACACAGATGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTAGTGACTCGGTACAGCTTGTGCTGCgccc                                                                                                                                                                >  2:83957/1‑247 (MQ=255)
                                                                                   gCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagc                                                                                                                                                             >  1:460454/1‑249 (MQ=255)
                                                                                      gCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCg                                                                                                                                                        >  1:357751/1‑251 (MQ=255)
                                                                                         ggTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCTGCATCTCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGTGATATTTTGCGGGCTTAACTGTTCCACCTCTTCGAGGAACTCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCATCGCCTTGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgaga                                                                                                                                                     <  1:454565/251‑1 (MQ=255)
                                                                                                   gcagGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATAttt                                                                                                                                           <  1:252500/251‑1 (MQ=255)
                                                                                                          acgcacgcGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCgcgc                                                                                                                                     >  1:239754/1‑250 (MQ=255)
                                                                                                            cctcgcgAGTCTGCCTCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCaga                                                                                                                                  <  2:412884/248‑1 (MQ=255)
                                                                                                              acgcGGGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTACGCACGTACCCAGGTGGAAACGTTCTCCAGCGCCTGGCGATAATGCTCTTCCTGGTGACGCGGTACTGCTTGTGCTGCGCCCAGCAGGCGAGAGCGAATATTTTCTCGCagat                                                                                                                                 >  2:25218/1‑250 (MQ=255)
                                                                                                                gcgAGTCTGCAGCATTTTTTCCATGATGGCCTGGCTTTGCAGGGGTTCCGGAAGATCCATCGAGAGATTTTGCGGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTTGTGGCATCATCAGTATTGTAGTAAGCACGTACTCATGTGGGGACGTTCGCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGTCGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCagataga                                                                                                                              <  1:305245/251‑1 (MQ=255)
                                                                                                                          aTCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGAtt                                                                                                                    >  1:231517/1‑251 (MQ=255)
                                                                                                                          aTCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTACGGAAGAACCATCGAGATCTTTTGCTGGCTTAACTGGTCCACATCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTATAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTACTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGCATATTTTCGCGCAGATAGATATCCTGAtt                                                                                                                    >  1:263216/1‑251 (MQ=255)
                                                                                                                             aGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCCGCAGGCGAGCGCGAATATTTTCGCGCAGATAGATATCCTGATTTgg                                                                                                                 >  1:181896/1‑251 (MQ=255)
                                                                                                                              gTTTTTCCAGCATCTCCTGGCTTTGCAGGGTTTCCGGATGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGGATTGTAGTAAGCACGTACCCAGGTTGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGgc                                                                                                                <  2:282131/251‑1 (MQ=255)
                                                                                                                                  ttCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTCAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCGTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGCGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGgcgcg                                                                                                             >  2:222579/1‑249 (MQ=255)
                                                                                                                                      aGCATCGCCTGGCTTTGCAGGGTGTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAg                                                                                                        <  1:167274/251‑1 (MQ=255)
                                                                                                                                             ccTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTAcg                                                                                                 <  2:175856/251‑1 (MQ=255)
                                                                                                                                             ccTGGCTTTGCAGGGATTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTTGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTAcg                                                                                                 <  2:33659/251‑1 (MQ=255)
                                                                                                                                                         gggTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGtcatca                                                                                     <  2:345936/251‑1 (MQ=255)
                                                                                                                                                                   aaGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGaa                                                                           <  2:231517/251‑1 (MQ=255)
                                                                                                                                                                    aGATGCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTGCTCCAGCGCCTGGCGATAAGTCTCTTCCTGTTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCATGACGGCGAAt                                                                          <  2:132597/251‑1 (MQ=255)
                                                                                                                                                                    aGATCCATCGGGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGGAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAAt                                                                          <  2:454341/251‑1 (MQ=255)
                                                                                                                                                                    aGATCCATCGAGATATTTTGCTGGCTTACCTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGgttgtatacacatggagcac                                                                          >  2:379067/1‑233 (MQ=255)
                                                                                                                                                                       tCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGGCCTCACGTCGGCGGATCg                                                                        >  1:397245/1‑250 (MQ=255)
                                                                                                                                                                        ccATCGAGATATTTTGCTGGATTAACTGGTTCACCTAGTAGAGGAACGCTTTGGTGGTGGCATCATCAGAATCGTAGTAAACCCGGACACAGGAGGAAACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGAGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGCAAGCTATCCTGATTTGGCGCTAACAGCGGTACGGCGGCGTGATCAACACGGCGaaccgg                                                                       >  1:85502/1‑246 (MQ=255)
                                                                                                                                                                        ccATCGAGATAGTTTGCTGGCTTTACTGGTCCACCGCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTAGTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTGCGCGCAGATAGAGATCCTGATTTGGCGCTAACAGCGGTACGGCTGTGTCATCACGACGGCGAATCGTa                                                                      <  1:163408/251‑1 (MQ=255)
                                                                                                                                                                              agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTTATCACCACGGCGCATcttttttacg                                                                >  1:136168/1‑242 (MQ=255)
                                                                                                                                                                                        cTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGCGGCGACGTTCTCCAGCGCCTGGAGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGATTCATCACGACcgcgcataggaacaaaactgtcccaa                                                       >  1:211160/1‑228 (MQ=255)
                                                                                                                                                                                          ggCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATTACGAGCGCGAATCGTAATGAAATTGTCCATaa                                                      >  1:178617/1‑249 (MQ=255)
                                                                                                                                                                                              tAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACGTTCCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTc                                                 <  2:460454/250‑1 (MQ=255)
                                                                                                                                                                                               aaCTTGTCCCCCTGGTCGAGGAACGCTTTGGTGGGGGCATCATCAGTATTGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTg                                               <  2:89348/251‑1 (MQ=255)
                                                                                                                                                                                               aaCTGGTCCCCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATTGTAGTAAGCACTTACCCAGGTGGAGACGTTCTACAGCGCCTGGCGATAAGTCTCTTCCTGGGGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCTAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCAAACAGCGGTACGGCGGTGTCATCACGACGGAGAATCGTAATGAAGTTGTCCATAAAGTTCt                                                <  2:397245/250‑1 (MQ=255)
                                                                                                                                                                                               aaCTGGTCAACCTAGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGGGTGATTAAGGCGGGGAATATTAATGAAGTTGTCTATAAAATTCtt                                               >  2:210381/1‑250 (MQ=255)
                                                                                                                                                                                                               cgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCGCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTGTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGa                               <  2:357751/251‑1 (MQ=255)
                                                                                                                                                                                                                                       catcaGTATTGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTCCTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCTGATGGaa       <  1:25218/251‑1 (MQ=255)
                                                                                                                                                                                                                                             tATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGGAATGAAGTTGTCCATAAAATTTTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCTGATGGAAATGGa  >  1:201526/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
TTTAGCATGGTTATTCTCCTTGCGGCGCGGCTGCCGGAGTATCAGCTTGCGGTGCAGGTGCAGGTTTAGCTTCCGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCGAATATTTTCGCGCAGATAGATATCCTGATTTGGCGCTAACAGCGGTACGGCGGTGTCATCACGACGGCGAATCGTAATGAAGTTGTCCATAAAGTTCTGCCAGCTTTTTTGCAGATTGATACGCCATTCGCTGATGGAACTGGA  >  NC_000913/3987875‑3988361

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