Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,522,716 T→A 7.1% 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 7.1% 669.3 / 30.2 192D117V (GAC→GTC) yncDputative iron outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (84/94);  new base A (10/4);  total (94/98)
Fisher's exact test for biased strand distribution p-value = 9.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-01

Only 100 of 193 total aligned reads displayed.

TTCAATGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCA  >  NC_000913/1522568‑1522863
                                                                                                                                                    |                                                                                                                                                   
    aTGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGCGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCAGTCggg                                                                                                                                                 >  4:449925/1‑149 (MQ=255)
       gttggtGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCATGCCCGACGGGc                                                                                                                                                <  1:392634/147‑1 (MQ=255)
         tggtggCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCAtgg                                                                                                                                            <  1:138353/149‑1 (MQ=255)
           gtggCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCAtggtg                                                                                                                                          >  4:66618/1‑149 (MQ=255)
              gCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCg                                                                                                                                       <  3:352624/149‑1 (MQ=255)
                tgttgtCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCggg                                                                                                                                     >  1:581325/1‑149 (MQ=255)
                 gttgtCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGa                                                                                                                                    <  1:830991/149‑1 (MQ=255)
                 gttgtCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGa                                                                                                                                    <  3:1263410/149‑1 (MQ=255)
                  ttgtCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGaa                                                                                                                                   <  3:1689495/149‑1 (MQ=255)
                        ggTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGt                                                                                                                             <  3:223641/149‑1 (MQ=255)
                        ggTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGt                                                                                                                             <  1:770028/149‑1 (MQ=255)
                         gTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTc                                                                                                                            <  2:715635/149‑1 (MQ=255)
                         gTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTc                                                                                                                            <  1:1518240/149‑1 (MQ=255)
                         gTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTc                                                                                                                            <  4:342561/149‑1 (MQ=255)
                           cTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCca                                                                                                                          <  1:1044255/149‑1 (MQ=255)
                                 ggtgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACa                                                                                                                    <  4:1445292/149‑1 (MQ=255)
                                 ggtgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACa                                                                                                                    <  2:1604959/149‑1 (MQ=255)
                                   tgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAgg                                                                                                                  <  4:1631595/149‑1 (MQ=255)
                                    gACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgc                                                                                                                 <  2:448935/149‑1 (MQ=255)
                                    gACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgc                                                                                                                 >  1:1175140/1‑149 (MQ=255)
                                       attcattACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATGCAGgcgaa                                                                                                              <  3:648051/149‑1 (MQ=255)
                                       attcattACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaa                                                                                                              >  2:122094/1‑149 (MQ=255)
                                             tACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgc                                                                                                        <  3:439218/149‑1 (MQ=255)
                                             tACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgc                                                                                                        <  3:66618/149‑1 (MQ=255)
                                                cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGTCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaa                                                                                                     >  4:1191334/1‑149 (MQ=255)
                                                cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTc                                                                                                                            >  1:1433340/1‑126 (MQ=255)
                                                cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTc                                                                                                                            <  2:1433340/126‑1 (MQ=255)
                                                cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaa                                                                                                     >  4:742074/1‑149 (MQ=255)
                                                cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaa                                                                                                     >  4:1224048/1‑149 (MQ=255)
                                                cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgc                                                                                                        >  1:1272280/1‑146 (MQ=255)
                                                 ccaccaGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaat                                                                                                    <  3:1177378/149‑1 (MQ=255)
                                                  caccaGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaata                                                                                                   >  3:1334620/1‑149 (MQ=255)
                                                  caccaGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaata                                                                                                   <  1:902552/149‑1 (MQ=255)
                                                      aGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGt                                                                                               >  4:506172/1‑149 (MQ=255)
                                                                 cATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCgg                                                                                    <  4:820990/149‑1 (MQ=255)
                                                                   tACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGa                                                                                  <  3:1458023/149‑1 (MQ=255)
                                                                   tACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGa                                                                                  <  3:506172/149‑1 (MQ=255)
                                                                       gggCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGCGGAGCGGGAGCCa                                                                              >  4:1360311/1‑149 (MQ=255)
                                                                        ggCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGc                                                                                      <  2:172185/140‑1 (MQ=255)
                                                                           agagAAGGGGCCAAGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATc                                                                          >  4:55307/1‑149 (MQ=255)
                                                                            gagaAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTgg                                                                                         >  4:1261353/1‑133 (MQ=255)
                                                                             agaAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaa                                                                                                              >  3:827520/1‑111 (MQ=255)
                                                                               aaGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCcgcg                                                                      >  3:1411864/1‑149 (MQ=255)
                                                                               aaGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCcgcg                                                                      >  4:928241/1‑149 (MQ=255)
                                                                               aaGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCcgcg                                                                      >  2:208749/1‑149 (MQ=255)
                                                                                   ggCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATc                                                                  <  4:1600681/149‑1 (MQ=255)
                                                                                    gCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCggg                                                                                   <  1:1622472/131‑1 (MQ=255)
                                                                                     ccACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGa                                                                <  2:1076637/149‑1 (MQ=255)
                                                                                      cACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGAc                                                               >  4:264593/1‑149 (MQ=255)
                                                                                           agcaCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGTCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGt                                                                                           <  4:855940/116‑1 (MQ=255)
                                                                                           agcaCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGTCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGt                                                                                           >  3:855940/1‑116 (MQ=255)
                                                                                              aCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTaa                                                       >  2:63208/1‑149 (MQ=255)
                                                                                                ttCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAAt                                                     >  2:238591/1‑149 (MQ=255)
                                                                                                  ccACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCgg                                                                                                                                      >  2:801684/1‑66 (MQ=255)
                                                                                                  ccACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCgg                                                                                                                                      <  1:801684/66‑1 (MQ=255)
                                                                                                    acaTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTg                                                 >  4:1655320/1‑149 (MQ=255)
                                                                                                    acaTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTg                                                 >  2:149443/1‑149 (MQ=255)
                                                                                                    acaTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTg                                                 <  4:102746/149‑1 (MQ=255)
                                                                                                     caTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTgc                                                <  2:1175140/149‑1 (MQ=255)
                                                                                                          tGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATa                                           <  1:197163/149‑1 (MQ=255)
                                                                                                               acTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTcc                                                                                                                           >  1:1024557/1‑64 (MQ=255)
                                                                                                               acTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTaa                                                                                             >  1:668729/1‑94 (MQ=255)
                                                                                                               acTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAg                                                             >  1:1204195/1‑126 (MQ=255)
                                                                                                                 tGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGc                                    >  3:286600/1‑149 (MQ=255)
                                                                                                                  gCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGcc                                   <  2:923873/149‑1 (MQ=255)
                                                                                                                    ttAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCgg                                 >  4:1264179/1‑149 (MQ=255)
                                                                                                                        aTCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtttt                             >  4:146320/1‑149 (MQ=255)
                                                                                                                        aTCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtttt                             <  2:436127/149‑1 (MQ=255)
                                                                                                                        aTCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtttt                             <  2:1357242/149‑1 (MQ=255)
                                                                                                                         tCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTg                            <  2:850731/149‑1 (MQ=255)
                                                                                                                         tCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTg                            <  1:874445/149‑1 (MQ=255)
                                                                                                                          cGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGt                           >  1:1594920/1‑149 (MQ=255)
                                                                                                                           gATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGa                                                                                  <  1:727512/93‑1 (MQ=255)
                                                                                                                           gATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtttt                             >  3:489104/1‑146 (MQ=255)
                                                                                                                           gATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtttt                             <  4:489104/146‑1 (MQ=255)
                                                                                                                           gATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTa                          >  3:1247020/1‑149 (MQ=255)
                                                                                                                            aTGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTAc                         <  2:697077/149‑1 (MQ=255)
                                                                                                                             tGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTAcc                        <  3:1601430/149‑1 (MQ=255)
                                                                                                                               ttGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCg                                                                                     >  1:757901/1‑86 (MQ=255)
                                                                                                                               ttGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCg                                                                                     <  2:757901/86‑1 (MQ=255)
                                                                                                                               ttGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTg                      >  4:1282658/1‑149 (MQ=255)
                                                                                                                                 ggATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGttt                              <  4:657616/139‑1 (MQ=255)
                                                                                                                                 ggATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCa                    >  1:1553024/1‑149 (MQ=255)
                                                                                                                                  gATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCaa                   >  1:1064885/1‑149 (MQ=255)
                                                                                                                                   aTGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCaaa                  <  2:947402/149‑1 (MQ=255)
                                                                                                                                      tttGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCcgcg                                                                      >  2:9001/1‑94 (MQ=255)
                                                                                                                                      tttGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTg                                                 >  4:1570516/1‑115 (MQ=255)
                                                                                                                                      tttGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCgg                                 >  1:710431/1‑131 (MQ=255)
                                                                                                                                       ttgcccttgcccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAg              <  2:529800/149‑1 (MQ=255)
                                                                                                                                       gccaccatgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAg              >  2:958328/8‑149 (MQ=255)
                                                                                                                                        tgcccttgcccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAgg             <  3:1282658/149‑1 (MQ=255)
                                                                                                                                         caccatgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGttt                              <  3:1583075/126‑1 (MQ=255)
                                                                                                                                         caccatgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGc            >  3:325292/6‑149 (MQ=255)
                                                                                                                                           ccttgcccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCAc          >  3:1599193/1‑149 (MQ=255)
                                                                                                                                            catgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtt                               <  4:947948/125‑1 (MQ=255)
                                                                                                                                            catgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtt                               >  3:947948/3‑127 (MQ=255)
                                                                                                                                            catgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACg         >  3:868078/3‑149 (MQ=255)
                                                                                                                                                 ccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGt    <  3:781827/149‑1 (MQ=255)
                                                                                                                                                  cGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGAc                                                               >  1:473134/1‑89 (MQ=255)
                                                                                                                                                   gTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCa  >  2:277125/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCAATGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCA  >  NC_000913/1522568‑1522863

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