Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,429,858 C→T 100% S270S (TCC→TCT 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,8580CT97.4% 4835.2 / 40.9 1711S270S (TCC→TCTstfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base T (872/787);  total (899/812)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 1871 total aligned reads displayed.

GCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGT  >  NC_000913/1429712‑1430005
                                                                                                                                                  |                                                                                                                                                   
gcctCCACCGCGGCCCCGAAAGCGTCAGAGGCCGCCACTTCAGCCCGAGCTGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCGTCCTCGGCCACGGCGGGAGAAAATTCtgc                                                                                                                                                   >  6:100371/4‑149 (MQ=255)
gcctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgc                                                                                                                                                   >  4:268127/4‑149 (MQ=255)
 cctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgc                                                                                                                                                   <  7:56927/146‑1 (MQ=255)
 cctCCACCGCGGCCACGAAAGCGGCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgcc                                                                                                                                                  <  4:68585/147‑2 (MQ=255)
   tCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccag                                                                                                                                                >  2:211174/1‑146 (MQ=255)
   tCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccag                                                                                                                                                >  1:82657/1‑146 (MQ=255)
   tCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTCGCCTCAAAAGAGGCAGCAAAATCATCAGAACCGAAAGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccag                                                                                                                                                >  2:126822/1‑146 (MQ=255)
    ccACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGCGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCTGCAAATTCtgccagg                                                                                                                                               >  8:361791/1‑145 (MQ=255)
    ccACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGCCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccagg                                                                                                                                               >  8:338661/1‑145 (MQ=255)
     cACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccaggg                                                                                                                                              >  1:100674/1‑144 (MQ=255)
     cACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccaggg                                                                                                                                              <  8:242611/149‑6 (MQ=255)
     cACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccaggg                                                                                                                                              <  6:325331/149‑6 (MQ=255)
       ccGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGggcg                                                                                                                                            <  7:2317/149‑1 (MQ=255)
        cgcgGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGggcgg                                                                                                                                           <  1:262150/149‑1 (MQ=255)
         gcgGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCCCGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGggcggc                                                                                                                                          <  4:330632/149‑1 (MQ=255)
             ccACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCaaa                                                                                                                                       >  3:261460/1‑148 (MQ=255)
               aCGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAAc                                                                                                                                    <  6:85990/149‑1 (MQ=255)
                  aaagcgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGATGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTc                                                                                                                                 <  2:155567/143‑1 (MQ=255)
                  aaagcgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGt                                                                                                                                  >  8:173470/7‑148 (MQ=255)
                  aaagcgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTc                                                                                                                                 <  6:272150/143‑1 (MQ=255)
                  aaagcgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTc                                                                                                                                 <  5:135644/143‑1 (MQ=255)
                   aagcgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTcc                                                                                                                                <  4:148527/144‑1 (MQ=255)
                   aagcgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTcc                                                                                                                                >  1:171818/6‑149 (MQ=255)
                    agcgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTcc                                                                                                                                >  3:206896/5‑148 (MQ=255)
                       gtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCgaga                                                                                                                            <  4:219441/148‑1 (MQ=255)
                       aTCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGACCGCATCATCAAGAGCCAGGAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAATCTTCCgaga                                                                                                                            >  1:395414/1‑149 (MQ=255)
                         cAGAGGCCGCCACTTCAGCACGAGATGCGGCGGCCTCAAAATAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGCTCGTGCAGCTTCCTCGGCAACGGCGGCATAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACg                                                                                                                          <  6:52717/149‑1 (MQ=255)
                          agagGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTc                                                                                                                                 <  5:269046/141‑1 (MQ=255)
                            aaGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGa                                                                                                                                                 <  6:243548/123‑1 (MQ=255)
                                       tCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGtcttc                                                                                                            >  7:89886/1‑149 (MQ=255)
                                         agcacgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGACCGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGgcatct                                                                                                           <  1:363185/141‑1 (MQ=255)
                                          gcacgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGa                                                                                                         <  2:188442/143‑1 (MQ=255)
                                           cacgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTcc                                                                                                                                <  4:125068/120‑1 (MQ=255)
                                            ccGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGaaa                                                                                                       <  4:302336/149‑1 (MQ=255)
                                            acgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGaaa                                                                                                       <  1:140318/145‑1 (MQ=255)
                                             ccagATGCGGGGGCCTCAAAAGAGGCAGCAACATCCTCAGAACCGAACGCATCATCAAGTGCCGGTCGCGCAGCTTCCTCCGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGGCATCTGaacc                                                                                                      <  8:314426/146‑3 (MQ=255)
                                              gagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGAc                                                                                                                           >  1:299334/3‑127 (MQ=255)
                                               agATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAcag                                                                                                    <  2:274559/148‑3 (MQ=255)
                                                gATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCGACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGc                                                                                                   >  8:114862/1‑149 (MQ=255)
                                                 aTGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCa                                                                                                  <  8:209202/149‑1 (MQ=255)
                                                  tGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAcagcag                                                                                                 <  6:179950/149‑1 (MQ=255)
                                                  tGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAcagcag                                                                                                 <  2:8591/149‑1 (MQ=255)
                                                   gcggtgGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACagcagc                                                                                                <  3:30739/149‑1 (MQ=255)
                                                    cggtggCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACagcagc                                                                                                <  2:302634/148‑1 (MQ=255)
                                                       tggcCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCggaa                                                                                            <  7:109617/148‑3 (MQ=255)
                                                       tggcCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCggaa                                                                                            >  7:247947/2‑147 (MQ=255)
                                                         gcCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAAcg                                                                                          <  3:48791/149‑2 (MQ=255)
                                                           cTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgc                                                                                                                                                   <  6:214929/90‑1 (MQ=255)
                                                           cTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAAcg                                                                                          >  1:347892/1‑146 (MQ=255)
                                                           cTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAAcg                                                                                          >  4:236512/1‑146 (MQ=255)
                                                            tCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGcc                                                                                                                    >  6:157393/1‑120 (MQ=255)
                                                            tCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAcagcag                                                                                                 >  7:395090/1‑139 (MQ=255)
                                                             cAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAgc                                                                                      >  2:352159/1‑149 (MQ=255)
                                                              aaaaGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAATGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAgc                                                                                      <  3:174103/148‑1 (MQ=255)
                                                              aaaaGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAgcg                                                                                     >  6:293070/1‑149 (MQ=255)
                                                              aaaaGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAgcg                                                                                     >  7:281868/1‑149 (MQ=255)
                                                               aaaGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGa                                                                                                                         >  2:327888/1‑112 (MQ=255)
                                                                aagaggCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACggag                                                                                       <  3:110368/145‑1 (MQ=255)
                                                                aagaggCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGcc                                                                                   >  1:318255/1‑149 (MQ=255)
                                                                 agaggCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCt                                                                                  <  2:50274/149‑1 (MQ=255)
                                                                   aggcagCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTct                                                                                >  3:233351/6‑148 (MQ=255)
                                                                     gcagCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAca                                                                                                     >  7:136405/4‑125 (MQ=255)
                                                                     gcagCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGc                                                                              <  5:221662/146‑1 (MQ=255)
                                                                     gcagCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGc                                                                              >  8:267714/4‑149 (MQ=255)
                                                                      cGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCt                                                                             <  3:141314/149‑1 (MQ=255)
                                                                       agCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGccg                                                                            >  7:309892/2‑147 (MQ=255)
                                                                        gCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGc                                                                           <  1:232999/149‑1 (MQ=255)
                                                                          aaaaTCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCgg                                                                         >  3:89638/1‑149 (MQ=255)
                                                                          aaaaTCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCgg                                                                         <  2:378653/149‑1 (MQ=255)
                                                                          aaaaTCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCgg                                                                         <  6:353397/149‑1 (MQ=255)
                                                                           aaaTCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTGCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGc                                                                        <  7:169300/149‑1 (MQ=255)
                                                                            aaTCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCa                                                                       >  5:55454/1‑149 (MQ=255)
                                                                              tcatcaGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAga                                                                     >  6:376672/1‑148 (MQ=255)
                                                                                atcaGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAgacg                                                                   >  1:58150/1‑146 (MQ=255)
                                                                                 tcaGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAgacg                                                                   >  5:159833/1‑145 (MQ=255)
                                                                                    gAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAgg                                                                                                                 >  5:106395/1‑99 (MQ=255)
                                                                                    gAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCaaa                                                               <  1:278646/149‑1 (MQ=255)
                                                                                    gAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCaa                                                                >  4:173759/1‑148 (MQ=255)
                                                                                       acgaacgCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAAc                                                            >  5:762/1‑149 (MQ=255)
                                                                                        cgaacgCCTCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACa                                                           <  8:141753/149‑1 (MQ=255)
                                                                                        cgaacgCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACa                                                           <  7:369837/149‑1 (MQ=255)
                                                                                         gaacgCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAg                                                          >  5:244465/1‑149 (MQ=255)
                                                                                           acgCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAAcg                                                                                          >  5:375364/1‑114 (MQ=255)
                                                                                            cgCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAGCgg                                                       <  5:318001/149‑1 (MQ=255)
                                                                                             gCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAGCg                                                        <  8:62739/147‑1 (MQ=255)
                                                                                             gCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAGCGGc                                                      >  3:324271/1‑149 (MQ=255)
                                                                                              catcatCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCgg                                                                         >  5:195793/1‑129 (MQ=255)
                                                                                               atcatcAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAGCGGcgg                                                    >  5:348648/1‑147 (MQ=255)
                                                                                               atcatcAAGTGCCGGTCGTGCACCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACCTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGACCCGAGCGCCTCTCCCGCGGCAGACGCAAAAACAGCGGccg                                                    >  5:295511/1‑147 (MQ=255)
                                                                                                tcatcaAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGccg                                                                            >  6:122274/1‑122 (MQ=255)
                                                                                                   tcaAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAgcggcggcggg                                                <  2:10484/149‑3 (MQ=255)
                                                                                                   tcaAGTGCCGGTCGGGCAGCTTCCTCGGCAACGGCGGCAGAAATTTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGATCGTAGCGCCTCTGCCGCGGCAGACGCAAAAACAgcggcggcggg                                                <  1:121562/149‑3 (MQ=255)
                                                                                                    caAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAgcggcggcggg                                                >  4:358310/1‑146 (MQ=255)
                                                                                                     aaGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAgcggcggcgggga                                              <  7:198563/149‑5 (MQ=255)
                                                                                                     aaGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAgcggcggcggg                                                >  2:246284/1‑145 (MQ=255)
                                                                                                      agtgcCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAg                                                                                                                  >  1:196169/4‑80 (MQ=255)
                                                                                                                      cTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAgacg                                                                   <  8:392267/111‑4 (MQ=255)
                                                                                                                           tCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAgcggcggcg                                                  >  3:346208/1‑123 (MQ=255)
                                                                                                                                    gcggcAGGAAATTCCGCGAAGGCGGCAAAAACGTGCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGTCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCa               <  3:282024/149‑1 (MQ=255)
                                                                                                                                             aaTTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAgcggcggcggg                                                <  7:337943/107‑3 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGT  >  NC_000913/1429712‑1430005

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