Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,039,480 T→C 5.7% V62A (GTC→GCC) ‡ cbdB → cytochrome bd‑II oxidase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,039,4800TC 5.7% 215.3 / 3.1 70V62A (GTC→GCC) ‡cbdBcytochrome bd‑II oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base C (2/2);  total (34/36)
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.63e-01

TGCTGATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGAAAAATCGCC  >  NC_000913/1039333‑1039628
                                                                                                                                                   |                                                                                                                                                    
tgctgATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTc                                                                                                                                                     >  2:439383/1‑149 (MQ=255)
    gATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGt                                                                                                                                                 <  2:582254/149‑1 (MQ=255)
    gATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGt                                                                                                                                                 <  1:41897/149‑1 (MQ=255)
      ttGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTg                                                                                                                                               <  4:671038/149‑1 (MQ=255)
       tGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgc                                                                                                                                              <  3:83855/149‑2 (MQ=255)
        ggCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGAt                                                                                                                                             <  3:625056/149‑1 (MQ=255)
             gATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCg                                                                                                                                        <  4:192777/149‑1 (MQ=255)
                 cTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCtgg                                                                                                                                    >  1:314744/1‑149 (MQ=255)
                     tggtCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTggg                                                                                                                                <  4:221480/149‑1 (MQ=255)
                     tggtCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTggg                                                                                                                                <  4:384879/149‑1 (MQ=255)
                              tGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTGATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCGGGTCTGGTTGATCCTCGCTGGTGGGGCATTAttt                                                                                                                       >  1:237981/1‑149 (MQ=255)
                              tGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTAttt                                                                                                                       >  2:111485/1‑149 (MQ=255)
                              tGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTAttt                                                                                                                       <  2:314744/149‑1 (MQ=255)
                               gATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTg                                                                                                                      >  4:35161/1‑149 (MQ=255)
                                aTCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTgc                                                                                                                     <  4:308811/149‑1 (MQ=255)
                                 tctcCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTgcc                                                                                                                    >  4:312939/1‑149 (MQ=255)
                                     cGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCt                                                                                                                <  3:675911/149‑1 (MQ=255)
                                      gACGGATTTGACATGAGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGTAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTg                                                                                                               <  2:203813/149‑1 (MQ=255)
                                       aCGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTgg                                                                                                              >  2:88848/1‑149 (MQ=255)
                                            tttGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCag                                                                                                         >  3:93738/1‑149 (MQ=255)
                                            tttGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCag                                                                                                         >  3:16452/1‑149 (MQ=255)
                                            tttGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCag                                                                                                         >  3:457832/1‑149 (MQ=255)
                                            tttGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAAGGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCag                                                                                                         >  3:363497/1‑149 (MQ=255)
                                               gACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgt                                                                                                      <  4:16452/149‑1 (MQ=255)
                                                 cATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgtgt                                                                                                    <  2:55047/149‑1 (MQ=255)
                                                  aTGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGTTAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTa                                                                                                   >  3:124884/1‑149 (MQ=255)
                                                  aTGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTa                                                                                                   <  1:1895/149‑1 (MQ=255)
                                                  aTGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTa                                                                                                   >  2:1895/1‑149 (MQ=255)
                                                  aTGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTa                                                                                                   >  4:628149/1‑149 (MQ=255)
                                                  aTGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATGATTTGCCGCCTGGCCCAGAGTGTa                                                                                                   >  4:388726/1‑149 (MQ=255)
                                                   tGGAGATCGGCTATCTCCTGCCGCTGGTCACGCGTAAAGTTAATGTACGCCGGATAGTTATAAACAGCGTTGCGACGCCCTGGTAAGGCAACCAGGTCTGGTTGATGCTCGCTGGTGGGGCATTATTTGCCGCCTAGCCCAGCGTGTAt                                                                                                  >  4:217494/1‑149 (MQ=255)
                                                                tCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGtttt                                                                                     <  4:608656/149‑1 (MQ=255)
                                                                 ctgctgCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTc                                                                                    <  2:91428/149‑1 (MQ=255)
                                                                  tgctgcCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTcc                                                                                   <  2:287072/149‑1 (MQ=255)
                                                                     tgcCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGc                                                                                <  1:111485/149‑1 (MQ=255)
                                                                      gcCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGATCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCt                                                                               <  4:557256/149‑1 (MQ=255)
                                                                              tggCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTAtgtg                                                                       <  1:507315/149‑1 (MQ=255)
                                                                              tggCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTAtgtg                                                                       <  3:257389/149‑1 (MQ=255)
                                                                                gcgcGTAATGATGATGGACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGc                                                                     >  1:229890/1‑149 (MQ=255)
                                                                                         gatgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCt                                                            >  1:278332/1‑149 (MQ=255)
                                                                                          atgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTg                                                           >  3:250391/1‑149 (MQ=255)
                                                                                          atgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTg                                                           <  3:312939/149‑1 (MQ=255)
                                                                                             atgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgtg                                                        <  3:219065/149‑1 (MQ=255)
                                                                                              tgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATACAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgtgc                                                       <  3:243122/149‑1 (MQ=255)
                                                                                                  cGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgtgctgtg                                                   <  4:363497/149‑1 (MQ=255)
                                                                                                       gATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACAAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTTTGCTCAc                                              <  3:201947/149‑1 (MQ=255)
                                                                                                        aTAGTGATAGACATCGTTGGTGCACACTGGGAAGGTAACGAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTAGCGATGATATTGGTGCTGTGCTCACt                                             <  2:347154/149‑1 (MQ=255)
                                                                                                          aGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGt                                           >  2:378141/1‑149 (MQ=255)
                                                                                                            tGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGttc                                         >  4:401691/1‑149 (MQ=255)
                                                                                                            tGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCATCGCTGGTGGGGCATTATTTGCCGCCTGGCCTAGAGTGTATGCAGCGGCGTTTTCCGGCATTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGttc                                         <  3:573078/149‑1 (MQ=255)
                                                                                                              aTAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGttctt                                       <  4:124884/149‑1 (MQ=255)
                                                                                                                aaaCAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTcc                                     <  1:162999/149‑1 (MQ=255)
                                                                                                                  aCAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTTCTCACTGTTCTTccgc                                   >  3:433774/1‑149 (MQ=255)
                                                                                                                       gTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCt                              >  3:542534/1‑149 (MQ=255)
                                                                                                                           gTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTAtgt                                                                        >  3:456796/1‑103 (MQ=255)
                                                                                                                           gTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTAtgt                                                                        <  4:456796/103‑1 (MQ=255)
                                                                                                                             gCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCtt                        >  4:447444/1‑149 (MQ=255)
                                                                                                                               acacTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTg                      >  2:308436/1‑149 (MQ=255)
                                                                                                                               acacTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTg                      >  2:386237/1‑149 (MQ=255)
                                                                                                                                cacTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGa                     >  1:677648/1‑149 (MQ=255)
                                                                                                                                cacTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGa                     <  1:386237/149‑1 (MQ=255)
                                                                                                                                 acTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGAt                    >  3:228740/1‑149 (MQ=255)
                                                                                                                                       ggatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCt                         >  2:530986/5‑138 (MQ=255)
                                                                                                                                       ggatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCt                         <  1:530986/134‑1 (MQ=255)
                                                                                                                                       ggatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATcgc              >  3:243203/5‑149 (MQ=255)
                                                                                                                                            aaCCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGaaa         >  3:20161/1‑149 (MQ=255)
                                                                                                                                                aGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGAAAAATc     >  4:333114/1‑149 (MQ=255)
                                                                                                                                                  gTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGAAAAATCGc   <  3:252506/149‑1 (MQ=255)
                                                                                                                                                  gTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGAAAAATCGc   <  1:308436/149‑1 (MQ=255)
                                                                                                                                                  gTCTGGTTGATCCTCGCTGGTGGGACATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGAAAAATCGc   >  3:177108/1‑149 (MQ=255)
                                                                                                                                                   tCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGAAAAATCGcc  <  3:405087/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGCTGATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGAAAAATCGCC  >  NC_000913/1039333‑1039628

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑