Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,522,716 T→A 7.6% 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.6% 226.8 / 8.3 66D117V (GAC→GTC) yncDputative iron outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (34/27);  new base A (2/3);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

TTCAATGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTC  >  NC_000913/1522568‑1522862
                                                                                                                                                    |                                                                                                                                                  
ttCAATGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGt                                                                                                                                                    >  3:59320/1‑149 (MQ=255)
    aTGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCggg                                                                                                                                                >  2:427656/1‑149 (MQ=255)
         tggtggCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGACGGGCAagg                                                                                                                                           <  1:155262/149‑1 (MQ=255)
         tggtggCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGACGGGCAagg                                                                                                                                           <  4:183363/149‑1 (MQ=255)
              gCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCg                                                                                                                                      >  1:617370/1‑149 (MQ=255)
                 gttgtCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGa                                                                                                                                   >  2:357996/1‑149 (MQ=255)
                    gtCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATa                                                                                                                                <  2:209950/149‑1 (MQ=255)
                         gTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTc                                                                                                                           <  3:223918/149‑1 (MQ=255)
                         gTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTc                                                                                                                           <  2:105148/149‑1 (MQ=255)
                               gtggtgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATa                                                                                                                     >  1:584765/1‑149 (MQ=255)
                                  gtgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAg                                                                                                                  <  4:59320/149‑1 (MQ=255)
                                   tgACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAgg                                                                                                                 <  4:515499/149‑1 (MQ=255)
                                     aCATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcg                                                                                                               >  4:103919/1‑149 (MQ=255)
                                       attcattACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaa                                                                                                             >  2:294404/1‑149 (MQ=255)
                                         tcattACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaata                                                                                                           <  4:529209/149‑1 (MQ=255)
                                                cccACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATATTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaa                                                                                                    <  3:515735/149‑1 (MQ=255)
                                                  caccaGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAgg                                                                                                                 >  3:331797/1‑134 (MQ=255)
                                                  caccaGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAgg                                                                                                                 <  4:331797/134‑1 (MQ=255)
                                                  caccaGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaata                                                                                                  >  1:397205/1‑149 (MQ=255)
                                                  caccaGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaataccgcgaata                                                                                                  >  2:232791/1‑149 (MQ=255)
                                                          gcgTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGt                                                                                          >  4:131364/1‑149 (MQ=255)
                                                                 cATACAGGGCAGAGAAGGGGTCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCGCTCGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAagtggagtgg                                                                                   <  3:349747/149‑1 (MQ=255)
                                                                 cATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATAGCGTAAGTGGAGCgg                                                                                   <  3:43137/149‑1 (MQ=255)
                                                                 cATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCgg                                                                                   >  4:471121/1‑149 (MQ=255)
                                                                  aTACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCggg                                                                                  <  3:131364/149‑1 (MQ=255)
                                                                  aTACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCggg                                                                                  <  4:21168/149‑1 (MQ=255)
                                                                   tACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAg                                                                                                                  <  3:467964/116‑1 (MQ=255)
                                                                   tACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAg                                                                                                                  >  4:467964/1‑116 (MQ=255)
                                                                   tACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGa                                                                                 <  1:427656/149‑1 (MQ=255)
                                                                        ggCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCaa                                                                            >  4:540685/1‑149 (MQ=255)
                                                                         gCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCaaa                                                                           >  1:297855/1‑149 (MQ=255)
                                                                           agagAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATc                                                                         >  4:257945/1‑149 (MQ=255)
                                                                            gagaAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATcc                                                                        >  1:568711/1‑149 (MQ=255)
                                                                             agaAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCcg                                                                       >  2:203990/1‑149 (MQ=255)
                                                                             agaAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCcg                                                                       >  1:351814/1‑149 (MQ=255)
                                                                                  gggCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAAt                                                                  >  3:250356/1‑149 (MQ=255)
                                                                                   ggCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATc                                                                 >  2:555621/1‑149 (MQ=255)
                                                                                   ggCCACGCAGCACTTCCACATTGTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCTGCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATc                                                                 <  4:427318/149‑1 (MQ=255)
                                                                                     ccACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGa                                                               >  1:10774/1‑149 (MQ=255)
                                                                                            gcaCTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGt                                                        <  3:103919/149‑1 (MQ=255)
                                                                                                ttCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAAt                                                    >  1:53167/1‑149 (MQ=255)
                                                                                                  ccACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCt                                                  >  4:369287/1‑149 (MQ=255)
                                                                                                    acaTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTg                                                <  2:360063/149‑1 (MQ=255)
                                                                                                     caTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTgc                                               <  2:10774/149‑1 (MQ=255)
                                                                                                     caTTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTgc                                               <  2:397205/149‑1 (MQ=255)
                                                                                                          tGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATa                                          <  1:320216/149‑1 (MQ=255)
                                                                                                                cTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTg                                    >  2:1766/1‑149 (MQ=255)
                                                                                                                  gCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGcc                                  <  1:294404/149‑1 (MQ=255)
                                                                                                                    ttAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCgg                                >  1:545582/1‑149 (MQ=255)
                                                                                                                    ttAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCgg                                >  4:156848/1‑149 (MQ=255)
                                                                                                                      aaaTCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtt                              <  3:135810/149‑1 (MQ=255)
                                                                                                                        aTCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtttt                            >  2:230306/1‑149 (MQ=255)
                                                                                                                        aTCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGtttt                            >  4:457490/1‑149 (MQ=255)
                                                                                                                         tCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGAGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTg                           <  2:617370/149‑1 (MQ=255)
                                                                                                                           gATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTa                         >  3:426228/1‑149 (MQ=255)
                                                                                                                           gATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTa                         >  1:481399/1‑149 (MQ=255)
                                                                                                                           gATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGACGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTa                         <  4:426228/149‑1 (MQ=255)
                                                                                                                            aTGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTAc                        <  2:584765/149‑1 (MQ=255)
                                                                                                                                   aTGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCTAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCaaa                 <  3:399678/149‑1 (MQ=255)
                                                                                                                                       ttgcccttgcccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAg             >  2:249925/1‑149 (MQ=255)
                                                                                                                                       gccaccatgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAg             >  4:102411/8‑149 (MQ=255)
                                                                                                                                        tgcccttgcccGTCGGGCATGGTGGCGGGAATACCGTCcaca                                                                                                                       >  2:70794/1‑42 (MQ=255)
                                                                                                                                        tgcccttgcccGTCGGGCATGGTGGCGGGAATACCGTCcaca                                                                                                                       <  1:70794/42‑1 (MQ=255)
                                                                                                                                         caccatgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGc           >  4:406642/6‑149 (MQ=255)
                                                                                                                                             atgcccGACGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGc       <  3:102411/148‑1 (MQ=255)
                                                                                                                                                cccGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCgg    >  4:553260/1‑149 (MQ=255)
                                                                                                                                                  cGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTc  >  1:433899/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTCAATGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTC  >  NC_000913/1522568‑1522862

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: ‑