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

GCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTT  >  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)
      gtCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCggg                                                                                                                                             <  1:1041831/150‑1 (MQ=255)
      gtCCGGGAAAATATCACCCTGGGCGGTCAGGCCGGTTATGCCCAGTGCGTCAGCGGCTGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCTTCGCCTCTCTGTGGTCCCGGTCATGATGACCggg                                                                                                                                             <  2:600848/150‑1 (MQ=255)
       tCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCgggc                                                                                                                                            >  2:431323/1‑149 (MQ=255)
         cGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCGGGAcc                                                                                                                                          <  1:122985/150‑1 (MQ=255)
         cGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCGGGAcc                                                                                                                                          <  1:123612/150‑1 (MQ=255)
             aaaaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCGCGGTCATCATGACCGGGACCCGGa                                                                                                                                      <  2:106524/150‑1 (MQ=255)
                atatCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTGATCATGACCGGGACCCGGACCg                                                                                                                                   >  1:1309698/1‑150 (MQ=255)
                   tCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGgcg                                                                                                                                <  1:739159/150‑1 (MQ=255)
                    cACCCTGGGCGTGCAGGCCGGTTATGCCCACAGCGGCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGgcgc                                                                                                                               <  1:62265/150‑1 (MQ=255)
                      cccTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGAGAGATCCATCGCCTCTCTGTGGTCCCGGTCTTCATTACCGGGACCCGGACCGGCGCaa                                                                                                                             >  1:1134861/1‑150 (MQ=255)
                       ccTGGGCGATCAGGCCGGTTATGCCCACAGCTTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCCGTCATAAAGACCAGGGCCCCGGAACGCGCaaa                                                                                                                            >  2:1194043/1‑149 (MQ=255)
                        cTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACCGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACg                                                                                                                           <  1:161423/150‑1 (MQ=255)
                          gggCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGa                                                                                                                         <  2:747512/150‑1 (MQ=255)
                            gCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGCCCGGCGCAACGGATc                                                                                                                       >  1:582765/1‑150 (MQ=255)
                            gCGGTCAGGCCGGTTATGCCCTCAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCAGACTGAATGTGTCGTTCTGACAGCACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCTGCGCAACGGATc                                                                                                                       <  1:37772/150‑1 (MQ=255)
                                   ggCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCCACGGATCTTCAAcc                                                                                                                >  2:1108720/1‑149 (MQ=255)
                                   ggCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACATTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGCGTCTTCAACg                                                                                                                >  2:218726/1‑150 (MQ=255)
                                     ccGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGGCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGcc                                                                                                              <  2:179414/150‑1 (MQ=255)
                                      cGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGACCCATCGCCTCTCTGTGGTCCCGGGCGTGATGCCCGGGCCCCGGCCCGGGGCCACGGATCTTCAACGccc                                                                                                             >  1:314866/1‑149 (MQ=255)
                                       ggTTAGGCCCACTGCGGCAGCGGCAGCAGCGCTTATGGTTTAATCTGTCAGGCCACACTTAATGTGACCTTCTCCCCGAACCATCGCCTCTCTGTGGCCCCGGTCACCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCcac                                                                                                            <  2:531465/150‑1 (MQ=255)
                                        gTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCCGGCGTGAAGTACGGGCCCCGGCCCGGCGCCGCGGGGCTTCAACGccccc                                                                                                           >  1:1079108/1‑147 (MQ=255)
                                          tATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGGCAGGCCACACTGAAGGTGACCTTCTGACACAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACAtt                                                                                                         <  1:366175/150‑1 (MQ=255)
                                            tGCCCACAGCGTCAGCGGCAGCAGCGTTTATGTGTATAACGTTCTGGCCACAATTAATGTGACATTCTGACAGACCCATCGCCTCTCTGTGGTCCCGGTCATCATGACGGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCg                                                                                                       <  2:1189780/150‑1 (MQ=255)
                                            tGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGGCATCATGACCGGGGCCCGGGCCGGGGCAACGGCTCGTCCACGcccccttgg                                                                                                       >  1:971146/1‑143 (MQ=255)
                                              cccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGGCCCGGCCCGGCGCCACGGTTCTTCAACGCCCCATTCggg                                                                                                     >  1:441406/1‑148 (MQ=255)
                                              cccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCt                                                                                                     >  2:1560/1‑150 (MQ=255)
                                               ccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGCACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGCCCCGGACCGGCGCCACGGATCTTCAACGCCACCTTCGCTg                                                                                                    >  2:290730/1‑150 (MQ=255)
                                                cacaGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTgg                                                                                                   >  2:89830/1‑150 (MQ=255)
                                                cacaGCGACAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTTACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTgg                                                                                                   <  2:1285151/150‑1 (MQ=255)
                                                    gCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCAtt                                                                                               >  2:563477/1‑150 (MQ=255)
                                                     cGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCCACGCCACCTTCGCTGGCATTa                                                                                              >  2:637068/1‑150 (MQ=255)
                                                            ggtagcagcGCTGAAGGGTATAACGGTCAGGCCCCACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGCCCCGGCCCGGCGCAACGGAGCTTCAACGCCACAGTCGCGGGCAGTACCAAGaa                                                                                       >  2:165334/4‑150 (MQ=255)
                                                                 cagcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGCCCGGGACCCGGACCGGCGCCACGGCTCTTCACCGCCACCTTCGGTGGCCTTAACAATacactga                                                                                  >  1:444271/1‑144 (MQ=255)
                                                                    ggcttatgggAATAAGGGGCAGCCCACACTCAAGGTTCCCTCCTCCCCGACCCAACGCCTCTCTGTGGCCCCGGTCAACATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGatat                                                                               <  1:1312607/143‑1 (MQ=255)
                                                                        gAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCTGCGCAACGGATCTTAAGCGCCACAGTCGCTGGGATTAACAAAAACATGATATtaat                                                                           >  2:196086/1‑150 (MQ=255)
                                                                         aaGGGTATACCGGTCAGGCCACACTGAATGGGACCTTCTGACAGAACCATCGCCGCTCTGTGGTCCCGGGCATCATGACCGGGACCCGGAACGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtcatc                                                                          <  2:85316/150‑1 (MQ=255)
                                                                               ataAGGGTCAGGCCACACTGAATGTGACGTTCTCACACAACCAGCGCCTCTCTGTGGCCCCGGCCATCATGACCGGTACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAg                                                                    <  1:427632/150‑1 (MQ=255)
                                                                                   cGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGGCCCGGACCGGCGCAACGGATCTTCAACGCCACCTTGGCTGGCATTAACAATAACATGATATTCATCACGGAGTGAc                                                                >  1:391943/1‑150 (MQ=255)
                                                                                    ggTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACt                                                               <  1:537264/150‑1 (MQ=255)
                                                                                           ccACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTAc                                                        <  1:322170/150‑1 (MQ=255)
                                                                                             acacTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAg                                                      <  1:561669/150‑1 (MQ=255)
                                                                                             acacTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGCACCGGCGCAACGGAACTTCAACCCCACAATCGCTGTCATTAACAATAACATGAAATTCATCACGGAGTGACTAAGTTACac                                                      >  2:53059/1‑149 (MQ=255)
                                                                                                 tGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGGCCCGGTCATCTTGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAg                                                  <  2:306598/150‑1 (MQ=255)
                                                                                                    aTGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCg                                               <  1:1280317/150‑1 (MQ=255)
                                                                                                    aTGTCACCTTCTCACACAACCATCGCCTCTCTGTGGACCCGGTCATCAGGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCg                                               <  1:605254/150‑1 (MQ=255)
                                                                                                     ggtcaCCTTCTGACCCACCCATCGCCCCTCTGGGGCCCCGGCCTACATGCCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCgg                                              <  1:1259951/146‑1 (MQ=255)
                                                                                                                    caaCCATCGCGGCTCTGGGGTCCGGGTATACATGCCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctgct                               <  2:1188030/149‑1 (MQ=255)
                                                                                                                        cATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATATTCAACGCCACATTAGCTGGCATTAACAATAACATGATATTCATAACGGAGTGACTATGTTTCAGATAGTCTGCGCGCAGATACAGCTGATc                           >  2:1216885/1‑150 (MQ=255)
                                                                                                                          tCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGc                        <  2:368577/150‑1 (MQ=255)
                                                                                                                             ccgcgctGTGGCCCCGGGCATCATGACCGGTACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAgg                     <  1:431323/145‑1 (MQ=255)
                                                                                                                             ccTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAgg                     <  1:214919/150‑1 (MQ=255)
                                                                                                                             ccTCTCTGTGGTCCCGGGCATCTTGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCTCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAgg                     <  2:33883/150‑1 (MQ=255)
                                                                                                                                ttctGTGGTCCCGTTCATCATGCCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGAtt                  <  2:1140633/149‑1 (MQ=255)
                                                                                                                                 tctGTGGCCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCGATCAGCAGGAtca                 <  1:579056/150‑3 (MQ=255)
                                                                                                                                   tgtgGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGc               <  2:22956/150‑1 (MQ=255)
                                                                                                                                   tgtgGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGc               <  1:738514/150‑1 (MQ=255)
                                                                                                                                   tgtgGTCCCCGTGATCATGACCGGGACCCGGACCGGCGCAACGGATCTTAAACGCCACATTCGCTGGCATTAAAAATAAAATGATATTAAAAACCGAGTGACTATTTTACACATAGTACGCGCGCTGATCCAGCTTATCGCAGGATTTGc               >  2:954575/1‑150 (MQ=255)
                                                                                                                                   ggggggcccGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGc               <  1:340266/144‑1 (MQ=255)
                                                                                                                                    gtgGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGAGTTGcc              >  1:873414/1‑150 (MQ=255)
                                                                                                                                     tgGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCa             >  2:940483/1‑150 (MQ=255)
                                                                                                                                      ggTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATCATATTCATCAAGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAt            >  2:726586/1‑150 (MQ=255)
                                                                                                                                       gTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGATGATCGCAGGATTTGCCAtt           >  1:521659/1‑150 (MQ=255)
                                                                                                                                        tCCCGGTCATCATGACCGGGACCCGGACCGGCGCTACGGATCTTCAACGCCACAGTCGCTGGCATTCACATTAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGGAGGATTTGCCCttc          >  2:390452/1‑150 (MQ=255)
                                                                                                                                          ccGGTCATCATGACCGGGACCCGGACCGGTGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttctt        >  2:924249/1‑150 (MQ=255)
                                                                                                                                           cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttcttc       >  1:20067/1‑150 (MQ=255)
                                                                                                                                           cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttcttc       <  2:675797/150‑1 (MQ=255)
                                                                                                                                           cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAgtcttc       >  2:870453/1‑150 (MQ=255)
                                                                                                                                              tcatcaTGACGGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATAAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGcc    <  2:486275/150‑1 (MQ=255)
                                                                                                                                              acatcatGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGcc    <  1:1289161/149‑1 (MQ=255)
                                                                                                                                                atcatGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTACGCCtt  >  2:38802/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTT  >  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: