Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,522,716 T→A 6.2% D117V (GAC→GTC)  yncD ← putative iron outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,522,7160TA 6.2% 389.3 / 13.4 112D117V (GAC→GTC) yncDputative iron outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (54/51);  new base A (3/4);  total (57/55)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01

Only 100 of 113 total aligned reads displayed.

TTCAATGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCA  >  NC_000913/1522568‑1522863
                                                                                                                                                    |                                                                                                                                                   
ttCAATTGTTGGTCGCTGTTGTCCGGTCTGGGTGGTGACATTCATTAACCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGTGCCACACAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGt                                                                                                                                                     >  4:583644/1‑149 (MQ=255)
   aaTGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCgg                                                                                                                                                  >  4:274607/1‑149 (MQ=255)
    aTGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCgg                                                                                                                                                  >  2:823972/1‑148 (MQ=255)
     tggttggtGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCggg                                                                                                                                                 <  1:592949/148‑1 (MQ=255)
        ttggtGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCAtg                                                                                                                                             <  3:700183/149‑1 (MQ=255)
             ggCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCAtggtgg                                                                                                                                         >  1:165691/1‑148 (MQ=255)
               cTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCgg                                                                                                                                      <  4:134282/149‑1 (MQ=255)
               cTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCg                                                                                                                                       >  2:459085/1‑148 (MQ=255)
                tgttgtCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCggg                                                                                                                                     >  3:98708/1‑149 (MQ=255)
                 gttgtCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCggg                                                                                                                                     >  1:514192/1‑148 (MQ=255)
                  ttgtCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTTCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCGCTTGCCCGTCGGGCATGGTGGCGGGaa                                                                                                                                   <  3:861919/149‑1 (MQ=255)
                    gtCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATa                                                                                                                                 >  3:567740/1‑149 (MQ=255)
                      ccGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATAc                                                                                                                                >  2:398354/1‑148 (MQ=255)
                       cGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCg                                                                                                                              >  3:648041/1‑149 (MQ=255)
                         gTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTc                                                                                                                            <  4:72526/149‑1 (MQ=255)
                             gggTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCcaca                                                                                                                        >  4:586076/1‑149 (MQ=255)
                              ggtggtGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACAt                                                                                                                       >  4:295113/1‑149 (MQ=255)
                              ggtggtGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGACGGGCaag                                                                                                                                             <  1:479273/127‑3 (MQ=255)
                              ggtggtGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGACGGGCaag                                                                                                                                             >  2:479273/1‑125 (MQ=255)
                                  gtgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCCCACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATATCGTCCACATACat                                                                                                                   <  4:584754/149‑2 (MQ=255)
                                  gtgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCAtggt                                                                                                                                           >  2:90799/1‑125 (MQ=255)
                                  gtgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCAtggt                                                                                                                                           <  1:90799/125‑1 (MQ=255)
                                  gtgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACa                                                                                                                    <  2:797465/148‑1 (MQ=255)
                                   tgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAgg                                                                                                                  <  3:274607/149‑1 (MQ=255)
                                     aCATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgc                                                                                                                 >  1:21610/1‑148 (MQ=255)
                                      cattcattACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcg                                                                                                                >  1:806780/1‑148 (MQ=255)
                                       attcattACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGTGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCg                                                                                                                                       <  4:722004/124‑1 (MQ=255)
                                       attcattACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCg                                                                                                                                       >  3:722004/1‑124 (MQ=255)
                                       attcattACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaa                                                                                                              >  4:147896/1‑149 (MQ=255)
                                        ttcattACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaat                                                                                                             >  3:464924/1‑149 (MQ=255)
                                            ttACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccg                                                                                                         <  4:24020/149‑1 (MQ=255)
                                            ttACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaatacc                                                                                                          >  1:475768/1‑148 (MQ=255)
                                             tACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccg                                                                                                         <  2:806780/148‑1 (MQ=255)
                                             tACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgagaataccgc                                                                                                        <  3:263606/149‑1 (MQ=255)
                                                cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTc                                                                                                                                                    <  1:906731/102‑1 (MQ=255)
                                                cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTc                                                                                                                                                    >  2:906731/1‑102 (MQ=255)
                                                cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaa                                                                                                     >  3:532063/1‑149 (MQ=255)
                                                  caccaGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaat                                                                                                    >  2:901175/1‑148 (MQ=255)
                                                        aCGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTa                                                                                              <  1:454114/148‑1 (MQ=255)
                                                        aCGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTa                                                                                              <  1:823972/148‑1 (MQ=255)
                                                        aCGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTa                                                                                              >  2:437702/1‑148 (MQ=255)
                                                        aCGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATAACGTaa                                                                                             >  4:101565/1‑149 (MQ=255)
                                                           cgTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGt                                                                                           <  2:514192/148‑1 (MQ=255)
                                                             ttGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTgg                                                                                         >  2:906856/1‑148 (MQ=255)
                                                             ttGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGa                                                                                        <  3:641006/149‑1 (MQ=255)
                                                             ttGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGa                                                                                        <  4:382233/149‑1 (MQ=255)
                                                              tGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAg                                                                                       <  3:133180/149‑1 (MQ=255)
                                                                ccATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGc                                                                                      <  1:107056/148‑1 (MQ=255)
                                                                ccATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCg                                                                                     >  4:821623/1‑149 (MQ=255)
                                                                 cATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCgg                                                                                    >  4:708557/1‑149 (MQ=255)
                                                                   tACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCggg                                                                                   <  1:906856/148‑1 (MQ=255)
                                                                   tACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGa                                                                                  <  3:101565/149‑1 (MQ=255)
                                                                   tACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGa                                                                                  <  3:20938/149‑1 (MQ=255)
                                                                     cAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAg                                                                                 <  1:22693/148‑1 (MQ=255)
                                                                        ggCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACAGCGAATACCGTAAGTGGAGCGGGAGCCaa                                                                             <  4:857387/149‑1 (MQ=255)
                                                                              gaAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACa                                                                                                                    >  3:324111/1‑104 (MQ=255)
                                                                              gaAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACa                                                                                                                    <  4:324111/104‑1 (MQ=255)
                                                                              gaAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCcgc                                                                       <  4:609970/149‑1 (MQ=255)
                                                                                   ggCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAAt                                                                   <  1:730152/148‑1 (MQ=255)
                                                                                     ccACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGa                                                                >  3:504459/1‑149 (MQ=255)
                                                                                      cACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGa                                                                <  1:208578/148‑1 (MQ=255)
                                                                                        cGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAg                                                             >  4:351900/1‑149 (MQ=255)
                                                                                        cGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGACGGGCAaggg                                                                                                                                           <  2:194502/71‑2 (MQ=255)
                                                                                          cagcaCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGc                                                            <  1:859171/148‑1 (MQ=255)
                                                                                          cagcaCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGc                                                            >  1:229900/1‑148 (MQ=255)
                                                                                           agcaCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTg                                                          >  4:518667/1‑149 (MQ=255)
                                                                                              aCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTaa                                                       >  3:805617/1‑149 (MQ=255)
                                                                                                  ccACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCt                                                   >  4:423206/1‑149 (MQ=255)
                                                                                                   cacaTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCtt                                                  >  4:802315/1‑149 (MQ=255)
                                                                                                   cacaTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCtt                                                  <  3:884075/149‑1 (MQ=255)
                                                                                                    acaTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGCCCACATACAGGCGAATACCGCGAATACCGTAAGTGGGGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTg                                                 <  3:708557/149‑1 (MQ=255)
                                                                                                     caTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTg                                                 <  1:693424/148‑1 (MQ=255)
                                                                                                     caTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTg                                                 <  1:459085/148‑1 (MQ=255)
                                                                                                      aTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTgc                                                >  1:561109/1‑148 (MQ=255)
                                                                                                        tttGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCg                                                                                     >  2:329897/1‑109 (MQ=255)
                                                                                                        tttGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCg                                                                                     <  1:329897/109‑1 (MQ=255)
                                                                                                             acacTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGtt                                        >  4:660636/1‑149 (MQ=255)
                                                                                                               acTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTc                                       <  2:487204/148‑1 (MQ=255)
                                                                                                                     tAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCggg                                                                                   >  2:451103/1‑98 (MQ=255)
                                                                                                                     tAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCggg                                                                                   <  1:451103/98‑1 (MQ=255)
                                                                                                                     tAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCgg                                 <  1:437702/148‑1 (MQ=255)
                                                                                                                     tAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCgg                                 >  1:449368/1‑148 (MQ=255)
                                                                                                                       aaTCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtt                               <  1:398354/148‑1 (MQ=255)
                                                                                                                        aTCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtttt                             >  4:608475/1‑149 (MQ=255)
                                                                                                                        aTCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtttt                             >  3:50276/1‑149 (MQ=255)
                                                                                                                         tCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTg                            <  4:805617/149‑1 (MQ=255)
                                                                                                                              gTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTAcc                        <  2:165691/148‑1 (MQ=255)
                                                                                                                                tGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTg                            >  1:220048/1‑142 (MQ=255)
                                                                                                                                 ggATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCa                    >  4:283411/1‑149 (MQ=255)
                                                                                                                                  gATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCa                    >  1:608311/1‑148 (MQ=255)
                                                                                                                                      tttGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCa               <  3:583644/149‑1 (MQ=255)
                                                                                                                                          cccttgcccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGc            <  2:561109/148‑1 (MQ=255)
                                                                                                                                          accatgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATGGTTc                                       >  2:393855/5‑121 (MQ=255)
                                                                                                                                           ccttgcccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCAc          <  3:351900/149‑1 (MQ=255)
                                                                                                                                           ccatgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCa           >  1:802825/4‑148 (MQ=255)
                                                                                                                                              tgcccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGc        <  2:449368/148‑1 (MQ=255)
                                                                                                                                                 ccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCgg     >  2:610628/1‑148 (MQ=255)
                                                                                                                                                  cGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTc   <  3:586076/149‑1 (MQ=255)
                                                                                                                                                   gTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCGAATCAGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCa  >  4:543958/1‑149 (MQ=255)
                                                                                                                                                   gTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTc   >  2:222563/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCAATGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCA  >  NC_000913/1522568‑1522863

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