Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,074,691 C→G 8.2% 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 8.2% 173.7 / 6.5 70intergenic (+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 (31/32);  new base G (2/4);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTT  >  NC_000913/2074545‑2074838
                                                                                                                                                  |                                                                                                                                                     
gCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGAGGTCCCGGTCATCATg                                                                                                                                                    >  2:296811/1‑150 (MQ=255)
 cccGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGa                                                                                                                                                   >  2:579056/1‑150 (MQ=255)
   cGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGCACCATCGCCTCTCTGTGGTCCCCGTCATCATGccc                                                                                                                                                 >  1:784110/1‑150 (MQ=255)
      gtCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACC‑ggg                                                                                                                                             <  1:1041831/150‑1 (MQ=255)
      gtCCGGGAAAATATCACCCTGGGCGGTCAGGCCGGTTATGCCCAGTGCGTCAGCGGCTGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCTTCGCCTCTCTGTGGTCCCGGTCATGATGACC‑ggg                                                                                                                                             <  2:600848/150‑1 (MQ=255)
       tCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑gggc                                                                                                                                            >  2:431323/1‑149 (MQ=255)
         cGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACC‑GGGAcc                                                                                                                                          <  1:122985/150‑1 (MQ=255)
         cGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACC‑GGGAcc                                                                                                                                          <  1:123612/150‑1 (MQ=255)
             aaaaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCGCGGTCATCATGACC‑GGGACCCGGa                                                                                                                                      <  2:106524/150‑1 (MQ=255)
                atatCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTGATCATGACC‑GGGACCCGGACCg                                                                                                                                   >  1:1309698/1‑150 (MQ=255)
                   tCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGgcg                                                                                                                                <  1:739159/150‑1 (MQ=255)
                    cACCCTGGGCGTGCAGGCCGGTTATGCCCACAGCGGCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGgcgc                                                                                                                               <  1:62265/150‑1 (MQ=255)
                      cccTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGAGAGATCCATCGCCTCTCTGTGGTCCCGGTCTTCATTACC‑GGGACCCGGACCGGCGCaa                                                                                                                             >  1:1134861/1‑150 (MQ=255)
                       ccTGGGCGATCAGGCCGGTTATGCCCACAGCTTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCCGTCATAAAGACCAGGGCCCCGGAAC‑GCGCaaa                                                                                                                            >  2:1194043/1‑149 (MQ=255)
                        cTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACCGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACg                                                                                                                           <  1:161423/150‑1 (MQ=255)
                          gggCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGa                                                                                                                         <  2:747512/150‑1 (MQ=255)
                            gCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGCCCGGCGCAACGGATc                                                                                                                       >  1:582765/1‑150 (MQ=255)
                            gCGGTCAGGCCGGTTATGCCCTCAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCAGACTGAATGTGTCGTTCTGACAGCACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCTGCGCAACGGATc                                                                                                                       <  1:37772/150‑1 (MQ=255)
                                   ggCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCCACGGATCTTCAAcc                                                                                                                >  2:1108720/1‑149 (MQ=255)
                                   ggCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACATTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGCGTCTTCAACg                                                                                                                >  2:218726/1‑150 (MQ=255)
                                     ccGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGGCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGcc                                                                                                              <  2:179414/150‑1 (MQ=255)
                                      cGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGACCCATCGCCTCTCTGTGGTCCCGGGCGTGATGCCC‑GGGCCCCGGCCCGGGGCCACGGATCTTCAACGccc                                                                                                             >  1:314866/1‑149 (MQ=255)
                                       ggTTAGGCCCACTGCGGCAGCGGCAGCAGCGCTTATGGTTTAATCTGTCAGGCCACACTTAATGTGACCTTCTCCCCGAACCATCGCCTCTCTGTGGCCCCGGTCACCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCcac                                                                                                            <  2:531465/150‑1 (MQ=255)
                                        gTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCCGGCGTGAAGTAC‑GGGCCCCGGCCCGGCGCCGCGGGGCTTCAACGccccc                                                                                                           >  1:1079108/1‑147 (MQ=255)
                                          tATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGGCAGGCCACACTGAAGGTGACCTTCTGACACAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACAtt                                                                                                         <  1:366175/150‑1 (MQ=255)
                                            tGCCCACAGCGTCAGCGGCAGCAGCGTTTATGTGTATAACGTTCTGGCCACAATTAATGTGACATTCTGACAGACCCATCGCCTCTCTGTGGTCCCGGTCATCATGACG‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCg                                                                                                       <  2:1189780/150‑1 (MQ=255)
                                            tGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGGCATCATGACC‑GGGGCCCGGGCCGGGGCAACGGCTCGTCCACGcccccttgg                                                                                                       >  1:971146/1‑143 (MQ=255)
                                              cccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGGCCCGGCCCGGCGCCACGGTTCTTCAACGCCCCATTCggg                                                                                                     >  1:441406/1‑148 (MQ=255)
                                              cccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCt                                                                                                     >  2:1560/1‑150 (MQ=255)
                                               ccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGCACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGCCCCGGACCGGCGCCACGGATCTTCAACGCCACCTTCGCTg                                                                                                    >  2:290730/1‑150 (MQ=255)
                                                cacaGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTgg                                                                                                   >  2:89830/1‑150 (MQ=255)
                                                cacaGCGACAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTTACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTgg                                                                                                   <  2:1285151/150‑1 (MQ=255)
                                                    gCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCAtt                                                                                               >  2:563477/1‑150 (MQ=255)
                                                     cGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCCACGCCACCTTCGCTGGCATTa                                                                                              >  2:637068/1‑150 (MQ=255)
                                                            ggtagcagcGCTGAAGGGTATAACGGTCAGGCCCCACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGCCCCGGCCCGGCGCAACGGAGCTTCAACGCCACAGTCGCGGGCAGTACCAAGaa                                                                                       >  2:165334/4‑150 (MQ=255)
                                                                 cagcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGCCC‑GGGACCCGGACCGGCGCCACGGCTCTTCACCGCCACCTTCGGTGGCCTTAACAATacactga                                                                                  >  1:444271/1‑144 (MQ=255)
                                                                    ggcttatgggAATAAGGGGCAGCCCACACTCAAGGTTCCCTCCTCCCCGACCCAACGCCTCTCTGTGGCCCCGGTCAACATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGatat                                                                               <  1:1312607/143‑1 (MQ=255)
                                                                        gAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCTGCGCAACGGATCTTAAGCGCCACAGTCGCTGGGATTAACAAAAACATGATATtaat                                                                           >  2:196086/1‑150 (MQ=255)
                                                                         aaGGGTATACCGGTCAGGCCACACTGAATGGGACCTTCTGACAGAACCATCGCCGCTCTGTGGTCCCGGGCATCATGACC‑GGGACCCGGAACGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtcatc                                                                          <  2:85316/150‑1 (MQ=255)
                                                                               ataAGGGTCAGGCCACACTGAATGTGACGTTCTCACACAACCAGCGCCTCTCTGTGGCCCCGGCCATCATGACC‑GGTACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAg                                                                    <  1:427632/150‑1 (MQ=255)
                                                                                   cGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGGCCCGGACCGGCGCAACGGATCTTCAACGCCACCTTGGCTGGCATTAACAATAACATGATATTCATCACGGAGTGAc                                                                >  1:391943/1‑150 (MQ=255)
                                                                                    ggTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACt                                                               <  1:537264/150‑1 (MQ=255)
                                                                                           ccACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTAc                                                        <  1:322170/150‑1 (MQ=255)
                                                                                             acacTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAg                                                      <  1:561669/150‑1 (MQ=255)
                                                                                             acacTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGCACCGGCGCAACGGAACTTCAACCCCACAATCGCTGTCATTAACAATAACATGAAATTCATCACGGAGTGACTAAGTTACac                                                      >  2:53059/1‑149 (MQ=255)
                                                                                                 tGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGGCCCGGTCATCTTGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAg                                                  <  2:306598/150‑1 (MQ=255)
                                                                                                    aTGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCg                                               <  1:1280317/150‑1 (MQ=255)
                                                                                                    aTGTCACCTTCTCACACAACCATCGCCTCTCTGTGGACCCGGTCATCAGGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCg                                               <  1:605254/150‑1 (MQ=255)
                                                                                                     ggtcaCCTTCTGACCCACCCATCGCCCCTCTGGGGCCCCGGCCTACATGCCC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCgg                                              <  1:1259951/146‑1 (MQ=255)
                                                                                                                    caaCCATCGCGGCTCTGGGGTCCGGGTATACATGCCC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctgct                               <  2:1188030/149‑1 (MQ=255)
                                                                                                                        cATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATATTCAACGCCACATTAGCTGGCATTAACAATAACATGATATTCATAACGGAGTGACTATGTTTCAGATAGTCTGCGCGCAGATACAGCTGATc                           >  2:1216885/1‑150 (MQ=255)
                                                                                                                          tCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑Gc                        <  2:368577/150‑1 (MQ=255)
                                                                                                                             ccgcgctGTGGCCCCGGGCATCATGACC‑GGTACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAgg                     <  1:431323/145‑1 (MQ=255)
                                                                                                                             ccTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAgg                     <  1:214919/150‑1 (MQ=255)
                                                                                                                             ccTCTCTGTGGTCCCGGGCATCTTGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCTCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAgg                     <  2:33883/150‑1 (MQ=255)
                                                                                                                                ttctGTGGTCCCGTTCATCATGCCC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGAtt                  <  2:1140633/149‑1 (MQ=255)
                                                                                                                                 tctGTGGCCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGC‑GATCAGCAGGAtca                 <  1:579056/150‑3 (MQ=255)
                                                                                                                                   tgtgGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGc               <  2:22956/150‑1 (MQ=255)
                                                                                                                                   tgtgGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGc               <  1:738514/150‑1 (MQ=255)
                                                                                                                                   tgtgGTCCCCGTGATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTAAACGCCACATTCGCTGGCATTAAAAATAAAATGATATTAAAAACCGAGTGACTATTTTACACATAGTACGCGCGCTGATCCAGCTTATC‑GCAGGATTTGc               >  2:954575/1‑150 (MQ=255)
                                                                                                                                   ggggggcccGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGc               <  1:340266/144‑1 (MQ=255)
                                                                                                                                    gtgGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGAGTTGcc              >  1:873414/1‑150 (MQ=255)
                                                                                                                                     tgGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCa             >  2:940483/1‑150 (MQ=255)
                                                                                                                                      ggTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATCATATTCATCAAGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAt            >  2:726586/1‑150 (MQ=255)
                                                                                                                                       gTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGATGATC‑GCAGGATTTGCCAtt           >  1:521659/1‑150 (MQ=255)
                                                                                                                                        tCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCTACGGATCTTCAACGCCACAGTCGCTGGCATTCACATTAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GGAGGATTTGCCCttc          >  2:390452/1‑150 (MQ=255)
                                                                                                                                          ccGGTCATCATGACC‑GGGACCCGGACCGGTGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAttctt        >  2:924249/1‑150 (MQ=255)
                                                                                                                                           cGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAttcttc       >  1:20067/1‑150 (MQ=255)
                                                                                                                                           cGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAttcttc       <  2:675797/150‑1 (MQ=255)
                                                                                                                                           cGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAgtcttc       >  2:870453/1‑150 (MQ=255)
                                                                                                                                              tcatcaTGACG‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATAAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGcc    <  2:486275/150‑1 (MQ=255)
                                                                                                                                              acatcatGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGcc    <  1:1289161/149‑1 (MQ=255)
                                                                                                                                                atcatGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTACGCCtt  >  2:38802/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACC‑GGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTT  >  NC_000913/2074545‑2074838

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