Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,513,111 C→G 6.1% G100A (GGC→GCC)  fecC ← ferric citrate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,513,1110CG 6.1% 349.2 / 11.9 99G100A (GGC→GCC) fecCferric citrate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (40/53);  new base G (3/3);  total (43/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

Only 100 of 101 total aligned reads displayed.

TCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGG  >  NC_000913/4512967‑4513256
                                                                                                                                                |                                                                                                                                                 
tCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGGCGCGGCGCt                                                                                                                                               <  2:429704/149‑1 (MQ=255)
  aTGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGt                                                                                                                                             <  2:788820/149‑1 (MQ=255)
  aTGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGt                                                                                                                                             <  3:747150/149‑1 (MQ=255)
  aTGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGt                                                                                                                                             <  1:444361/149‑1 (MQ=255)
   tGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGtt                                                                                                                                            >  2:78521/1‑149 (MQ=255)
      gTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAACTCACGCCGCCCCCGCAGGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTAAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGCTAAt                                                                                                                                         <  2:866554/149‑1 (MQ=255)
        aTGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGGCGCGGCGCTggctatgg                                                                                                                                       <  2:914812/149‑8 (MQ=255)
         tGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGcc                                                                                                                                      <  4:203000/149‑1 (MQ=255)
         tGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGcc                                                                                                                                      <  1:410817/149‑1 (MQ=255)
          gACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCg                                                                                                                                     >  1:797454/1‑149 (MQ=255)
          gACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCg                                                                                                                                     >  1:170321/1‑149 (MQ=255)
          gACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCg                                                                                                                                     >  1:815024/1‑149 (MQ=255)
          gACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCg                                                                                                                                     <  4:330441/149‑1 (MQ=255)
          gACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCg                                                                                                                                     >  4:445798/1‑149 (MQ=255)
           aCGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGa                                                                                                                                    <  3:788272/149‑1 (MQ=255)
           aCGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGa                                                                                                                                    <  1:749601/149‑1 (MQ=255)
            cGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAg                                                                                                                                   >  1:19307/1‑149 (MQ=255)
             gAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGc                                                                                                                                  >  3:725573/1‑149 (MQ=255)
             gAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGc                                                                                                                                  >  2:897059/1‑149 (MQ=255)
              aaaTCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCa                                                                                                                                 <  1:446347/149‑1 (MQ=255)
                aTCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGt                                                                                                                               <  3:590133/149‑1 (MQ=255)
                aTCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGt                                                                                                                               <  1:357621/149‑1 (MQ=255)
                  ccGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGa                                                                                                                             >  3:519930/1‑149 (MQ=255)
                    gCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAg                                                                                                                           <  2:418235/149‑1 (MQ=255)
                    gCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAg                                                                                                                           >  4:923920/1‑149 (MQ=255)
                        cctGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGgaga                                                                                                                       <  2:815024/149‑1 (MQ=255)
                        cctGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGgaga                                                                                                                       <  3:834025/149‑1 (MQ=255)
                         ctGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCTAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGagag                                                                                                                      <  2:501464/149‑1 (MQ=255)
                          tGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAgg                                                                                                                     >  2:860813/1‑149 (MQ=255)
                           gCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGc                                                                                                                    <  2:298519/149‑1 (MQ=255)
                               tCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTGAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCAtt                                                                                                                >  4:46637/1‑149 (MQ=255)
                               tCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCAtt                                                                                                                <  3:212454/149‑1 (MQ=255)
                               tCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCAtt                                                                                                                >  1:68826/1‑149 (MQ=255)
                               tCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCAtt                                                                                                                >  2:559488/1‑149 (MQ=255)
                               tCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCAtt                                                                                                                >  1:316721/1‑149 (MQ=255)
                                cATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTg                                                                                                               >  1:189909/1‑149 (MQ=255)
                                    aCCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGtt                                                                                                           <  4:546277/149‑1 (MQ=255)
                                       agcagcCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTtgtg                                                                                                        <  3:765899/149‑1 (MQ=255)
                                        gcagcCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTgg                                                                                                       <  1:687303/149‑1 (MQ=255)
                                           gccagcTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTc                                                                                                    >  3:263236/1‑149 (MQ=255)
                                           gccagcTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTc                                                                                                    >  3:235461/1‑149 (MQ=255)
                                            ccagcTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTtgt                                                                                                         >  2:585611/1‑143 (MQ=255)
                                            ccagcTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTtgt                                                                                                         <  1:585611/143‑1 (MQ=255)
                                              agcTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAgg                                                                                                 >  4:728664/1‑149 (MQ=255)
                                              agcTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAgg                                                                                                 <  1:454768/149‑1 (MQ=255)
                                                 tCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAAt                                                                                                                                         <  4:490024/106‑1 (MQ=255)
                                                 tCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAAt                                                                                                                                         >  3:490024/1‑106 (MQ=255)
                                                  cACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGttt                                                                                             >  1:542129/1‑149 (MQ=255)
                                                   aCGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCa                                                                                                                                 <  4:884915/112‑1 (MQ=255)
                                                   aCGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCa                                                                                                                                 >  3:884915/1‑112 (MQ=255)
                                                   aCGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTg                                                                                            <  2:96182/149‑1 (MQ=255)
                                                        gccCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTgcagca                                                                                       <  1:643465/149‑1 (MQ=255)
                                                            ccGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGt                                                                                   <  1:860813/149‑1 (MQ=255)
                                                                      cgATGAATGACAGAGAATAGCCTGCAATCGGCGTGGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCCGCGTGCCCGCGAgc                                                                         <  2:846617/149‑1 (MQ=255)
                                                                        atgaatgaCAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGtt                                                                                              >  1:873163/1‑126 (MQ=255)
                                                                        atgaatgaCAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGtt                                                                                              <  2:873163/126‑1 (MQ=255)
                                                                          gaatgaCAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGTTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCa                                                                     <  4:468820/149‑1 (MQ=255)
                                                                          gaatgaCAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCa                                                                     <  2:575740/149‑1 (MQ=255)
                                                                           aatgaCAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAg                                                                    <  1:559488/149‑1 (MQ=255)
                                                                                cAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTg                                                               >  2:72092/1‑149 (MQ=255)
                                                                                cAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTg                                                               <  2:189909/149‑1 (MQ=255)
                                                                                 agagAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTgc                                                              >  1:757815/1‑149 (MQ=255)
                                                                                  gagaATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTgcg                                                             >  3:527088/1‑149 (MQ=255)
                                                                                   agaATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTgcgc                                                            <  1:810428/149‑1 (MQ=255)
                                                                                    gaATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGcc                                                           <  3:428196/149‑1 (MQ=255)
                                                                                        aGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATc                                                       <  4:843035/149‑1 (MQ=255)
                                                                                             gCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAAc                                                  <  1:633118/149‑1 (MQ=255)
                                                                                              cAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACg                                                 <  1:78521/149‑1 (MQ=255)
                                                                                               aaTCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGt                                                                                   <  1:70144/114‑1 (MQ=255)
                                                                                               aaTCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGt                                                                                   >  2:70144/1‑114 (MQ=255)
                                                                                                aTCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGc                                               <  1:600401/149‑1 (MQ=255)
                                                                                                 tCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCg                                              >  4:250720/1‑149 (MQ=255)
                                                                                                          gACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCtt                                     >  4:727067/1‑149 (MQ=255)
                                                                                                          gACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCtt                                     >  1:350560/1‑149 (MQ=255)
                                                                                                             tCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGt                                  <  1:72092/149‑1 (MQ=255)
                                                                                                             tCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGt                                  >  4:456208/1‑149 (MQ=255)
                                                                                                             tCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGt                                  >  2:923405/1‑149 (MQ=255)
                                                                                                              cAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTg                                 <  3:741219/149‑1 (MQ=255)
                                                                                                                       ggTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGaa                        <  2:858413/149‑1 (MQ=255)
                                                                                                                       ggTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGaa                        <  1:28400/149‑1 (MQ=255)
                                                                                                                        gTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAg                       <  1:897059/149‑1 (MQ=255)
                                                                                                                        gTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAg                       <  2:19307/149‑1 (MQ=255)
                                                                                                                           aGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGaa                        <  2:692158/145‑1 (MQ=255)
                                                                                                                           aGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGaa                        >  1:692158/1‑145 (MQ=255)
                                                                                                                           aGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGtt                    >  4:355238/1‑149 (MQ=255)
                                                                                                                            gcgcCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGttt                   >  3:681946/1‑149 (MQ=255)
                                                                                                                             cgcCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGtttt                  <  4:203430/149‑1 (MQ=255)
                                                                                                                              gcCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTg                 <  2:316721/149‑1 (MQ=255)
                                                                                                                               ccATAGCCAGCGCCGCGCTGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGc                >  2:269875/1‑149 (MQ=255)
                                                                                                                                cATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCa               >  3:508483/1‑149 (MQ=255)
                                                                                                                                 aTAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCAc              <  1:233778/149‑1 (MQ=255)
                                                                                                                                  ttaacAGCGGCGCGGCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCAcc             >  3:671816/5‑149 (MQ=255)
                                                                                                                                  tAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGtttt                  <  2:441894/144‑1 (MQ=255)
                                                                                                                                  tAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGtttt                  >  1:441894/1‑144 (MQ=255)
                                                                                                                                    gCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCACGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAg           >  3:815882/1‑149 (MQ=255)
                                                                                                                                     acAGCGGCGCGGCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACgg                                                >  3:259599/2‑111 (MQ=255)
                                                                                                                                     acAGCGGCGCGGCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACgg                                                <  4:259599/110‑1 (MQ=255)
                                                                                                                                     acAGCGGCGCGGCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAgc          >  3:116764/2‑149 (MQ=255)
                                                                                                                                          gccgcgccgcTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTc     >  3:907273/1‑149 (MQ=255)
                                                                                                                                             gcgccgcTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTgg  >  4:440168/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGG  >  NC_000913/4512967‑4513256

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: