Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,072,491 C→T 5.4% T318M (ACG→ATG)  flu → CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,072,4910CT 5.4% 233.6 / 2.6 74T318M (ACG→ATG) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base C (32/38);  new base T (2/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

TTGTGAACAGTGACGGCTGGCAGATTGTCAAAAACGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAA  >  NC_000913/2072344‑2072637
                                                                                                                                                   |                                                                                                                                                  
ttgtGAACAGTGACGGTTGGCAGATTGTCAAAAACGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACg                                                                                                                                                   >  1:130315/1‑149 (MQ=255)
ttgtGAACAGTGACGGCTGGCAGATTGTCAAAAACGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTAACCCTGAAGCAGGGCGGCGCTCTGGTTACCAGTACg                                                                                                                                                   >  5:302542/1‑149 (MQ=255)
    gAACAGTGACGGCTGGCAGATTGTCAAAAACGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTg                                                                                                                                               >  7:213969/1‑149 (MQ=255)
    gAACAGTGACGGCTGGCAGATTGTCAAAAACGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACACTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTg                                                                                                                                               >  7:350367/1‑149 (MQ=255)
          tGACGGCTGGCAGATTGTCAAAAACGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCg                                                                                                                                         >  2:49287/1‑149 (MQ=255)
           gACGGCTGGCAGATTGTCAAAAACGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGAACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTCGTTACCAGTACGGCTGCAACCGt                                                                                                                                        >  5:155158/1‑149 (MQ=255)
                   gCAGATTGTCAAAAACGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCTGTACGGCTGCAACCGTTACCGGc                                                                                                                                 <  1:141695/148‑1 (MQ=255)
                   gCAGATTGTCAAAAACGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGc                                                                                                                                 <  1:353649/148‑1 (MQ=255)
                      gATTGTCAAAAACGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATa                                                                                                                              >  4:276091/1‑148 (MQ=255)
                       aTTGTCAAAAACGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATaaa                                                                                                                            <  2:130315/149‑1 (MQ=255)
                           tCAAAAACGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGAGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGc                                                                                                                        >  5:41887/1‑149 (MQ=255)
                             aaaaaCGGGGGTGTGGCCGGGAATACCTCCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTAGCGGCATAAACCGCCt                                                                                                                      <  8:84732/149‑1 (MQ=255)
                             aaaaaCGGGGGTGTGGCCGGGAATACCACCGTTAATCAGCAGGGCAGACTGCAGGTGGACGCCGGGGGTACAGCCACGAATGGCACCCTGAAGCAGGGCGGCTCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCt                                                                                                                      <  5:319654/149‑1 (MQ=255)
                              aaaaCGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCTGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCta                                                                                                                     <  6:157336/149‑2 (MQ=255)
                              aaaaCGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTg                                                                                                                     <  7:17240/149‑1 (MQ=255)
                               aaaCGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCGCTGGTTACCAGTACGGCGGCAACCGTTACCGGCATAAACCGCCTgg                                                                                                                    >  2:228186/1‑149 (MQ=255)
                                aaCGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTggg                                                                                                                   <  4:178062/149‑1 (MQ=255)
                                 aCGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGa                                                                                                                  <  5:275898/149‑1 (MQ=255)
                                        tgtgGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATtct                                                                                                           >  7:250538/1‑149 (MQ=255)
                                        tgtgGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATtct                                                                                                           >  1:229574/1‑149 (MQ=255)
                                         ggggCCGGGAATACCACCGTTAATCATAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGGCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGAATAAACCGCCTGGGAGCATtctc                                                                                                          <  6:28690/147‑1 (MQ=255)
                                               gggAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACATCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCtgttgt                                                                                                    >  7:44182/1‑149 (MQ=255)
                                                ggAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTtgtg                                                                                                   <  1:138438/149‑1 (MQ=255)
                                                 gAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGCTACAGACACGAATGTCACCCTGAAGCAGGGCGGCGCACAGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTtgtg                                                                                                   >  7:304875/1‑148 (MQ=255)
                                                  aattacactGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGa                                                                                                 <  8:258980/144‑1 (MQ=255)
                                                    tACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTATGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAgg                                                                                               >  7:103141/1‑149 (MQ=255)
                                                    tACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTATGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAgg                                                                                               >  7:103143/1‑149 (MQ=255)
                                                      ccaccGTTAATCAGAAGGGCAGACTGCAGGTGAACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGt                                                                                             <  6:155158/149‑1 (MQ=255)
                                                       caccGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCtgtt                                                                                                      >  2:84258/1‑139 (MQ=255)
                                                       caccGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCtgtt                                                                                                      <  1:84258/139‑1 (MQ=255)
                                                       caccGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGt                                                                                             >  5:327042/1‑148 (MQ=255)
                                                               aTCAGATGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTTGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGa                                                                                     <  2:10119/148‑1 (MQ=255)
                                                               aTCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGAt                                                                                    <  3:102868/149‑1 (MQ=255)
                                                                tCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACGCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATTAACCGACTGGGAGCATTCTATGTTGTGGAGGGTAAAGCTGAt                                                                                    >  7:124041/1‑148 (MQ=255)
                                                                  aGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAAt                                                                                 <  2:92169/149‑1 (MQ=255)
                                                                   gAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATg                                                                                <  1:49287/149‑1 (MQ=255)
                                                                   gAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCAGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGTCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATg                                                                                <  2:122479/149‑1 (MQ=255)
                                                                    aaGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGt                                                                               <  3:276091/149‑1 (MQ=255)
                                                                      gggCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCAAACTCAAGAAGGGCGGCGCAATGGTTACCAGTACTGCTTCAACCGTTACCGGCATAAAGCGCCTGGGAGCATTATCTGTTTGGGAGGGTAAAGCTGATAATGTCg                                                                             >  2:96424/1‑149 (MQ=255)
                                                                        gCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTa                                                                           <  6:334967/149‑1 (MQ=255)
                                                                         cAGACTGCAGTTGTACCCCGGAGTTACAGCCACGAATGTCACCCTGAAGCAGGCCGGCGCAATGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTAc                                                                          >  2:299293/1‑149 (MQ=255)
                                                                          aGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACt                                                                         <  2:229574/149‑1 (MQ=255)
                                                                            aCTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTgg                                                                       <  6:49968/149‑1 (MQ=255)
                                                                                cAGGTGGACGCCTGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGaaaa                                                                   <  6:41887/149‑1 (MQ=255)
                                                                                  ggTGGACGCCGGGGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATg                                                                 <  8:289392/149‑1 (MQ=255)
                                                                                    tGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGc                                                               <  8:213969/149‑1 (MQ=255)
                                                                                      gACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCg                                                              >  1:254749/1‑148 (MQ=255)
                                                                                       aCGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTATGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGa                                                            <  8:103141/149‑1 (MQ=255)
                                                                                       aCGCCGGTGGTACAGCCACGAATGGCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTATGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGa                                                            <  8:103143/149‑1 (MQ=255)
                                                                                                  aCAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGAtgt                                                 >  6:338241/1‑149 (MQ=255)
                                                                                                    aGCCACGCATGTCCCCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGGAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGc                                               >  4:13981/1‑149 (MQ=255)
                                                                                                    aGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGc                                               <  6:327042/149‑1 (MQ=255)
                                                                                                        aCGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGAc                                           >  7:255853/1‑149 (MQ=255)
                                                                                                        aCGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGAc                                           >  2:134992/1‑149 (MQ=255)
                                                                                                          gAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGACGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCg                                         >  4:236618/1‑149 (MQ=255)
                                                                                                            aTGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTTCCGGAATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGa                                       <  8:350367/149‑1 (MQ=255)
                                                                                                                cACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGacac                                    <  8:304875/148‑1 (MQ=255)
                                                                                                                cACCCTGAAGCAGGGCGGCGCACTGGGTACCAGTACGGCGGCAACCGTTACCGGCATAAACCGCCTGCGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGTAAAATGGCTGACGCCTGGATGTGCTGACCGGacaca                                   <  1:96424/149‑1 (MQ=255)
                                                                                                                  cccTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGTGGACGCCTGGATGTGCTGACCGGacaca                                   >  7:129895/1‑147 (MQ=255)
                                                                                                                  cccTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGacaca                                   >  3:295894/1‑147 (MQ=255)
                                                                                                                  cccTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGacaca                                   >  4:384122/1‑147 (MQ=255)
                                                                                                                   ccTGAAGCGGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAg                                >  1:384514/1‑149 (MQ=255)
                                                                                                                      gAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCa                             >  2:342851/1‑149 (MQ=255)
                                                                                                                       aaGCAGGGAGGCGCACTGGTTACCAGTACGGATGCACCCGTTACCTGCATAAACCTCCGGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATTGCGGACGCCTGGATGTGCTGACCGGACACACAGCCa                             <  6:302542/148‑1 (MQ=255)
                                                                                                                          cAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTaa                         <  8:250538/149‑1 (MQ=255)
                                                                                                                              gcggcgCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATAcc                     <  2:254749/149‑1 (MQ=255)
                                                                                                                                 gcgcACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCcgc                  <  4:295894/149‑1 (MQ=255)
                                                                                                                                          ttACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGgatgat         <  1:134992/149‑1 (MQ=255)
                                                                                                                                             ccAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGAc                                           <  3:163290/112‑1 (MQ=255)
                                                                                                                                             ccAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGAc                                           >  4:163290/1‑112 (MQ=255)
                                                                                                                                             ccAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGc      <  1:228186/149‑1 (MQ=255)
                                                                                                                                                gTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCgg    >  6:205917/1‑148 (MQ=255)
                                                                                                                                                gTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCgg    >  7:255199/1‑148 (MQ=255)
                                                                                                                                                 tACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGaa  <  1:299293/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGTGAACAGTGACGGCTGGCAGATTGTCAAAAACGGGGGTGTGGCCGGGAATACCACCGTTAATCAGAAGGGCAGACTGCAGGTGGACGCCGGTGGTACAGCCACGAATGTCACCCTGAAGCAGGGCGGCGCACTGGTTACCAGTACGGCTGCAACCGTTACCGGCATAAACCGCCTGGGAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAA  >  NC_000913/2072344‑2072637

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