Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,073,397 C→T 5.6% A620V (GCA→GTA)  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,073,3970CT 5.6% 229.1 / 2.5 71A620V (GCA→GTA) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base C (36/30);  new base T (2/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

AGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGC  >  NC_000913/2073253‑2073541
                                                                                                                                                |                                                                                                                                                
aGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTCCGCCCCGATATGGca                                                                                                                                                 >  5:264335/1‑146 (MQ=255)
   cGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGaa                                                                                                                                            >  2:369744/1‑148 (MQ=255)
        gtggATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATg                                                                                                                                       <  3:90447/148‑1 (MQ=255)
         tggATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCt                                                                                                                                     <  8:301635/149‑1 (MQ=255)
            aTGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGAc                                                                                                                                  <  1:369744/149‑1 (MQ=255)
                ccTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGAc                                                                                                                              >  1:10649/1‑149 (MQ=255)
                 cTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACt                                                                                                                             <  8:78329/149‑1 (MQ=255)
                  tCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTAGTACCGGAAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGCCAACATCAGCCTACGAGTACGCCCCGATAAGGAAAATAAAAATGATGACAGACTg                                                                                                                            >  7:276300/1‑149 (MQ=255)
                    aGCCATGCCGGAAAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGACCGGATATGGCACGTAACAATGCTGACAGACAGGt                                                                                                                          >  6:906/1‑149 (MQ=255)
                      ccATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGTCGAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGCCTGGACa                                                                                                                        >  4:43421/1‑149 (MQ=255)
                      ccATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCa                                                                                                                        >  1:115551/1‑149 (MQ=255)
                        aTGCCGGACAGATTCATTTCACCTCCACCCACACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCAtt                                                                                                                      <  2:218838/149‑1 (MQ=255)
                                   aTTCATTTCACCTCCACCCGCAGAGGGAAGTTCGTACCGGCAACGCTGAAAGTGAAAAACCTGAACGGGCAGATTGGCTCACTCAGCCTGCGTGTACGCACGGATATGGCACAGCCCAATGCTGACAGACTGCTCATGGACGGCGGCAg                                                                                                           >  8:210740/1‑149 (MQ=255)
                                   aTTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACCCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAg                                                                                                           >  1:193821/1‑149 (MQ=255)
                                       aTTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGc                                                                                                        >  7:322819/1‑148 (MQ=255)
                                          tCACCTCCACACGCACAGGGAAGTTCGTACCGGCAACCCTGAAATTGAAAAACCTGAACGGAAAGAATGGCCCCATCAGCCGGCGTGGACGCCCGGATATGGCACAGCACAATGCTGCCAGCCTGGGCCTTGCCGGCGGCAGGGCCAcc                                                                                                    >  6:34427/1‑149 (MQ=255)
                                            aCCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCg                                                                                                   >  1:354418/1‑148 (MQ=255)
                                                ccccccGCACCGGGAAGTTCGTACCGGCAACCCTCAAAGGGAAAACCCTGAACGGACCGACAGGCACCATCAGCCTGCGTGGACGCCCGGATATGGCACAGCACACTGCTGACACACTGCGCATTGACGGCGGCAGGGCAACCGGaaaa                                                                                              <  6:137004/146‑1 (MQ=255)
                                                 cACCCGCACAGGGAAGTGCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGaaaaa                                                                                             <  3:344569/149‑1 (MQ=255)
                                                            ggAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAAc                                                                                  <  2:115551/149‑1 (MQ=255)
                                                               aGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTg                                                                               <  2:354418/149‑1 (MQ=255)
                                                                gTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTgg                                                                              >  8:221876/1‑149 (MQ=255)
                                                                   cGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTg                                                                            <  4:360475/148‑1 (MQ=255)
                                                                        cGGCAACCCTGAAAGTGAAAAACCTCCACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGACATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGcc                                                                      >  3:162000/1‑149 (MQ=255)
                                                                         ggCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCg                                                                     >  8:171494/1‑149 (MQ=255)
                                                                         ggCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCg                                                                     >  8:49116/1‑149 (MQ=255)
                                                                          gCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGTATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGTTGAACGCCgg                                                                    <  2:10649/149‑1 (MQ=255)
                                                                          gCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCgg                                                                    >  5:26296/1‑149 (MQ=255)
                                                                          gCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCgg                                                                    >  3:248875/1‑149 (MQ=255)
                                                                              cccTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCaa                                                                 >  6:239609/1‑148 (MQ=255)
                                                                               ccTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAAGAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGATAAACCATCCTGAACCTGGTGAACGCCGGCAACa                                                               <  5:906/149‑1 (MQ=255)
                                                                               ccTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCaa                                                                 >  5:220119/1‑147 (MQ=255)
                                                                               ccTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCaa                                                                 >  1:255259/1‑147 (MQ=255)
                                                                                cTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAggg                                                                                                         >  1:240604/1‑106 (MQ=255)
                                                                                cTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAggg                                                                                                         <  2:240604/106‑1 (MQ=255)
                                                                                 tGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGt                                                             >  4:191386/1‑149 (MQ=255)
                                                                                 tGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGt                                                             >  3:264944/1‑149 (MQ=255)
                                                                                   aaaGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGt                                                             <  7:171494/147‑1 (MQ=255)
                                                                                   aaaGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGt                                                             <  7:49116/147‑1 (MQ=255)
                                                                                       tGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCg                                                       <  4:248875/149‑1 (MQ=255)
                                                                                         aaaaaCCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGCAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCggg                                                     >  5:339122/1‑149 (MQ=255)
                                                                                          aaaaCCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCgggg                                                    >  3:200605/1‑149 (MQ=255)
                                                                                          aaaaCCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCgggg                                                    <  8:322819/149‑1 (MQ=255)
                                                                                           aaaCCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGc                                                   >  1:145926/1‑149 (MQ=255)
                                                                                            aaCCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCt                                                  >  3:141474/1‑149 (MQ=255)
                                                                                              ccTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTgg                                                >  5:238093/1‑149 (MQ=255)
                                                                                                  aaCGGACAGAATGGCACCATCAGCCTGCGGGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCg                                              <  3:43421/147‑1 (MQ=255)
                                                                                                  aaCGGACAGAAGGGCAGCCTCCGCCTGCGTGTACGCCCCGATATGGCACTGCACAATTCTGACCGACTGGTCAGTGACGGCGGCCGGGCTACCGGAAACACCATCCTCAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGa                                             <  4:162000/148‑1 (MQ=255)
                                                                                                       aCAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGGCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCg                                       <  5:34427/149‑1 (MQ=255)
                                                                                                           aaTGGCACCATCAGCCTGCGTGTACGCCCGGATATGGTACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGCACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGt                                     >  4:335278/1‑147 (MQ=255)
                                                                                                           aaTGGCACCATCAGCCTGCGTGTACGCCCGGATATGGTACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGt                                     >  4:329972/1‑147 (MQ=255)
                                                                                                           aaTGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGACCCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGt                                     >  4:211939/1‑147 (MQ=255)
                                                                                                           aaTGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTaa                                   <  6:220119/149‑1 (MQ=255)
                                                                                                             ttgCACCATCAGCCTGCGTGTACGCCCGGATCTGGCACAGAACAATGCTGACAGACTGGTAATTGACGGCGGCAGGGCAACCTGAAAAACCATCCTGAACCTTGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAgg                                 <  6:264335/147‑1 (MQ=255)
                                                                                                             tGGCACCATCAGCCTGCGTGTACGCCCGGATATGGTACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAgg                                 <  3:329972/149‑1 (MQ=255)
                                                                                                             tGGCACCATCAGCCTGCGTGTACGCCCGGATATGGTACAGAACAATGCTGACAGACTGGGCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAgg                                 <  3:335278/149‑1 (MQ=255)
                                                                                                             tGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAgg                                 <  5:239609/149‑1 (MQ=255)
                                                                                                                     tCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAg                         >  4:163218/1‑149 (MQ=255)
                                                                                                                        gCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAggtg                      >  3:338075/1‑149 (MQ=255)
                                                                                                                          cTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAggtggt                    <  3:163218/149‑1 (MQ=255)
                                                                                                                              gtgtACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGaa                <  2:145926/149‑1 (MQ=255)
                                                                                                                                 tACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGCACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGc              >  7:297812/1‑148 (MQ=255)
                                                                                                                                    gCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCAt           >  6:268674/1‑148 (MQ=255)
                                                                                                                                    gCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCAt           >  4:158234/1‑148 (MQ=255)
                                                                                                                                    gCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCAt           >  6:195455/1‑148 (MQ=255)
                                                                                                                                      ccGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTaa        >  8:44273/1‑149 (MQ=255)
                                                                                                                                        ggATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACg      <  6:26296/149‑1 (MQ=255)
                                                                                                                                          atatGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGt    <  7:221876/149‑1 (MQ=255)
                                                                                                                                          atatGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGt    <  2:255259/149‑1 (MQ=255)
                                                                                                                                           tatGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGcc             >  7:242742/1‑139 (MQ=255)
                                                                                                                                           tatGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGcc             <  8:242742/139‑1 (MQ=255)
                                                                                                                                             tGGCACAGAACAATGCAGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGATCGTGGTGTACGCCGGCAACAGTGCGTCGGGGCTGGCGACGAGCGGTAAGGGTATTCAGGTGGTGGTAGCCATTAACGGTGc  <  2:193821/148‑1 (MQ=255)
                                                                                                                                                cacaGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAAc                                                                                                     <  2:373808/46‑1 (MQ=255)
                                                                                                                                                cacaGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAAc                                                                                                     >  1:373808/1‑46 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGAATGGCACCATCAGCCTGCGTGTACGCCCGGATATGGCACAGAACAATGCTGACAGACTGGTCATTGACGGCGGCAGGGCAACCGGAAAAACCATCCTGAACCTGGTGAACGCCGGCAACAGTGCGTCGGGGCTGGCGACCAGCGGTAAGGGTATTCAGGTGGTGGAAGCCATTAACGGTGC  >  NC_000913/2073253‑2073541

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: