Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,074,691 C→G 5.7% intergenic (+33/‑88) flu → / → yeeR CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter/CP4‑44 prophage; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,074,6910CG 5.7% 302.5 / 8.4 88intergenic (+33/‑88)flu/yeeRCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter/CP4‑44 prophage; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base C (43/40);  new base G (2/3);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

GTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCG  >  NC_000913/2074551‑2074834
                                                                                                                                            |                                                                                                                                               
gtCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCg                                                                                                                                          >  7:333486/1‑148 (MQ=255)
   cGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCGGGa                                                                                                                                       <  6:187511/148‑1 (MQ=255)
    gggAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGCCCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGAc                                                                                                                                      >  6:365291/1‑148 (MQ=255)
       aaaaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGGGGTCCCGGTCATGATGACCGGGAccac                                                                                                                                   <  3:24102/148‑3 (MQ=255)
          atatCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGAc                                                                                                                                >  2:160194/1‑148 (MQ=255)
            atCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCg                                                                                                                              >  1:348148/1‑148 (MQ=255)
             tCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCgg                                                                                                                             <  4:151648/148‑1 (MQ=255)
               aCCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGgcg                                                                                                                           <  8:333329/148‑1 (MQ=255)
                cccTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGgcg                                                                                                                           >  3:87913/1‑147 (MQ=255)
                      gCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGa                                                                                                                    <  6:362619/148‑1 (MQ=255)
                          tCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGAcc                                                                                                                               <  3:312820/133‑1 (MQ=255)
                          tCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGAcc                                                                                                                               >  4:312820/1‑133 (MQ=255)
                           cAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCtt                                                                                                                <  1:170248/147‑1 (MQ=255)
                           cAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCtt                                                                                                                >  7:349081/1‑147 (MQ=255)
                           cAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCtt                                                                                                                <  6:206941/147‑1 (MQ=255)
                            aGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCa                                                                                                              <  8:358785/148‑1 (MQ=255)
                              gCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCg                                                                                                                              <  6:3501/130‑1 (MQ=255)
                              gCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCg                                                                                                                              >  5:3501/1‑130 (MQ=255)
                                cGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGc                                                                                                          >  8:401346/1‑148 (MQ=255)
                                   ttATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGCCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCca                                                                                                        >  8:172702/1‑147 (MQ=255)
                                       gCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTc                                                                                                   >  5:257291/1‑148 (MQ=255)
                                            caGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTgg                                                                                              <  7:401346/148‑1 (MQ=255)
                                             aGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTgg                                                                                              >  1:225867/1‑147 (MQ=255)
                                                 tCAGCGGCAGCAGCGCTGAAGTGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCATCTCTGTGGTCCCGGTCATCATTACCGTGACCCGGACAGTCGCAACGGATCTTCAACGCCACATTCGCTGGCATTa                                                                                         <  7:416772/148‑1 (MQ=255)
                                                 tCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTa                                                                                         <  5:345179/148‑1 (MQ=255)
                                                 tCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTa                                                                                         <  1:264257/148‑1 (MQ=255)
                                                 tCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTa                                                                                         <  7:416770/148‑1 (MQ=255)
                                                  cAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTaa                                                                                        <  2:402770/148‑1 (MQ=255)
                                                   aGCGGCAGCAGCTTTGTAGGGTATAACGGTCAGGGCACAATGAATGTGACCTTCTGACAGATCCATCGCCTCTCTGTTGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCaggtat                                                                                       <  8:333486/148‑6 (MQ=255)
                                                     cGGCTGCTGCGCTGAAGGGTATAACGGTCAGGCCACCCTGAAGGTGACCTTCTGACAGACCCATCGCCTCTCTGTGGACCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACaa                                                                                     <  7:161609/148‑1 (MQ=255)
                                                      ggCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTa                                                                                         <  2:42058/143‑1 (MQ=255)
                                                      ggCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTa                                                                                         >  1:42058/1‑143 (MQ=255)
                                                        cagcagCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGCCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCTCTGGCATTAACAATaa                                                                                  <  4:258476/148‑1 (MQ=255)
                                                          gcagcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGCCCGGGCCCCGGCCCGGCGCCACGGATCTTCAACGCCCCAGTCGCTGGCATTAACAATAACa                                                                                >  4:237142/1‑148 (MQ=255)
                                                           cagcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACAt                                                                               <  8:349081/148‑1 (MQ=255)
                                                            agcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATg                                                                              <  2:252681/148‑1 (MQ=255)
                                                              cgcTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGa                                                                             >  7:154117/1‑147 (MQ=255)
                                                                cTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGTGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGatat                                                                          <  1:63770/148‑1 (MQ=255)
                                                                    aGGGTATAACGGTCAGGCCACACTGAACGTGACCTTCTGACAGAACCATCGCCGCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtcat                                                                      <  6:235747/148‑1 (MQ=255)
                                                                         ataACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACgg                                                                 >  7:22899/1‑148 (MQ=255)
                                                                          taACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGa                                                                <  6:257291/148‑1 (MQ=255)
                                                                           aaCGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGat                                                                            <  2:405344/135‑1 (MQ=255)
                                                                           aaCGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGat                                                                            >  1:405344/1‑135 (MQ=255)
                                                                           aaCGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAg                                                               >  3:192773/1‑148 (MQ=255)
                                                                           aaCGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGGCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtcatc                                                                     >  7:376498/1‑142 (MQ=255)
                                                                           aaCGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGGAATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtcatc                                                                     <  8:376498/142‑1 (MQ=255)
                                                                            aCGGTCAGGCCACACTGAATGTGACCTTATGACAGAACCATCGACTGTCTGTGGTCCAGGTCATCATGAACGGTAGCCTGACCTGCGCAAAGGATCTTCAAAGACACATTCGGTGGCATTAAGAATAAAATCATATTAATAACTGAGt                                                              >  5:167115/1‑148 (MQ=255)
                                                                               gTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGAc                                                           >  1:227029/1‑148 (MQ=255)
                                                                                     ccACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGCCCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACGATGt                                                      >  2:181148/1‑147 (MQ=255)
                                                                                     ccACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGt                                                      >  5:245134/1‑147 (MQ=255)
                                                                                       acacTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTAc                                                   <  8:154117/148‑1 (MQ=255)
                                                                                        cacTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTa                                                    >  4:205261/1‑146 (MQ=255)
                                                                                            gAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGa                                                >  4:95801/1‑146 (MQ=255)
                                                                                             aaTGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcatc                                                                                                                                                 >  6:74716/1‑48 (MQ=255)
                                                                                             aaTGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcatc                                                                                                                                                 <  5:74716/48‑1 (MQ=255)
                                                                                              aTGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGGGATCATGACCGGGACCCGGGCCGGGGGCACGGGGCGTGAACGCCCCAGTCGCGGGGGATAAAAATAACCTGATATTTATCCCGGAGTGGCGTTGTTTAAGATAGt                                            >  8:55853/1‑148 (MQ=255)
                                                                                                 tgACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCgg                                         >  3:149878/1‑148 (MQ=255)
                                                                                                    ccTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGgcgc                                      >  1:20424/1‑148 (MQ=255)
                                                                                                      ttCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGgcgcgc                                    >  4:333698/1‑148 (MQ=255)
                                                                                                           aCAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGAt                                <  3:95801/147‑1 (MQ=255)
                                                                                                           aCAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATc                               <  3:99841/148‑1 (MQ=255)
                                                                                                            cAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATcc                              >  4:335359/1‑148 (MQ=255)
                                                                                                            cAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATcc                              >  4:134009/1‑148 (MQ=255)
                                                                                                            cAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATcc                              >  6:11745/1‑148 (MQ=255)
                                                                                                             aGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCct                             <  3:333698/148‑1 (MQ=255)
                                                                                                              gcgccctCGCCTCTGTGGGGGCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCACCGCCACATTCGCTGGCATTAACAATAACCTGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctg                            <  6:354514/142‑1 (MQ=255)
                                                                                                              gAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctg                            >  1:190178/1‑148 (MQ=255)
                                                                                                              gAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctg                            <  3:134009/148‑1 (MQ=255)
                                                                                                                   aTCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAATGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtcatc                                                                     <  8:133218/102‑1 (MQ=255)
                                                                                                                   aTCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAATGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtcatc                                                                     >  7:133218/1‑102 (MQ=255)
                                                                                                                   aTCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCACCGCCCCCTTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGCCTATGTGACAGAGAGTCGGCGCGCTGATCctgc                           >  5:354514/1‑144 (MQ=255)
                                                                                                                   aTCGCCTCTCTGTGGTCCCGGGCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTCACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGAt                       <  3:237142/148‑1 (MQ=255)
                                                                                                                    tCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATc                      <  8:22899/148‑1 (MQ=255)
                                                                                                                     cGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCg                     >  1:43337/1‑148 (MQ=255)
                                                                                                                       ccTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCg                     >  5:66538/1‑146 (MQ=255)
                                                                                                                        ctctctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAg                  >  1:115054/1‑148 (MQ=255)
                                                                                                                         tctctGTGGTCCCGGTCATCATGACCGGGACCCGGCCCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAg                  >  5:395185/1‑147 (MQ=255)
                                                                                                                          ctctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGa                                                                <  3:236107/100‑1 (MQ=255)
                                                                                                                          ctctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGa                                                                >  4:236107/1‑100 (MQ=255)
                                                                                                                          ctctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGa                >  1:211717/1‑148 (MQ=255)
                                                                                                                              gtgGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTg            >  3:375651/1‑148 (MQ=255)
                                                                                                                                ggTCCCGGTCATGATGCCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctgc                           >  5:375028/1‑131 (MQ=255)
                                                                                                                                ggTCCCGGTCATGATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctgc                           <  6:375028/131‑1 (MQ=255)
                                                                                                                                 gTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGcc          <  6:66538/147‑1 (MQ=255)
                                                                                                                                 gTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCACTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGcc          <  4:150472/147‑1 (MQ=255)
                                                                                                                                 gTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTAGATGGCATTAACAATAACATGATATTCATAACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGATGATCGCAGGATTTGCCa         >  7:298258/1‑148 (MQ=255)
                                                                                                                                     cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct     <  2:20424/148‑1 (MQ=255)
                                                                                                                                     cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct     <  2:190178/148‑1 (MQ=255)
                                                                                                                                      ggTCATGATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttctt    >  5:411806/1‑148 (MQ=255)
                                                                                                                                       gTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACTGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttcttc   <  3:21681/148‑1 (MQ=255)
                                                                                                                                        tcatcaTGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCg  >  2:260229/1‑148 (MQ=255)
                                                                                                                                            catGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGcc                                                                                                         <  5:308007/41‑1 (MQ=255)
                                                                                                                                            catGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGcc                                                                                                         >  6:308007/1‑41 (MQ=255)
                                                                                                                                            |                                                                                                                                               
GTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCG  >  NC_000913/2074551‑2074834

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