Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,072,659 G→A G374D (GGC→GAC)  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,6590GA100.0% 303.6 / NA 86G374D (GGC→GAC) fluCP4‑44 prophage, antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (42/44);  total (42/44)

GAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGG  >  NC_000913/2072524‑2072790
                                                                                                                                       |                                                                                                                                   
gAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGacac                                                                                                                                  <  2:3000/139‑1 (MQ=255)
 aGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACAcc                                                                                                                                 >  2:233561/1‑139 (MQ=255)
   cATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGc                                                                                                                               <  2:63954/139‑1 (MQ=255)
   cATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGc                                                                                                                               <  1:233561/139‑1 (MQ=255)
     ttCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGcca                                                                                                                             <  2:534925/139‑1 (MQ=255)
     ttCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGcca                                                                                                                             <  2:38438/139‑1 (MQ=255)
       ctctGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGccacc                                                                                                                           <  2:473324/139‑1 (MQ=255)
          tgttgtGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCaccacc                                                                                                                        <  1:49210/139‑1 (MQ=255)
                ggAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGcca                                                                                                                             <  1:32688/128‑1 (MQ=255)
                ggAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGcca                                                                                                                             >  2:32688/1‑128 (MQ=255)
                    ggTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATgg                                                                                                              >  1:551491/1‑139 (MQ=255)
                      tAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGa                                                                                                            >  1:61036/1‑139 (MQ=255)
                        aaGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGaaa                                                                                                          <  2:61036/139‑1 (MQ=255)
                           cTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATgg                                                                                                       <  1:306787/139‑1 (MQ=255)
                              aTAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGcca                                                                                                                             <  1:273219/114‑1 (MQ=255)
                              aTAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGcca                                                                                                                             >  2:273219/1‑114 (MQ=255)
                              aTAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCgg                                                                                                    >  2:585317/1‑139 (MQ=255)
                              aTAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCgg                                                                                                    >  1:466061/1‑139 (MQ=255)
                                   gTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTAc                                                                                               >  2:416950/1‑139 (MQ=255)
                                     cGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTActg                                                                                             >  1:555475/1‑139 (MQ=255)
                                     cGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTActg                                                                                             <  2:609423/139‑1 (MQ=255)
                                       tACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTActgct                                                                                           >  1:113216/1‑139 (MQ=255)
                                       tACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTActgct                                                                                           >  1:528032/1‑139 (MQ=255)
                                        aCTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGccac                                                                                                                            <  2:428370/105‑1 (MQ=255)
                                        aCTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGccac                                                                                                                            >  1:428370/1‑105 (MQ=255)
                                         cTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGacac                                                                                                                                  <  1:571761/98‑1 (MQ=255)
                                         cTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGacac                                                                                                                                  >  2:571761/1‑98 (MQ=255)
                                         cTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTgg                                                                                         <  2:503205/139‑1 (MQ=255)
                                          tGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGc                                                                                        <  2:113216/139‑1 (MQ=255)
                                          tGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGc                                                                                        >  2:10607/1‑139 (MQ=255)
                                           ggAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGcc                                                                                       <  2:142786/139‑1 (MQ=255)
                                              aaaTGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGAt                                                                                    >  2:84381/1‑139 (MQ=255)
                                              aaaTGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGAt                                                                                    >  2:225716/1‑139 (MQ=255)
                                                 tGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTcc                                                                                 >  1:388646/1‑139 (MQ=255)
                                                 tGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTcc                                                                                 >  2:556760/1‑139 (MQ=255)
                                                      gACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGcc                                                                                                                              >  1:344992/1‑89 (MQ=255)
                                                      gACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGcc                                                                                                                              <  2:344992/89‑1 (MQ=255)
                                                       aCGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGcc                                                                           <  2:389663/139‑1 (MQ=255)
                                                        cGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCg                                                                          <  1:10607/139‑1 (MQ=255)
                                                        cGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCg                                                                          <  1:531969/139‑1 (MQ=255)
                                                              gATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCa                                                                    <  2:528032/139‑1 (MQ=255)
                                                               aTGTGCTGACCGGACATACAGCCAATAATGCCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAAGGGTGACTCCGCCACCACCGTATCCATGTGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTAAg                                                                   <  1:8956/139‑1 (MQ=255)
                                                                   gCTGACCGGACACACAGTCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGaaa                                                                                                          >  2:381936/1‑96 (MQ=255)
                                                                   gCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGaaa                                                                                                          <  1:381936/96‑1 (MQ=255)
                                                                   gCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCt                                                                        <  2:121655/130‑1 (MQ=255)
                                                                   gCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCt                                                                        >  1:121655/1‑130 (MQ=255)
                                                                    cTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTa                                                              >  2:597952/1‑139 (MQ=255)
                                                                     tGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTAc                                                             >  1:42430/1‑139 (MQ=255)
                                                                       aCCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCg                                                                                                     <  1:47907/97‑1 (MQ=255)
                                                                       aCCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCg                                                                                                     >  2:47907/1‑97 (MQ=255)
                                                                          ggACACACAGCCACTAATACCCGCGTGGATTATGGCGGAACGCTGGTTGTCCGCAACGGTGACACCGccac                                                                                                                            <  2:372365/71‑1 (MQ=255)
                                                                          ggACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGccac                                                                                                                            >  1:372365/1‑71 (MQ=255)
                                                                           gACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAg                                                       <  2:2136/139‑1 (MQ=255)
                                                                            acacacAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGc                                                      <  1:150805/139‑1 (MQ=255)
                                                                                      aCTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGcc                                                                           >  1:85720/1‑108 (MQ=255)
                                                                                      aCTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGcc                                                                           <  2:85720/108‑1 (MQ=255)
                                                                                          aTACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAg                                                       >  2:603183/1‑124 (MQ=255)
                                                                                          aTACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAg                                                       <  1:603183/124‑1 (MQ=255)
                                                                                          aTACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCAtt                                        <  1:585317/139‑1 (MQ=255)
                                                                                           tACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAANGNTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGCAAGGCATTc                                       <  2:326289/139‑1 (MQ=255)
                                                                                                gcgTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACAcc                                                                                                                                 <  1:192192/44‑1 (MQ=255)
                                                                                                gcgTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACAcc                                                                                                                                 >  2:192192/1‑44 (MQ=255)
                                                                                                gcgTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTAt                                  >  1:160648/1‑139 (MQ=255)
                                                                                                 cgTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATc                                 >  2:312962/1‑139 (MQ=255)
                                                                                                  gTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCg                                >  2:424697/1‑139 (MQ=255)
                                                                                                  gTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCg                                >  2:539497/1‑139 (MQ=255)
                                                                                                   tGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGc                                           >  1:99857/1‑127 (MQ=255)
                                                                                                   tGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGc                                           <  2:99857/127‑1 (MQ=255)
                                                                                                   tGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCGGAGCGACGGAAAGGCATTCAGTATCGGa                              <  1:597952/139‑1 (MQ=255)
                                                                                                    ggATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGt                                                                      >  1:59456/1‑99 (MQ=255)
                                                                                                    ggATGATGGCGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGt                                                                      <  2:59456/99‑1 (MQ=255)
                                                                                                             cGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAg                     >  1:145083/1‑139 (MQ=255)
                                                                                                             cGGAACGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAg                     >  1:281690/1‑139 (MQ=255)
                                                                                                                aaCGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAGATGGCGgtg                                                                                                  >  1:166578/1‑59 (MQ=39)
                                                                                                                aaCGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAGATGGCGgtg                                                                                                  <  2:166578/59‑1 (MQ=39)
                                                                                                                aaCGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTc                                                                     <  2:94368/88‑1 (MQ=255)
                                                                                                                aaCGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTc                                                                     >  1:94368/1‑88 (MQ=255)
                                                                                                                aaCGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCa                                                                    <  1:284866/89‑1 (MQ=255)
                                                                                                                aaCGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCa                                                                    >  2:284866/1‑89 (MQ=255)
                                                                                                                aaCGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCg                  >  1:391320/1‑139 (MQ=255)
                                                                                                                  cGCTGGATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGa                <  2:488099/139‑1 (MQ=255)
                                                                                                                       gATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTAt                                  <  1:236880/116‑1 (MQ=255)
                                                                                                                       gATGTCCGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTAt                                  >  2:236880/1‑116 (MQ=255)
                                                                                                                             cGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGt                                                                      <  2:463391/74‑1 (MQ=255)
                                                                                                                             cGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGt                                                                      >  1:463391/1‑74 (MQ=255)
                                                                                                                             cGCAACGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGc     <  2:466061/139‑1 (MQ=255)
                                                                                                                               cAACGGTGACACCGCCACCACCGTATCCATGGGAAATgg                                                                                                       <  2:509957/39‑1 (MQ=37)
                                                                                                                               cAACGGTGACACCGCCACCACCGTATCCATGGGAAATgg                                                                                                       >  1:509957/1‑39 (MQ=37)
                                                                                                                                aaCGGTGACACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTgg  <  1:84381/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
GAGCATTCTCTGTTGTGGAGGGTAAAGCTGATAATGTCGTACTGGAAAATGGCGGACGCCTGGATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGG  >  NC_000913/2072524‑2072790

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