Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,429,851 G→A 6.1% G268E (GGA→GAA)  stfR → Rac prophage; putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,8510GA 6.1% 369.3 / 9.6 116G268E (GGA→GAA) stfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (57/51);  new base A (5/2);  total (63/53)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 116 total aligned reads displayed.

ACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCT  >  NC_000913/1429706‑1429997
                                                                                                                                                 |                                                                                                                                                  
acaTCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGaaa                                                                                                                                                 >  4:60890/1‑149 (MQ=255)
 caTCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAAt                                                                                                                                                >  3:121097/1‑149 (MQ=255)
  aTCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAAtt                                                                                                                                               >  5:399253/1‑149 (MQ=255)
      gCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCAGCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCcgc                                                                                                                                           >  6:276758/1‑149 (MQ=255)
       cATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCCGCTTCCTCGGCAACGGCGGCAGGAAATTCcgc                                                                                                                                           >  6:30914/1‑148 (MQ=255)
       cATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCcgc                                                                                                                                           >  5:193043/1‑148 (MQ=255)
       cATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCcgc                                                                                                                                           >  1:395058/1‑148 (MQ=255)
       cATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCcg                                                                                                                                            <  3:186832/147‑1 (MQ=255)
       cATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGACCGGGATGCGGCGGCCTCAAAAGTAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTc                                                                                                                                              <  8:247418/145‑1 (MQ=255)
         tCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGaa                                                                                                                                        <  1:88974/149‑1 (MQ=255)
                   ccACGACGGCATCAGAAGCTGCGACCCCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCaaaa                                                                                                                              <  5:107526/149‑1 (MQ=255)
                   ccACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCaaaa                                                                                                                              <  5:276758/149‑1 (MQ=255)
                     aCGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAAc                                                                                                                            <  7:405144/149‑1 (MQ=255)
                     aCGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAAc                                                                                                                            >  1:400332/1‑149 (MQ=255)
                         aGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTcc                                                                                                                        >  7:236557/1‑149 (MQ=255)
                              tCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGAc                                                                                                                   >  1:314961/1‑149 (MQ=255)
                              tCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGAc                                                                                                                   <  2:400332/149‑1 (MQ=255)
                                 gAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGaa                                                                                                                <  5:279579/149‑1 (MQ=255)
                                 gAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGaa                                                                                                                >  4:154587/1‑149 (MQ=255)
                                     cTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGc                                                                                                             >  5:136716/1‑148 (MQ=255)
                                     cTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGc                                                                                                             >  2:371879/1‑148 (MQ=255)
                                        cGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAgg                                                                                                         <  7:21701/149‑1 (MQ=255)
                                         gACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGt                                                                                                        >  2:157841/1‑149 (MQ=255)
                                         gACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGt                                                                                                        <  2:354031/149‑1 (MQ=255)
                                         gACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGt                                                                                                        >  7:201162/1‑149 (MQ=255)
                                            cTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCTAAATCTTCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGtctt                                                                                                     <  5:396940/149‑1 (MQ=255)
                                               ggCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGtcttc                                                                                                    <  3:60890/147‑1 (MQ=255)
                                                gCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTg                                                                                                  >  6:307647/1‑148 (MQ=255)
                                                 cccGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGtcttc                                                                                                    >  6:44141/1‑145 (MQ=255)
                                                 cccGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGtcttc                                                                                                    <  5:44141/145‑1 (MQ=255)
                                                  ccGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGACCGCCAGGTCTTCTGaaa                                                                                               >  7:56691/1‑149 (MQ=255)
                                                          cggcggCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGc                                                                                        <  6:193043/148‑1 (MQ=255)
                                                          cggcggCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCg                                                                                       <  3:154587/149‑1 (MQ=255)
                                                           ggcggcCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCgg                                                                                      >  6:202058/1‑149 (MQ=255)
                                                           ggcggcCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCgg                                                                                      >  7:306615/1‑149 (MQ=255)
                                                            gcggcCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCgg                                                                                      >  8:10005/1‑148 (MQ=255)
                                                               gcCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAACACGTCCGAGACGAATGCCAGGGCATCTGAACCAGCAGCGGAAcg                                                                                  >  2:118610/1‑148 (MQ=255)
                                                                  tCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCTGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGCACGCCAGGTCTTCTGCAACGGCAGCGGGACagag                                                                               >  1:294095/1‑149 (MQ=255)
                                                                    aaaaGAGGCAGCAAAATCAACAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGTCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAgcg                                                                             >  4:38392/1‑149 (MQ=255)
                                                                    aaaaGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGtcttct                                                                                                   >  2:236623/1‑127 (MQ=255)
                                                                      aagaagCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCACCGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGCAACGGCAGCGGGACAGAGCGcc                                                                           >  7:76530/1‑149 (MQ=255)
                                                                      aagaagCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCCGGACAGAgcgc                                                                            >  5:244819/1‑148 (MQ=255)
                                                                        gaagCGGCAAAATCATCAGACTCGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGTCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTc                                                                         <  8:22142/149‑1 (MQ=255)
                                                                             agCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGccg                                                                    >  7:173611/2‑147 (MQ=255)
                                                                              gCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGc                                                                   >  4:210400/1‑149 (MQ=255)
                                                                              gCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGc                                                                   <  1:371879/149‑1 (MQ=255)
                                                                              gCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGc                                                                   >  2:250657/1‑149 (MQ=255)
                                                                              gCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGACAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACagag                                                                               >  1:1139/1‑137 (MQ=255)
                                                                              gCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGACAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACagag                                                                               <  2:1139/137‑1 (MQ=255)
                                                                                aaaaTCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCgg                                                                 >  7:78669/1‑149 (MQ=255)
                                                                                aaaaTCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCgg                                                                 >  4:113862/1‑149 (MQ=255)
                                                                                 aaaTCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGc                                                                <  3:210400/149‑1 (MQ=255)
                                                                                  aaTCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGTAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCa                                                               <  4:40696/149‑1 (MQ=255)
                                                                                    tcatcaGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTACGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAgg                                                             <  2:395058/149‑1 (MQ=255)
                                                                                       tcaGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCt                                                           >  4:80397/1‑148 (MQ=255)
                                                                                         aGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCa                                                         >  4:256978/1‑148 (MQ=255)
                                                                                         aGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCa                                                         >  3:34918/1‑148 (MQ=255)
                                                                                          gAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCaa                                                        >  2:382839/1‑148 (MQ=255)
                                                                                            aacgaacgCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAg                                                              >  2:128970/1‑140 (MQ=255)
                                                                                            aacgaacgCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAg                                                              <  1:128970/140‑1 (MQ=255)
                                                                                             acgaacgCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCaaaaa                                                     >  6:192088/1‑148 (MQ=255)
                                                                                               gaacgCATCATCCAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAg                                                  <  1:118610/149‑1 (MQ=255)
                                                                                               gaacgCATCATCACGCGCCAGTAGTGCAGCTTCCTCGGCAACTGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGCACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTGAAAAACAg                                                  >  8:352859/1‑149 (MQ=255)
                                                                                                 acgCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGATCGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCg                                                <  8:236557/149‑1 (MQ=255)
                                                                                                 acgCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCg                                                <  5:30914/149‑1 (MQ=255)
                                                                                                      tcatcaAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCgaga                                                                                                                    <  6:109030/76‑1 (MQ=255)
                                                                                                      tcatcaAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCgaga                                                                                                                    >  5:109030/1‑76 (MQ=255)
                                                                                                      gcatcaAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGc                                           <  3:310246/148‑1 (MQ=255)
                                                                                                        atcaAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGcgt                                         >  4:150288/1‑149 (MQ=255)
                                                                                                        atcaAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGcgt                                         >  4:305337/1‑149 (MQ=255)
                                                                                                        atcaAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGcgt                                         <  1:250657/149‑1 (MQ=255)
                                                                                                         tcaAGTGCCGCTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAg                                                  <  4:101190/139‑1 (MQ=255)
                                                                                                         tcaAGTGCCGCTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAg                                                  >  3:101190/1‑139 (MQ=255)
                                                                                                         tcaAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAg                                                                                                          <  5:289422/83‑1 (MQ=255)
                                                                                                         tcaAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGcgtc                                        >  5:332409/1‑149 (MQ=255)
                                                                                                          caAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGcgtcg                                       >  7:402906/1‑149 (MQ=255)
                                                                                                           aaGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGcgtcgt                                      >  2:94006/1‑149 (MQ=255)
                                                                                                             gcgcCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCt                                    <  8:306615/149‑1 (MQ=255)
                                                                                                               gcCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGc                                  <  5:192088/149‑1 (MQ=255)
                                                                                                                  agtagtGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAg                               >  7:191747/1‑149 (MQ=255)
                                                                                                                    tagtGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTg                             >  6:396429/1‑149 (MQ=255)
                                                                                                                      gtGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCa                           <  8:402906/149‑1 (MQ=255)
                                                                                                                      gtGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCa                           <  6:399253/149‑1 (MQ=255)
                                                                                                                         cAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCg                        >  3:86667/1‑149 (MQ=255)
                                                                                                                          aGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGt                       >  3:188137/1‑149 (MQ=255)
                                                                                                                          aGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCGCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGt                       <  7:10005/149‑1 (MQ=255)
                                                                                                                               ccTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTcaac                   >  5:315532/1‑148 (MQ=255)
                                                                                                                               ccTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGGcaa                    <  6:332409/147‑1 (MQ=255)
                                                                                                                                cTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGTGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAAc                                                    <  4:265299/114‑1 (MQ=255)
                                                                                                                                 tCGGCAACGGCGGCAGGAAATTCCGCGACGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACACAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCGGCGTCGTCTGCCAGTGCAGCGTCCACAAg                >  6:134240/1‑149 (MQ=255)
                                                                                                                                 tCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCg                                                >  4:339698/1‑117 (MQ=255)
                                                                                                                                 tCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAg                >  6:176573/1‑149 (MQ=255)
                                                                                                                                   ggCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGt               <  3:88689/148‑1 (MQ=255)
                                                                                                                                          gcggcAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCa       <  4:34918/149‑1 (MQ=255)
                                                                                                                                             gcAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGc    <  6:274483/149‑1 (MQ=255)
                                                                                                                                             gcAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGc    <  1:94006/149‑1 (MQ=255)
                                                                                                                                             gcAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGc    <  5:202058/149‑1 (MQ=255)
                                                                                                                                              cAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGcc   >  1:164540/1‑149 (MQ=255)
                                                                                                                                              cAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGcc   <  5:396429/149‑1 (MQ=255)
                                                                                                                                               aGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCt  <  1:157841/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCT  >  NC_000913/1429706‑1429997

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑