Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 46,277 T→C C157C (TGT→TGC yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,2770TC96.3% 154.7 / ‑4.1 54C157C (TGT→TGCyaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (22/30);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGACCTGTGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTT  >  NC_000913/46031‑46517
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
tatTTCCGCTAAAGTCGGACGGCGCAAAATGTTCCTCATTGTTATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCCTCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGGGTTATCCAGTACCCGCCAGCGGGCGTTTTCCATCTGCTTTATCGCCGCCCTGTGGTATGTCGGCGCGACCTGCGCCg                                                                                                                                                                                                                                              >  2:162496/1‑251 (MQ=255)
 atTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCg                                                                                                                                                                                                                                              >  2:16523/1‑250 (MQ=255)
      cGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTgg                                                                                                                                                                                                                                        >  1:277750/1‑251 (MQ=255)
       gATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGCTGTTTGTTTCATCCCCCGTCGGACTGTTGGTGATGCGGGTACTTAACGGCTTTGTCATCGGTGCAGATTATCCACTCGCCACCTCAATGATCACCGAGTTCTCCAGTCCCCGTCAGCGGGCGTTTTCCATCCGATTTGTCGCCGCCCGGTGGTATGTCGACGCGACCTGCGCCGATCTggg                                                                                                                                                                                                                                       >  2:156947/1‑250 (MQ=255)
            aGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCt                                                                                                                                                                                                                                  >  1:329079/1‑251 (MQ=255)
                 gACGGCTCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCAGCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTGAGCGGGCTTTATCCATCAGCTTTAGCGCCGCCATGTGGTGTGTCGGCGCGTCCTGCGCCGATCTGGTCGGCTACTgg                                                                                                                                                                                                                             <  1:159362/251‑1 (MQ=255)
                 gACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACTGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCTATCTGGTCGGCTACTgg                                                                                                                                                                                                                             <  1:130766/251‑1 (MQ=255)
                   cGGCGCAAAATGTTCCTCATGGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCt                                                                                                                                                                                                                           <  2:277750/251‑1 (MQ=255)
                      cgcAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCTCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCttt                                                                                                                                                                                                                         >  1:392263/1‑250 (MQ=255)
                        cAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTCCTTATCGGCATTGTCATCGGTGCCGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCTGCTTTATCGCCGCCATGTGGTATGTCGGGCCGGCCTGCGCCGATCTGGTCGGCTGCTGGCtttttt                                                                                                                                                                                                                      >  2:133016/1‑248 (MQ=255)
                        cAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTAACCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTatg                                                                                                                                                                                                                      >  1:457891/1‑251 (MQ=255)
                         aaaaTGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTatg                                                                                                                                                                                                                      >  1:74243/1‑250 (MQ=255)
                            aTGTTCCTCATTGATATCATCGCCATCGGCGTGATAGCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCATTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATTTGGTATGTCGGCGCGCCCTGCGCCGATCTGGTCGGCTACGGGCTTTATGAtgt                                                                                                                                                                                                                  <  1:425452/251‑1 (MQ=255)
                                 ccTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGCCGATGTTTGTTTCATCCCCCGTCGGACTGTTGGTGATGCGGGAACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGAGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGGATGTCGACGAGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGa                                                                                                                                                                                                               >  1:388897/1‑249 (MQ=255)
                                    cATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAgg                                                                                                                                                                                                            <  2:387194/249‑1 (MQ=255)
                                                                 cGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCgcggc                                                                                                                                                                              <  1:455822/250‑1 (MQ=255)
                                                                         acgaTGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGCTTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGCCCTGCGCCGATCTGGTc                                                                                                                                                                                                                                      >  2:210776/1‑186 (MQ=255)
                                                                         acgaTGTTTGGTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGAGCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTc                                                                                                                                                                                                                                      <  1:210776/186‑1 (MQ=255)
                                                                                 tgtttCATCCCCCGTCGCACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTttgttg                                                                                                                                                              >  1:415023/1‑250 (MQ=255)
                                                                                     tCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGttgattt                                                                                                                                                          <  1:16523/250‑1 (MQ=255)
                                                                                     tCATCCCCCGGCGAACTGTTGGTGATGCGCGTTCTTATCGGCATTGTCATCGGTGCAGATTATCCCACCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCGTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTTGGTAGCGCGGCGATCCCCTGTTTGttgattt                                                                                                                                                          <  2:420028/250‑1 (MQ=255)
                                                                                        tCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGttgattttgat                                                                                                                                                      >  2:458258/1‑251 (MQ=255)
                                                                                         cccccGTCGAACTGTTGGGGAGGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGttgattttgatt                                                                                                                                                     <  2:139549/251‑1 (MQ=255)
                                                                                           cccGTCGAACTGTTGGAGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGCTTgg                                                                                                                                                   >  1:350564/1‑251 (MQ=255)
                                                                                                                   gTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCtgtttgt                                                                                                                                                                >  2:326425/1‑214 (MQ=255)
                                                                                                                   gTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCtgtttgt                                                                                                                                                                <  1:326425/214‑1 (MQ=255)
                                                                                                                      cTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCg                                                                                                                        >  1:479399/1‑251 (MQ=255)
                                                                                                                                 tGTCAGCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGttatta                                                                                                              <  1:349428/250‑1 (MQ=255)
                                                                                                                                     aTCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGCTGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTCTTACGcc                                                                                                          >  1:234488/1‑249 (MQ=255)
                                                                                                                                        ggTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAgg                                                                                                      <  1:151204/251‑1 (MQ=255)
                                                                                                                                           gCAGATTTTCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTGCCACCAGCTTGATTGCCGTGATGTGGTATGTCGGCGCGACCTGTGCCGTTCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGc                                                                                                    <  2:137081/250‑1 (MQ=255)
                                                                                                                                             aGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTTCCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAg                                                                                                 <  1:453858/251‑1 (MQ=255)
                                                                                                                                                  aTCCCATCGCCACCTCAATGATCACCGCGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTaaaa                                                                                            <  1:408697/251‑1 (MQ=255)
                                                                                                                                                    cccATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAag                                                                                           <  1:352706/250‑1 (MQ=255)
                                                                                                                                                    cccAGCGCTACCTCAATGACCACCGTGTTCTCCAGGATCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGGTGTAGGTCGGCGCGACCTGCGCCGAACTGGGCGGCTACTGGCGTTATGATGTGGAAGGCGGCTGGCGCTGGATGGTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACCGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAaga                                                                                          <  1:7667/251‑1 (MQ=255)
                                                                                                                                                         cGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGc                                                                                      <  1:458221/250‑1 (MQ=255)
                                                                                                                                                          gCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCCAAGGCCGAGTAAAAGAGTGCg                                                                                     >  1:221670/1‑250 (MQ=255)
                                                                                                                                                            cACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTCCCATCAGCTTTATCGCCGCCAGGTGGTATGTCGGCGCGACCTGCGCCGTTCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGAGTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAag                                                                                  <  2:82061/251‑1 (MQ=255)
                                                                                                                                                              ccTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAagag                                                                                <  1:39539/251‑1 (MQ=255)
                                                                                                                                                                 aaaTGATCACCGAGTTCTCCAGTTCCCGTCAGCGGGCGTTTTCCATCAGCTGTATCGCCGCCATGTGGTATGGCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACTCAAAGGGCGAGTAAAAGAGTGCGAAGAgat                                                                              <  2:415023/249‑1 (MQ=255)
                                                                                                                                                                      aTCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCg                                                                                                                                                >  1:133161/1‑179 (MQ=255)
                                                                                                                                                                      aTCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCg                                                                                                                                                <  2:133161/179‑1 (MQ=255)
                                                                                                                                                                        cACCGAGTTCTACAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGCTGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCTCCGCTGGTTATTACGCAAGGGGCGAGAAACAGTGGGCGTAGAgatgatgata                                                                       >  2:10585/1‑249 (MQ=255)
                                                                                                                                                                                      gTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGACCTGTGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGa                                                         >  1:239655/1‑250 (MQ=255)
                                                                                                                                                                                                 cGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCtgtttg                                                                                                                                                                 <  1:89882/135‑1 (MQ=255)
                                                                                                                                                                                                 cGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCtgtttg                                                                                                                                                                 >  2:89882/1‑135 (MQ=255)
                                                                                                                                                                                                 cGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGCAGAGATGATGATaaaa                                                                    >  2:317761/1‑228 (MQ=255)
                                                                                                                                                                                                 cGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATaaaa                                                                    <  1:317761/228‑1 (MQ=255)
                                                                                                                                                                                                  gggCGTTTTCCATCAGCTTTATCGCCGGCATGTTGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTc                                            <  2:329079/251‑1 (MQ=255)
                                                                                                                                                                                                      gTTTTCCATCAGCTTTATCGCCGCCATGGGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATGGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTAGTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATg                                        <  1:162496/251‑1 (MQ=255)
                                                                                                                                                                                                           caaTCAGCTTGATCGCCGCCATGTGGTATGTCGGCGCGACCGGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAagag                                   <  1:399451/249‑1 (MQ=255)
                                                                                                                                                                                                                agctttatcgctgccctgtgGTATGGCTGCGCTACCGGCGCCGATCTGGTCGGCTACTGGCTTTATGATGGGGAATGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGATTCTCCCCGCTGGTTATTACGCAAAGGGCGATTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGATCAGcc                              <  2:333159/235‑1 (MQ=255)
                                                                                                                                                                                                                           cgccaTGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAAccc                   <  1:197826/247‑1 (MQ=255)
                                                                                                                                                                                                                                            gcgcGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGttt  <  1:160702/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
TATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGACCTGTGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTT  >  NC_000913/46031‑46517

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