Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,073,094 T→C L519P (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,0940TC98.8% 314.4 / ‑4.6 84L519P (CTG→CCG) fluCP4‑44 prophage, antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (44/39);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGA  >  NC_000913/2072960‑2073228
                                                                                                                                      |                                                                                                                                      
ccGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCacgccgacg                                                                                                                                    >  1:96119/1‑139 (MQ=255)
 cGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCacgccgacgc                                                                                                                                   >  1:585365/1‑139 (MQ=255)
   tGAAGGCGATGCACTCATGCAGGGCGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCAAACTCACCCTTAAAGCCGTCAACCTGAATGAAGGCacgccgacgatg                                                                                                                                 <  2:313246/139‑1 (MQ=255)
     aaGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGaa                                                                                                                               >  2:584306/1‑139 (MQ=255)
     aaGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGaa                                                                                                                               >  1:6078/1‑139 (MQ=255)
     aaGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACACAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACCCTGaa                                                                                                                               >  2:777143/1‑139 (MQ=255)
        gCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGa                                                                                                                            >  2:513206/1‑139 (MQ=255)
         cGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAAc                                                                                                                              <  1:680700/136‑1 (MQ=255)
         cGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAAc                                                                                                                              >  2:680700/1‑136 (MQ=255)
                     gCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTcacc                                                                                                               <  1:84759/139‑1 (MQ=255)
                     gCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTcacc                                                                                                               <  1:154716/139‑1 (MQ=255)
                      cAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTcacca                                                                                                              >  1:65588/1‑139 (MQ=255)
                      cAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTcacca                                                                                                              >  1:159396/1‑139 (MQ=255)
                          gAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGa                                                                                                          >  1:180300/1‑139 (MQ=255)
                           aGGCTCTCTCACCGGTAACGGCAGCTCGGAACAATCACGACGTGGCACACTCACTGTCAACAACACCACACTCACTCAGACAGCCGGCAACCTGAATGAAGGCACGCCGACGCTAAACGACAGTACCGTCACCACggct                                                                                                         >  2:305135/1‑137 (MQ=255)
                           aGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGAt                                                                                                         >  1:332891/1‑139 (MQ=255)
                             gCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGt                                                                                                       >  2:48456/1‑139 (MQ=255)
                             gCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGt                                                                                                       >  1:108432/1‑139 (MQ=255)
                                   tcACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGc                                                                                                 >  1:135958/1‑139 (MQ=255)
                                     aCCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTc                                                                                               >  2:712998/1‑139 (MQ=255)
                                     aCCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTc                                                                                               >  1:533324/1‑139 (MQ=255)
                                       cGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAg                                                                                             >  2:22874/1‑139 (MQ=255)
                                         gTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCATGGATGTCATTGCTCAgcg                                                                                           <  2:585365/139‑1 (MQ=255)
                                         gTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcg                                                                                           <  2:180300/139‑1 (MQ=255)
                                         gTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcg                                                                                           <  2:630228/139‑1 (MQ=255)
                                         gTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcg                                                                                           <  2:742902/139‑1 (MQ=255)
                                         gTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcg                                                                                           <  2:762407/139‑1 (MQ=255)
                                             cGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGt                                                                                       >  1:689245/1‑139 (MQ=255)
                                              ggCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTa                                                                                      <  2:6078/139‑1 (MQ=255)
                                               gCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTAc                                                                                     <  1:498930/139‑1 (MQ=255)
                                               gCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTAc                                                                                     <  2:135958/139‑1 (MQ=255)
                                                cAGCGTGGAAAAATCAGGAAGTGGCACACCCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcg                                                                                           >  2:14331/1‑132 (MQ=255)
                                                cAGCGTGGAAAAATCAGGAAGTGGCACACCCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcg                                                                                           <  1:14331/132‑1 (MQ=255)
                                                   cGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGa                                                                                                          >  1:58108/1‑114 (MQ=255)
                                                   cGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGa                                                                                                          <  2:58108/114‑1 (MQ=255)
                                                   cGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcg                                                                                           >  2:653898/1‑129 (MQ=255)
                                                   cGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcg                                                                                           <  1:653898/129‑1 (MQ=255)
                                                   cGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGcc                                                                                 <  2:108432/139‑1 (MQ=255)
                                                       gAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGa                                                                             >  2:384275/1‑139 (MQ=255)
                                                            aTCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTg                                                                        <  2:331722/139‑1 (MQ=255)
                                                              cAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGAc                                                                      <  2:96119/139‑1 (MQ=255)
                                                               aGGAAGTGGCAGACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCTACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcg                                                                                           >  2:467090/1‑117 (MQ=255)
                                                               aGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcg                                                                                           <  1:467090/117‑1 (MQ=255)
                                                                 gAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTAc                                                                                                                      >  1:717995/1‑88 (MQ=255)
                                                                 gAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTAc                                                                                                                      >  1:774900/1‑88 (MQ=255)
                                                                 gAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTAc                                                                                                                      <  2:774900/88‑1 (MQ=255)
                                                                 gAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTAc                                                                                                                      <  2:717995/88‑1 (MQ=255)
                                                                 gAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCgg                                                                   >  1:402157/1‑139 (MQ=255)
                                                                     tGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagc                                                               <  2:445328/139‑1 (MQ=255)
                                                                      ggCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              <  1:384275/139‑1 (MQ=255)
                                                                      ggCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              <  1:351877/139‑1 (MQ=255)
                                                                      ggCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              <  2:159396/139‑1 (MQ=255)
                                                                      ggCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              <  2:59007/139‑1 (MQ=255)
                                                                               actGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGa                                                     >  1:364862/1‑139 (MQ=255)
                                                                               actGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTCAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGa                                                     >  2:497576/1‑139 (MQ=255)
                                                                                ctGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGaa                                                    <  2:65588/139‑1 (MQ=255)
                                                                                 tGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAAc                                                   >  2:228461/1‑139 (MQ=255)
                                                                                     aGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGt                                                                                                                   >  1:474898/1‑71 (MQ=255)
                                                                                     aGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGt                                                                                                                   <  2:474898/71‑1 (MQ=255)
                                                                                        aaCACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCa                                            <  2:313142/139‑1 (MQ=255)
                                                                                           accacACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAg                                                                                                                         >  1:352418/1‑59 (MQ=255)
                                                                                           accacACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAg                                                                                                                         <  2:352418/59‑1 (MQ=255)
                                                                                             cacaCTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              <  2:170086/116‑1 (MQ=255)
                                                                                             cacaCTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              >  1:170086/1‑116 (MQ=255)
                                                                                             cacaCTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGAc                                       >  1:215324/1‑139 (MQ=255)
                                                                                             cacaCTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGAc                                       >  2:336019/1‑139 (MQ=255)
                                                                                             cacaCTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGAc                                       >  2:617122/1‑139 (MQ=255)
                                                                                               cacTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGAccc                                     >  2:698434/1‑139 (MQ=255)
                                                                                                 cTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCa                                   >  2:221128/1‑139 (MQ=255)
                                                                                                    aCCCAGAAAGCCGTCAACCTGAATGAAGGCacgccgacgctg                                                                                                                                 <  2:33000/42‑1 (MQ=255)
                                                                                                    aCCCAGAAAGCCGTCAACCTGAATGAAGGCacgccgacgctg                                                                                                                                 >  1:33000/1‑42 (MQ=255)
                                                                                                    aCCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACg                                                                                                            >  2:665466/1‑63 (MQ=255)
                                                                                                    aCCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACg                                                                                                            <  1:665466/63‑1 (MQ=255)
                                                                                                    aCCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGa                                >  2:404116/1‑139 (MQ=255)
                                                                                                    aCCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGa                                >  1:273969/1‑139 (MQ=255)
                                                                                                    aCCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGa                                >  1:524883/1‑139 (MQ=255)
                                                                                                                 tCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTg                                                                        >  1:374657/1‑86 (MQ=255)
                                                                                                                 tCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTg                                                                        <  2:374657/86‑1 (MQ=255)
                                                                                                                    aCCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              <  2:115966/93‑1 (MQ=255)
                                                                                                                    aCCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              >  1:115966/1‑93 (MQ=255)
                                                                                                                    aCCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTc                >  1:248315/1‑139 (MQ=255)
                                                                                                                    aCCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTc                <  1:22874/139‑1 (MQ=255)
                                                                                                                        gaatgaaGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGt            <  2:273969/139‑1 (MQ=255)
                                                                                                                         aatgaaGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTg           <  2:600305/139‑1 (MQ=255)
                                                                                                                          atgaaGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGc          <  1:513206/139‑1 (MQ=255)
                                                                                                                           tgaaGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGcc         >  2:40565/1‑139 (MQ=255)
                                                                                                                               ggCACGCCGACGCTGAACGACAGTACCGTCACCACGGa                                                                                                          >  1:49963/1‑38 (MQ=255)
                                                                                                                               ggCACGCCGACGCTGAACGACAGTACCGTCACCACGGa                                                                                                          <  2:49963/38‑1 (MQ=255)
                                                                                                                                gCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTa                                                                                      <  1:533858/57‑1 (MQ=255)
                                                                                                                                gCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTa                                                                                      >  2:533858/1‑57 (MQ=255)
                                                                                                                                  acgccgacgctgAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGa  <  2:768680/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
CCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGA  >  NC_000913/2072960‑2073228

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: