Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,073,157 T→C L540P (CTG→CCG)  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,1570TC99.3% 496.8 / NA 136L540P (CTG→CCG) fluCP4‑44 prophage, antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (63/72);  total (63/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

GAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTC  >  NC_000913/2073025‑2073291
                                                                                                                                    |                                                                                                                                      
     ggCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGgcagca                                                                                                                             >  1:610685/1‑139 (MQ=255)
     ggCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGgcagca                                                                                                                             <  1:558373/139‑1 (MQ=255)
     ggCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGgcagca                                                                                                                             <  2:496068/139‑1 (MQ=255)
      gCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCAc                                                                                                                            >  1:1115652/1‑139 (MQ=255)
      gCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCAc                                                                                                                            >  1:264208/1‑139 (MQ=255)
        acacTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACtg                                                                                                                          <  2:114282/139‑1 (MQ=255)
          actcactGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACtgtg                                                                                                                        >  1:90710/1‑139 (MQ=255)
               ctGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGaa                                                                                                                   <  2:213500/139‑1 (MQ=255)
                      cAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGgcagca                                                                                                                             >  1:63556/1‑122 (MQ=255)
                      cAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGgcagca                                                                                                                             <  2:63556/122‑1 (MQ=255)
                      cAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGcc                                                                                                            >  2:1276902/1‑139 (MQ=255)
                      cAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGcc                                                                                                            >  1:317697/1‑139 (MQ=255)
                        acacCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCAt                                                                                                          >  1:267499/1‑139 (MQ=255)
                        acacCACACTCACCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCAtt                                                                                                         >  2:266831/1‑139 (MQ=255)
                        acacCACACTCACCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCAtt                                                                                                         <  1:266831/139‑1 (MQ=255)
                         caccacACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCAtt                                                                                                         <  1:1777266/139‑1 (MQ=255)
                            cacaCTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGAc                                                                                                      >  1:1017337/1‑139 (MQ=255)
                                cTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCa                                                                                                  >  1:1445504/1‑139 (MQ=255)
                                 tCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCAc                                                                                                 <  2:610685/139‑1 (MQ=255)
                                   aCCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGa                                                                                               >  1:1420746/1‑139 (MQ=255)
                                    cccAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGaa                                                                                              >  2:1409257/1‑139 (MQ=255)
                                     ccAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAAt                                                                                             >  2:199340/1‑139 (MQ=255)
                                      cAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATg                                                                                            >  1:1462278/1‑139 (MQ=255)
                                             ccGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGgcagca                                                                                                                             <  2:1469036/99‑1 (MQ=255)
                                             ccGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGgcagca                                                                                                                             >  1:1469036/1‑99 (MQ=255)
                                             ccGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCActct                                                                                     >  1:1783158/1‑139 (MQ=255)
                                             ccGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCActct                                                                                     >  1:153250/1‑139 (MQ=255)
                                               gTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCg                                                                                   >  1:1386680/1‑139 (MQ=255)
                                               gTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCg                                                                                   >  1:1493667/1‑139 (MQ=255)
                                                   aCCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGGCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTc                                                                               <  2:156950/139‑1 (MQ=255)
                                                        aatgaaGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTg                                                                          <  1:1715205/139‑1 (MQ=255)
                                                           gaaGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCa                                                                       <  2:1445504/139‑1 (MQ=255)
                                                            aaGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCAc                                                                      <  1:1409257/139‑1 (MQ=255)
                                                             aGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCAcc                                                                     <  1:86197/139‑1 (MQ=255)
                                                              ggCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCt                                                                    >  1:415220/1‑139 (MQ=255)
                                                              ggCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCt                                                                    >  2:268608/1‑139 (MQ=255)
                                                              ggCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCt                                                                    >  2:1766796/1‑139 (MQ=255)
                                                              ggCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCt                                                                    >  1:775602/1‑139 (MQ=255)
                                                                cACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTgg                                                                  >  2:1649809/1‑139 (MQ=255)
                                                                 acgctgacgctgAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGa                                                                 <  2:467450/139‑1 (MQ=255)
                                                                  cgctgacgctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTg                                                                                                                     >  2:349885/1‑86 (MQ=255)
                                                                  cgctgacgctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAATGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGt                                                                           <  1:1441331/128‑1 (MQ=255)
                                                                         gctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATcccc                                                         >  1:1562449/1‑139 (MQ=255)
                                                                           tgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGa                                                       <  2:841135/139‑1 (MQ=255)
                                                                             aaCGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATa                                                     >  1:1209075/1‑139 (MQ=255)
                                                                              aCGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATaa                                                    <  2:1562449/139‑1 (MQ=255)
                                                                                 aCAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGc                                                 <  2:353924/139‑1 (MQ=255)
                                                                                  cAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGcc                                                <  2:1783158/139‑1 (MQ=255)
                                                                                   aGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCa                                               <  2:260556/139‑1 (MQ=255)
                                                                                    gTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCAc                                              <  2:1462278/139‑1 (MQ=255)
                                                                                     tACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCgg                                                                                                                                  <  1:69799/54‑1 (MQ=255)
                                                                                     tACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCgg                                                                                                                                  >  2:69799/1‑54 (MQ=255)
                                                                                     tACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACg                                             <  2:1420746/139‑1 (MQ=255)
                                                                                        cGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTg                                          <  2:1479577/139‑1 (MQ=255)
                                                                                        cGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTg                                          <  1:1125874/139‑1 (MQ=255)
                                                                                          tCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTTTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTg                                          <  1:1718156/137‑1 (MQ=255)
                                                                                          tCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCa                                        >  2:349125/1‑139 (MQ=255)
                                                                                            accacGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGt                                      <  1:1133003/139‑1 (MQ=255)
                                                                                            accacGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGt                                      <  1:218476/139‑1 (MQ=255)
                                                                                               acGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCgg                                   <  2:90710/139‑1 (MQ=255)
                                                                                                cGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTg                                                                          >  1:770307/1‑99 (MQ=255)
                                                                                                cGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTg                                                                          <  2:770307/99‑1 (MQ=255)
                                                                                                 ggATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTg                                                                          <  2:749435/98‑1 (MQ=255)
                                                                                                 ggATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTg                                                                          >  1:749435/1‑98 (MQ=255)
                                                                                                   aTGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCggtggt                               >  1:88968/1‑139 (MQ=255)
                                                                                                    tGTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGgtggtg                              >  2:287718/1‑139 (MQ=255)
                                                                                                     gTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGtggtgg                             >  1:1597730/1‑139 (MQ=255)
                                                                                                     gTCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGtggtgg                             <  2:1493667/139‑1 (MQ=255)
                                                                                                      tCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                            <  1:335081/139‑1 (MQ=255)
                                                                                                      tCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                            <  1:138547/139‑1 (MQ=255)
                                                                                                      tCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                            <  1:917271/139‑1 (MQ=255)
                                                                                                      tCATTGCTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                            <  2:1017337/139‑1 (MQ=255)
                                                                                                            cTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCa                                                                                                           >  2:1606577/1‑54 (MQ=255)
                                                                                                            cTCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCa                                                                                                           <  1:1606577/54‑1 (MQ=255)
                                                                                                             tCAGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCggtg                                 >  1:342500/1‑127 (MQ=255)
                                                                                                               aGCGCGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTg                                                                          <  1:180516/84‑1 (MQ=255)
                                                                                                                gcgcGGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGc                  <  1:199340/139‑1 (MQ=255)
                                                                                                                 cgcgGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGt                                                                           >  2:90154/1‑81 (MQ=255)
                                                                                                                 cgcgGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGcc                 >  1:650273/1‑139 (MQ=255)
                                                                                                                  gcgGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCt                                                                                <  1:567764/75‑1 (MQ=255)
                                                                                                                  gcgGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCt                                                                                >  2:567764/1‑75 (MQ=255)
                                                                                                                   cgGTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCAt               <  2:775602/139‑1 (MQ=255)
                                                                                                                    ggTACAGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATg              <  2:264208/139‑1 (MQ=255)
                                                                                                                         aGCCCTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGa         <  2:1584603/139‑1 (MQ=255)
                                                                                                                         aGCCCCGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCAGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGa         <  2:1115652/139‑1 (MQ=255)
                                                                                                                           cccTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGAc                                                                                                      >  2:1638254/1‑44 (MQ=255)
                                                                                                                           cccTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGAc                                                                                                      <  1:1638254/44‑1 (MQ=255)
                                                                                                                           cccTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGt                                                                           >  1:594818/1‑71 (MQ=255)
                                                                                                                           cccTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGt                                                                           <  2:594818/71‑1 (MQ=255)
                                                                                                                            ccTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAg      <  2:1379725/139‑1 (MQ=255)
                                                                                                                             cTGAAGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGa     <  2:712793/139‑1 (MQ=255)
                                                                                                                                aaGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCAct                                                                                       >  2:114103/1‑54 (MQ=255)
                                                                                                                                aaGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCAct                                                                                       <  1:114103/54‑1 (MQ=255)
                                                                                                                                aaGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCtt                                                                               >  2:833078/1‑61 (MQ=255)
                                                                                                                                aaGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGt                                                                           >  1:176461/1‑66 (MQ=255)
                                                                                                                                aaGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGt                                                                           >  2:838464/1‑66 (MQ=255)
                                                                                                                                aaGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGt                                                                           <  1:838464/66‑1 (MQ=255)
                                                                                                                                aaGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGt                                                                           <  2:176461/66‑1 (MQ=255)
                                                                                                                                aaGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTc  >  2:327166/1‑139 (MQ=255)
                                                                                                                                aaGCCGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTc  >  1:997722/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                      
GAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTC  >  NC_000913/2073025‑2073291

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