Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 46,247 T→C I147I (ATT→ATC)  yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,2470TC96.7% 183.4 / ‑3.9 60I147I (ATT→ATC) yaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (27/31);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

GTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGACCTGTGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAAC‑TGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCC  >  NC_000913/46004‑46489
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
gTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATcgccgcca                                                                                                                                                                                                                                              <  1:209393/251‑3 (MQ=255)
    gTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCAtgtg                                                                                                                                                                                                                                          <  1:356385/251‑1 (MQ=255)
            aTCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCg                                                                                                                                                                                                                                  >  1:384570/1‑251 (MQ=255)
                    ttGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGAc                                                                                                                                                                                                                           <  1:600524/250‑1 (MQ=255)
                               tcctaTAAAGTCGGAAGGCGCAAAAGGGTGCTCATTGATATCATCGACAGCGGCGTGATAACGGTGGCGACGATGGTTGTGTCATCCCCCGCCGAACTGTTGGTGATGCGTGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCATTGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTGATCGCCGCCAGGTGGTATGTCGGCGCGACCTGCGCCGATc                                                                                                                                                                                                                <  2:120487/246‑1 (MQ=255)
                                cctgtcaaGTCGGACGGCGCAAAATGTTCCTCATTTATATCATCGCCATCGGGGTGATATCGGTGGCGACGATTTTTGTTTCATCCCCCGTCTAACTTTTGGTGATGCGGGTACGTATCGGCTTTGTCATCGGTGCAGATTATGCCATCGCCACCTCAATGATCAC‑CTAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTAGGGCGGCGCGACCTGCGACGATCt                                                                                                                                                                                                               <  2:588030/244‑1 (MQ=255)
                                     aaaGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGGGGCGGCGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTctg                                                                                                                                                                                                         >  1:408958/1‑249 (MQ=255)
                                       aGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCt                                                                                                                                                                                                       <  1:568699/251‑1 (MQ=255)
                                             aCGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCGTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGcc                                                                                                                                                                                                                    <  1:606994/232‑1 (MQ=255)
                                             aCGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGAGGAGACGAAGTTTGTTTCATACCCCGTCGAACTGGTGGTGATGCGGGTACTTATCGGCATTGTCCTCGGTGCAGATTATACCATCGCCCCCTCACTGAT‑ACACGCGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATACGCTTTATCGCCGCCCTGTGGGTTGGCGGCGCGACCTGCGcccgatcggaaaggcgtcg                                                                                                                                                                                                   >  2:606994/1‑232 (MQ=255)
                                                                attgattacagcGCCATCGGCGTGATATCGGGGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGGTGGTGATGCGGGTACTTATCGGCATTGTCAGCGGTGCAGATTATCCCAGCGCCTCCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGTTTTATCGCCGCCATGTGGTATGTCGGCGCGACCGGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAggcg                                                                                                                                                                               <  2:607070/239‑1 (MQ=255)
                                                                    atatCAACGCCAACGGGGTGATTGTGGTGGGGACGGTGGTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGGTCTTATCGGCATTGGCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTATTGGCTTTATGATGTGGAAGGCGGCTg                                                                                                                                                                           <  1:539730/250‑1 (MQ=255)
                                                                                     gTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGt                                                                                                                                                           >  1:341508/1‑249 (MQ=255)
                                                                                       gATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAgcg                                                                                                                                                       >  1:622451/1‑251 (MQ=255)
                                                                                        atatCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAgcgc                                                                                                                                                      >  1:266733/1‑251 (MQ=255)
                                                                                              gTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGACACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGAt                                                                                                                                                >  1:604232/1‑251 (MQ=255)
                                                                                               gggCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCGGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATc                                                                                                                                               <  2:478774/250‑1 (MQ=255)
                                                                                                     gtaTGTGTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTCCTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCCCCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCtgt                                                                                                                                         <  2:63072/249‑1 (MQ=255)
                                                                                                         gtttgtttCATCCCCCGTCGAACTGTCGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCtgtttgt                                                                                                                                     <  2:266733/251‑1 (MQ=255)
                                                                                                           ttgtttCATCCCCCTTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCtgtttgtt                                                                                                                                    <  1:658930/250‑1 (MQ=255)
                                                                                                             gtttCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCCATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGGCGGCTACTGGCTTTATGATGTGGACGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGttgat                                                                                                                                 >  2:50130/1‑251 (MQ=255)
                                                                                                                  aTCCCCCGTCGAACTGTTGGGGATGCGGGTACTTATCGGCATTTTCCTCGGTGCAGATTAGCCCATCGCCACCACAATGATCCC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCTGCTTTATCGCCGTCATGTTGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGttgattttga                                                                                                                            <  1:509438/251‑1 (MQ=255)
                                                                                                                        cGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATACCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTgggc                                                                                                                      >  1:237485/1‑249 (MQ=255)
                                                                                                                            gAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGAtt                                                                                                                  >  1:181573/1‑251 (MQ=255)
                                                                                                                             aaCTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTc                                                                                                                 <  1:109080/251‑1 (MQ=255)
                                                                                                                              aCTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCg                                                                                                                >  1:328457/1‑251 (MQ=255)
                                                                                                                              aCTGTTGGTGATGCGGGTACTTATCGGCATTGACATAGGTGCAGATTATCCCATCGCCACCTCAATGCTCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGAGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTCATGGCTTTATGATGTGGATGGCGGCTGGCGCTGGATACTGGGTAGCGCGGCGATCCCATGGTTGCTGATTTTGATTGGTCGAGTCg                                                                                                                >  2:329949/1‑251 (MQ=255)
                                                                                                                                   tGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACt                                                                                                            <  1:428597/250‑1 (MQ=255)
                                                                                                                                              gTACTTTTCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTcc                                                                                                <  1:162936/251‑1 (MQ=255)
                                                                                                                                               tACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTATCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGACCTGTGCCGGTCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGGAGCGCTGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGCATCTccc                                                                                               >  2:128324/1‑251 (MQ=255)
                                                                                                                                                aCTTAGCGGCATTGTCATCGGTGCAGATTATCCCAGCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTcccc                                                                                              <  1:242186/251‑1 (MQ=255)
                                                                                                                                                   tATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCt                                                                                           >  1:79230/1‑251 (MQ=255)
                                                                                                                                                    aTCGTCATTGTCATCGGTGCAGATGATCCCATCGCCACCTCAATGAGCAC‑CGAGTGCGCCAGGTCCCGTCAGCTGGCGTTGTCCATCAGCTTTACCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATTATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGGTCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTg                                                                                          <  2:237485/251‑1 (MQ=255)
                                                                                                                                                     tCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAA‑CGAGTTCTCCAGTACCCGTCCGAGGGCGTTTTACATCAGCATTATCGCCGCCATGTGGTATGTCGGCGCGAACAGCGCCGAGCTGGTCGGCTAATGGCTTTATGATGTGGAAGGCGGCAGGAGCTTGATGCGGGGTAGCGCGGGGTTGCCCTGTGTGTTGTTTTTGATTGGTAGATTCGACATGACTGAATATCCCAGCTgg                                                                                         >  2:224068/1‑251 (MQ=255)
                                                                                                                                                          aTTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGttatt                                                                                    >  1:198851/1‑251 (MQ=255)
                                                                                                                                                               cATCGGTGCAGATTATCCCAGCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTTGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCa                                                                               <  1:388127/251‑1 (MQ=255)
                                                                                                                                                                aTCGGGGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACGGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCaa                                                                              <  2:98193/251‑1 (MQ=255)
                                                                                                                                                                  cGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCaaa                                                                             <  2:181573/250‑1 (MQ=255)
                                                                                                                                                                  cGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAg                                                                            <  1:178434/251‑1 (MQ=255)
                                                                                                                                                                   ggTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCGCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAgg                                                                           <  2:341508/251‑1 (MQ=255)
                                                                                                                                                                      gCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCg                                                                        >  1:12692/1‑251 (MQ=255)
                                                                                                                                                                       cAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGaa                                                                                                     <  1:413222/221‑1 (MQ=255)
                                                                                                                                                                       cAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGAGCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTgga                                                                                                     >  2:413222/1‑219 (MQ=255)
                                                                                                                                                                        aGATTATCCCGTCGCCAACTCAATGATCAC‑AGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTTTCGCCGCACTGTGGTATGTCGGCGCGACCTGAGCCGCTCTGGTCGGCTACTGGCTTTATGCTGTGGAAGGCGGCTGGCGCTGGATGCTGGGGAGCGCGGCGTTCCCCTGTTTGTTGATTTTTATTGGTCGGTTCGCACTGCCTGACTCTCACCTCTGGTTCTTACGCCAAggcgg                                                                        >  2:616268/1‑246 (MQ=255)
                                                                                                                                                                        aGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTTTTACGCAAAGGGCg                                                                        >  2:630090/1‑249 (MQ=255)
                                                                                                                                                                         gATTATCCCATCTCCACCTCCATGTCCAC‑CGTGTTCTCCAGTTCCCGTCAGCGGGCTTTTTCCATCAGCTTTATCTCCGCCATGTGTTATGTCTTCGCGACCTGCTCCGATCTGGTCTGCTACTTGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGTTAGCGCTGCTATCCCCTGTTTTTTGATTTTGATTGTTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGt                                                                     <  1:289010/251‑1 (MQ=255)
                                                                                                                                                                         gATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGGATGTCGGCGCGACCTGCGCCGATCTGGTCGGCAACTGGCTTTATGATGTGGCAGGCGGCTGGCGCTGGAGGCGGGGTGGCGCGGCGATCCGCTGTTTGTTGATTTTGATTGGTCGATTCGCACTTCATGAACCCCCCCCCTGGttattagccacagggcgggg                                                                     >  2:142829/1‑237 (MQ=255)
                                                                                                                                                                            tATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTaa                                                                   <  1:178085/250‑1 (MQ=255)
                                                                                                                                                                                     gCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTCCTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTCGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGACCTGCCTGAATCTCCCCTCTGGTTATTACGCAAAGGGCGAGTAAAAGAGAGCgg                                                         >  2:178797/1‑250 (MQ=255)
                                                                                                                                                                                        aCCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTTGGACTGCCTGAATCCCCCCGCTGGTTATTACACAAAGGGCGGGTAAAATAGTGCGacgg                                                      >  2:4070/1‑248 (MQ=255)
                                                                                                                                                                                         ccTCAATGATCAC‑CGAGTTCTCCAGGACCCGGCAGCGGGCGTGTGCCATCAGCATTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGGAGCGCGTCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCGGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAagag                                                     <  1:289695/251‑1 (MQ=255)
                                                                                                                                                                                         ccTCAAGGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAagag                                                     <  1:259034/251‑1 (MQ=255)
                                                                                                                                                                                           tCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAag                                                       >  1:583514/1‑247 (MQ=255)
                                                                                                                                                                                            cAATGATCAC‑CGATTTATCCAGTACCCGTCAGCGGGCGTTTTCCATCATCTTTATTGCCGCGATGTGGTATGTCGGCGCGACCTGTGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCTGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAgat                                                   <  1:128324/250‑1 (MQ=255)
                                                                                                                                                                                              aTGAACAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAAGCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAgatgat                                                <  2:583514/251‑1 (MQ=255)
                                                                                                                                                                                                                     cccGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAAC‑TGTTTGGCGAACCg                         <  2:12692/251‑1 (MQ=255)
                                                                                                                                                                                                                     cccGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGAGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGACATGATGATAAAAC‑TGTTTGGAGAACCg                         >  2:469152/1‑251 (MQ=255)
                                                                                                                                                                                                                              gggCGTTTTCCATCTGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGTTCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGATATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGAT‑AAACATGTTTGGCGAACAGGAGGCtttg                >  2:542839/1‑250 (MQ=255)
                                                                                                                                                                                                                                   ttttCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTCGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAACGAGGGCGAAGAGATGATGATAAAAC‑TGTTTGGCGAACGGGTGGCTTTCgttg            >  2:218096/1‑247 (MQ=255)
                                                                                                                                                                                                                                         aTCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAAC‑TGTTTGGCGAACCGGTGGCTTTCGATGAAGagc      >  1:357437/1‑250 (MQ=255)
                                                                                                                                                                                                                                            aGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCACTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGAGTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAAC‑TGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGcc  <  1:652487/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
GTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCAC‑CGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGACCTGTGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAAC‑TGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCC  >  NC_000913/46004‑46489

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