Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 328,377 C→G 7.3% L119L (CTG→CTC) ‡ betB ← betaine aldehyde dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913328,3770CG 7.3% 564.5 / 28.8 164L119L (CTG→CTC) ‡betBbetaine aldehyde dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base C (82/70);  new base G (8/4);  total (90/74)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 165 total aligned reads displayed.

CGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATAC  >  NC_000913/328230‑328525
                                                                                                                                                   |                                                                                                                                                    
cgccAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCa                                                                                                                                                     <  2:943773/149‑1 (MQ=255)
cgccAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCa                                                                                                                                                     <  1:817853/149‑1 (MQ=255)
cgccAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCa                                                                                                                                                     <  2:741799/149‑1 (MQ=255)
   cAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGcc                                                                                                                                                  >  4:517427/1‑149 (MQ=255)
     gcgcCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCg                                                                                                                                                >  2:628722/1‑149 (MQ=255)
      cgcCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCgg                                                                                                                                               >  2:365199/1‑149 (MQ=255)
          ggggCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCgta                                                                                                                                           <  2:168759/149‑1 (MQ=255)
          ggggCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCgta                                                                                                                                           <  3:748957/149‑1 (MQ=255)
           gggCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCgtag                                                                                                                                          >  4:1524185/1‑149 (MQ=255)
             gCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCgtagta                                                                                                                                        >  1:295881/1‑149 (MQ=255)
               ggATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCCGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACt                                                                                                                                      >  1:975487/1‑149 (MQ=255)
                    tCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGc                                                                                                                                 <  2:1415217/149‑1 (MQ=255)
                      cacaGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGCTGATCCCGGCgct                                                                                                                                           >  3:686646/1‑135 (MQ=255)
                      cacaGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCAc                                                                                                                               >  2:1491137/1‑149 (MQ=255)
                                 cTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCgg                                                                                                                    >  4:1403408/1‑149 (MQ=255)
                                   ggATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTa                                                                                                                  >  4:1298859/1‑149 (MQ=255)
                                    gATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTaa                                                                                                                 >  2:465334/1‑149 (MQ=255)
                                     aTCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAAc                                                                                                                <  1:365199/149‑1 (MQ=255)
                                     aTCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAAc                                                                                                                <  4:577286/149‑1 (MQ=255)
                                        gggTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTTGTACTCCAGCACGTCCGCACCGGTAACGat                                                                                                             >  4:634582/1‑149 (MQ=255)
                                        gggTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGat                                                                                                             >  3:1345972/1‑149 (MQ=255)
                                           tAGTTCCAAGCGCCAATCCCTGCCACTACGCCCAGCTGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCTGGGATCAGCCCGGCGTAGTATTCCAGCACGTCCGCACCGGTAACGATATc                                                                                                          <  2:797548/149‑1 (MQ=255)
                                             gTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGGTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCGGCCCGGCGTAGTACTCCAGCACGGCCGCACCGGTAACGATATcaa                                                                                                        >  2:774925/1‑147 (MQ=255)
                                              ttCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCCCCGGTAACGATATCGAc                                                                                                       >  3:624554/1‑149 (MQ=255)
                                                ccATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACgg                                                                                                     <  1:1538422/149‑1 (MQ=255)
                                                 cATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGCTGATCCCGGCgct                                                                                                                                           <  4:251031/110‑3 (MQ=255)
                                                 cATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGCTGATCCCGGCgct                                                                                                                                           >  3:251031/1‑108 (MQ=255)
                                                  aTGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGtt                                                                                                   <  1:1491137/149‑1 (MQ=255)
                                                    gcgcCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGa                                                                                                 <  4:338471/149‑1 (MQ=255)
                                                     cgcCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAg                                                                                                <  3:94046/149‑1 (MQ=255)
                                                     cgcCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAg                                                                                                <  2:106020/149‑1 (MQ=255)
                                                     cgcCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAg                                                                                                <  3:1298859/149‑1 (MQ=255)
                                                      gcCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAgg                                                                                               >  1:238243/1‑149 (MQ=255)
                                                       ccAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGt                                                                                              >  3:642602/1‑149 (MQ=255)
                                                        cAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACgg                                                                                                     >  1:1001720/1‑141 (MQ=255)
                                                        cAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAgg                                                                                               <  4:45540/147‑1 (MQ=255)
                                                        cAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAgg                                                                                               >  3:45540/1‑147 (MQ=255)
                                                        cAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGtt                                                                                             >  1:440470/1‑149 (MQ=255)
                                                        cAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGt                                                                                              >  1:1499183/1‑148 (MQ=255)
                                                          aTCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTc                                                                                           >  1:686287/1‑149 (MQ=255)
                                                          aTCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGACGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTc                                                                                           <  1:628722/149‑1 (MQ=255)
                                                              cTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGaa                                                                                       <  3:1403408/149‑1 (MQ=255)
                                                              cTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGaa                                                                                       >  2:27588/1‑149 (MQ=255)
                                                               tGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAat                                                                                      >  2:325664/1‑149 (MQ=255)
                                                               tGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAat                                                                                      <  3:517427/149‑1 (MQ=255)
                                                               tGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAat                                                                                      >  2:1281953/1‑149 (MQ=255)
                                                                 ccACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCCGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAatat                                                                                    <  2:975487/149‑1 (MQ=255)
                                                                   aCTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGc                                                                                  >  4:563764/1‑149 (MQ=255)
                                                                   aCTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGc                                                                                  >  3:692764/1‑149 (MQ=255)
                                                                        gCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTc                                                                             >  2:908716/1‑149 (MQ=255)
                                                                         ctcAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTcc                                                                            <  4:1151380/147‑1 (MQ=255)
                                                                          ccAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCg                                                                           >  2:1551649/1‑149 (MQ=255)
                                                                           cAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCgg                                                                          <  3:1524185/149‑1 (MQ=255)
                                                                            aGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGgt                                                                         <  4:1152041/149‑1 (MQ=255)
                                                                                  tCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAg                                                                   <  4:792630/149‑1 (MQ=255)
                                                                                   cgcgGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAgg                                                                  <  2:1174831/149‑1 (MQ=255)
                                                                                      ggcggGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGtt                                                               <  3:1539452/149‑1 (MQ=255)
                                                                                          ggTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCa                                                           <  2:686287/149‑1 (MQ=255)
                                                                                            tataCACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGt                                                         >  1:1283370/1‑149 (MQ=255)
                                                                                            tataCACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGt                                                         >  1:1013942/1‑149 (MQ=255)
                                                                                              taCACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAg                                                                   <  3:166622/137‑1 (MQ=255)
                                                                                              taCACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAg                                                                   >  4:166622/1‑137 (MQ=255)
                                                                                                 acaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCg                                                                                                                     <  4:191768/84‑1 (MQ=255)
                                                                                                 acaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGc                                                    >  1:1588263/1‑149 (MQ=255)
                                                                                                 acaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGc                                                    >  3:681538/1‑149 (MQ=255)
                                                                                                 acaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGc                                                    >  4:1265035/1‑149 (MQ=255)
                                                                                                   aaaGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGa                                                  >  3:898497/1‑149 (MQ=255)
                                                                                                       gACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGgt                                                                         <  2:875503/122‑1 (MQ=255)
                                                                                                       gACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCTATTTTTGCGAGTTc                                              <  3:116549/149‑1 (MQ=255)
                                                                                                       gACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTc                                              >  2:1230047/1‑149 (MQ=255)
                                                                                                       gACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTc                                              <  4:1528233/149‑1 (MQ=255)
                                                                                                               aCGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGtttt                                                      <  4:671766/133‑1 (MQ=255)
                                                                                                                 gCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGa                                                                                                 >  1:1360065/1‑88 (MQ=255)
                                                                                                                 gCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGc                                      >  4:227381/1‑147 (MQ=255)
                                                                                                                 gCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTGTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGt                                    >  3:141049/1‑149 (MQ=255)
                                                                                                                 gCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGGCCGCACCGGTAACGATATCGACGGTTGa                                                                                                 <  2:1360065/88‑1 (MQ=255)
                                                                                                                  cAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTc                                                                      >  2:1499390/1‑114 (MQ=255)
                                                                                                                  cAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGtt                                                               >  2:1426190/1‑121 (MQ=255)
                                                                                                                   aaCGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTCCTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAg                                                                                                <  4:909736/87‑1 (MQ=255)
                                                                                                                   aaCGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTc                                  >  2:527575/1‑149 (MQ=255)
                                                                                                                   aaCGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTc                                  >  4:836636/1‑149 (MQ=255)
                                                                                                                    aCGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGttt                                                       >  1:835968/1‑127 (MQ=255)
                                                                                                                         aTCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAg                            <  2:1040086/149‑1 (MQ=255)
                                                                                                                          tCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCa                                                           >  3:346634/1‑117 (MQ=255)
                                                                                                                           cTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGaa                          <  3:217480/149‑1 (MQ=255)
                                                                                                                                tGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGTGTTCTTCATTGCGTTCACGCAGAATATCa                     >  4:580529/1‑149 (MQ=255)
                                                                                                                                tGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCa                     <  2:998112/149‑1 (MQ=255)
                                                                                                                                    ttCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGttt                                                       <  3:1192534/111‑1 (MQ=255)
                                                                                                                                    ttCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGttt                                                       >  4:1192534/1‑111 (MQ=255)
                                                                                                                                    ttCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACgg                 >  3:904558/1‑149 (MQ=255)
                                                                                                                                    ttCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACgg                 >  1:759999/1‑149 (MQ=255)
                                                                                                                                    ttCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGCTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACgg                 >  1:566275/1‑149 (MQ=255)
                                                                                                                                    ttCCAGCCCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATTGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTGTCCTGATTTGCGAGTGCGTCATTGCGATCACGCAGAATATCAACgg                 >  4:381594/1‑149 (MQ=255)
                                                                                                                                     tCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGc                >  3:558447/1‑149 (MQ=255)
                                                                                                                                     tCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGc                <  4:721706/149‑1 (MQ=255)
                                                                                                                                        tacgcCGGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAg                                                 >  1:1212779/3‑113 (MQ=255)
                                                                                                                                        tacgcCGGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCg             >  2:1157140/3‑149 (MQ=255)
                                                                                                                                        aGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCg             >  2:82592/1‑149 (MQ=255)
                                                                                                                                          cgcCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGAc           <  2:1322395/149‑1 (MQ=255)
                                                                                                                                             cGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCa        >  4:485378/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATAC  >  NC_000913/328230‑328525

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