Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 46,253 G→C A149A (GCG→GCC yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,2530GC96.3% 168.3 / ‑2.8 54A149A (GCG→GCCyaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base G (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 = 8.42e-01

TGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGACCTGTGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAG  >  NC_000913/46010‑46484
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       
tggaccaTCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCAGTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTgg                                                                                                                                                                                                                                   <  2:11752/245‑1 (MQ=255)
  gCACATCGCTGTTTGGCTATATTACCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTGCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTAt                                                                                                                                                                                                                                <  2:199566/251‑1 (MQ=255)
         gtggTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCAGCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGgcg                                                                                                                                                                                                                         <  2:363467/248‑1 (MQ=255)
                     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                                                                                                                               <  1:326425/214‑1 (MQ=255)
                                                                                                                                        gTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCtgtttgt                                                                                                                               >  2:326425/1‑214 (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)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       
TGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGACCTGTGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAG  >  NC_000913/46010‑46484

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: