Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,074,691 C→G 6.5% 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 6.5% 335.2 / 10.6 117intergenic (+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 (52/55);  new base G (4/4);  total (58/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

Only 100 of 123 total aligned reads displayed.

GAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTT  >  NC_000913/2074542‑2074838
                                                                                                                                                      |                                                                                                                                                     
  aGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTcatcat                                                                                                                                                     >  2:899516/1‑150 (MQ=255)
    cccGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGa                                                                                                                                                   <  1:96318/150‑1 (MQ=255)
          gccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATGATGACCGGGa                                                                                                                                             <  1:422497/149‑1 (MQ=255)
           ccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATGATGACCGGGAc                                                                                                                                            <  1:1117767/150‑1 (MQ=255)
            cGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATGATGACCGGGAcc                                                                                                                                           <  1:222546/150‑1 (MQ=255)
              ggAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGCCC‑Cg                                                                                                                                        >  2:452662/1‑150 (MQ=255)
              ggAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑Cg                                                                                                                                        <  2:1533096/150‑1 (MQ=255)
                  aaTATCACCGTGGGCGGTCAGGCCGGTTATGCCCACTTCGCCTTCGGCAGCAGCGCTTAAGGTTATAACGTTCAGGCCACTCTGAATGTGACCTTCTGACAGAACCATGGCCGGTCTGTGGG‑CCCGGTCATCATCACCGGGACC‑CGGAcc                                                                                                                                    <  2:828219/150‑1 (MQ=255)
                     atCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGCCC‑CGGCCCCgc                                                                                                                                 >  1:565369/1‑150 (MQ=255)
                        aCCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCa                                                                                                                              <  2:1025693/150‑1 (MQ=255)
                         cccTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCaa                                                                                                                             >  2:853631/1‑150 (MQ=255)
                             gggCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGa                                                                                                                         <  2:965171/150‑1 (MQ=255)
                             gggCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGa                                                                                                                         <  2:119262/150‑1 (MQ=255)
                              ggCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGAt                                                                                                                        >  2:1070772/1‑150 (MQ=255)
                               gCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGCTc                                                                                                                       >  2:545861/1‑150 (MQ=255)
                               gCGTGCAGGCCGGTTATGCCCACCGCGTCAGCGGCAGCCGCGCTGAAGGTTATAACGGTCAGGCCACACTCAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGC‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATc                                                                                                                       <  2:656677/150‑1 (MQ=255)
                                 gTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCtt                                                                                                                     >  2:633420/1‑150 (MQ=255)
                                 gTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCACCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCtt                                                                                                                     >  1:202516/1‑150 (MQ=255)
                                     aGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAAc                                                                                                                 <  2:480745/150‑1 (MQ=255)
                                      ggCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGCCCGGCCCCACCGGTCTTAAACg                                                                                                                >  1:1548845/1‑150 (MQ=255)
                                      ggCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACg                                                                                                                <  1:988842/150‑1 (MQ=255)
                                      ggCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACg                                                                                                                >  2:1581869/1‑150 (MQ=255)
                                       gCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGc                                                                                                               >  1:1115907/1‑150 (MQ=255)
                                        ccGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGGCAGCATGACCGGGACC‑CGGACCGGCGCCACGGATCTTCAACGcc                                                                                                              >  1:1252190/1‑150 (MQ=255)
                                         cGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCca                                                                                                             <  2:1625896/150‑1 (MQ=255)
                                           gTTATGCCCACAGCGCCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCcaca                                                                                                           <  1:349810/150‑1 (MQ=255)
                                              aTGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTc                                                                                                        <  1:466633/150‑1 (MQ=255)
                                                gCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGCCCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGCCCGGCGCAACGGATCTTCAACGCCCCATTCGc                                                                                                      >  2:1440758/1‑150 (MQ=255)
                                                gCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGCCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGc                                                                                                      >  1:1053491/1‑150 (MQ=255)
                                                 cccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCt                                                                                                     <  2:313895/150‑1 (MQ=255)
                                                 cccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCt                                                                                                     >  1:18045/1‑150 (MQ=255)
                                                  ccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGCCC‑CGGACCGGGGCAGCGGGTCTTCGACGGCACATTGGCTg                                                                                                    >  2:1588291/1‑150 (MQ=255)
                                                  ccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTg                                                                                                    >  2:1003936/1‑150 (MQ=255)
                                                  ccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCAGACAGAACCATCGCCACTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGAGCTTCAACGCCACATTCCCTg                                                                                                    >  1:1111610/1‑150 (MQ=255)
                                                   cacaGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCCCATTCGCTgg                                                                                                   >  2:369819/1‑150 (MQ=255)
                                                    acaGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGc                                                                                                  <  1:1039724/150‑1 (MQ=255)
                                                       gCGTCAGCGGCAGCAGCGCTTAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATTGCCTCTCTGTGGG‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCAtt                                                                                               <  1:774243/150‑1 (MQ=255)
                                                        cGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGATC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTa                                                                                              <  1:736290/150‑1 (MQ=255)
                                                          tCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACATAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAAc                                                                                            <  1:610219/150‑1 (MQ=255)
                                                          tCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCTACGGATCTTCAACGACACATTCGCTGGCATTaaa                                                                                            >  2:1165614/1‑149 (MQ=255)
                                                          tCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAAc                                                                                            >  2:1151596/1‑150 (MQ=255)
                                                            aGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACaa                                                                                          >  1:1223820/1‑150 (MQ=255)
                                                                gcagcaGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAAc                                                                                      <  1:37788/150‑1 (MQ=255)
                                                                    cagcGCTGAAGGTTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGa                                                                                  <  2:1596943/150‑1 (MQ=255)
                                                                       cgcTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGatat                                                                               >  1:1155709/1‑150 (MQ=255)
                                                                       cgcTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGatat                                                                               >  1:20668/1‑150 (MQ=255)
                                                                                ttataAGGGTCTGGCCACACTTAATGTTACCTTCTCACATAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACgg                                                                      <  2:565369/149‑1 (MQ=255)
                                                                                ttaaaagggTCAGCCCAAACTTATTGTGACCTTCTCACAGAACCATCGCCTCTATGTGGC‑CCCGGTCATCATCACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACgg                                                                      <  1:736255/143‑1 (MQ=255)
                                                                                 tataACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGa                                                                     <  1:975604/150‑1 (MQ=255)
                                                                                  ataACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCAGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACAGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGAAATTCATCACGGAg                                                                    >  2:1012293/1‑150 (MQ=255)
                                                                                   taACGGTCAGGCCACACTGAATGGGCCCTTCTGCCCGAACCTCCCCCTCTCGGTGGT‑CCCGGGCATCATGCCCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGt                                                                   <  1:823499/150‑1 (MQ=255)
                                                                                    aaCGGTCAGGCCACACTGAATCTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGGCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTg                                                                  <  2:1115888/150‑1 (MQ=255)
                                                                                      cGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCCACGGATCTTCAACGCCACATTAGCTGGCATTAACAATTACATGATATTCATCCCGGCGTGAc                                                                >  1:518068/1‑150 (MQ=255)
                                                                                          cAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATg                                                            <  1:304956/150‑1 (MQ=255)
                                                                                               cacaCTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACa                                                       >  2:1175134/1‑150 (MQ=255)
                                                                                               cacaCTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACa                                                       <  1:1070772/150‑1 (MQ=255)
                                                                                               caaagttaaTGTGACCTTCTCACAGAACCATCGCCTCTCTGGGGC‑CCCGGTCATCTTGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACa                                                       <  1:630555/143‑1 (MQ=255)
                                                                                                 cacTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATAATGACCGGGACC‑CGGACCGGCGCAACGGATCTTAAACGCCACAATAGCTGGCATTAACAATAACAAGATATTAATCACGGAGTGACTAAGTTACAGa                                                     >  1:1514960/1‑150 (MQ=255)
                                                                                                     gAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGt                                                 >  2:1194495/1‑150 (MQ=255)
                                                                                                        tgttACCTTCTGACTGAACCATCGCCTCTCTGTGTG‑CCCGGTCATCATTACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCgg                                              <  1:853631/150‑1 (MQ=255)
                                                                                                         gtgACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGgc                                             >  2:835427/1‑150 (MQ=255)
                                                                                                          tgACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGgcg                                            >  1:746126/1‑150 (MQ=255)
                                                                                                          tgACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGgcg                                            >  2:1641928/1‑150 (MQ=255)
                                                                                                          tgACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGgcg                                            <  2:1620637/150‑1 (MQ=255)
                                                                                                           gACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCCGGTAAACATGCCCGGGACC‑CGACCCGGCCCAACGGTTATTCACCCCCCCATTCCCGGCCATTAACAATAACATAATATTCAACACGGAGGGACTATGTTAAAGATAGACGgcc                                            >  2:1631713/1‑149 (MQ=255)
                                                                                                           gACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGgcgc                                           <  2:72824/150‑1 (MQ=255)
                                                                                                            aCCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATGTTCATCACGGAGTGACTATGTTACAGATAGTCGgcgcg                                          <  2:18045/150‑1 (MQ=255)
                                                                                                              cTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCt                                        <  1:899516/150‑1 (MQ=255)
                                                                                                              cTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACAGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCt                                        >  1:1200470/1‑150 (MQ=255)
                                                                                                               ttCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTg                                       <  1:1175134/150‑1 (MQ=255)
                                                                                                                  tGACAGAACCATCGCGTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATc                                    <  2:411799/150‑1 (MQ=255)
                                                                                                                   gACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATcc                                   >  1:899750/1‑150 (MQ=255)
                                                                                                                    aCAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACGTTCGCTGGCATTAACAATAACATGATATTCATCACGGATTGTCTATGTTTCAGATATTCGGCGCGCTGATCct                                  >  2:444262/1‑150 (MQ=255)
                                                                                                                    aCAGAACCATCGCCTCTCTGTGGT‑CCCGGCCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCct                                  <  2:1053491/150‑1 (MQ=255)
                                                                                                                     cAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctg                                 >  1:237880/1‑150 (MQ=255)
                                                                                                                     cAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctg                                 >  1:986191/1‑150 (MQ=255)
                                                                                                                     aaGAACCATCGCCTCTCTGTGTC‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctg                                 <  2:1332102/149‑1 (MQ=255)
                                                                                                                       gAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctgct                               <  2:758015/150‑1 (MQ=255)
                                                                                                                        aaCCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CTGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAAGAATAACATGATATTCCTCACGGAGTGACTATGTTACAGATAGTCGTCGCGCTGATCctgctg                              <  1:255391/150‑1 (MQ=255)
                                                                                                                          ccATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGAt                            <  2:900/150‑1 (MQ=255)
                                                                                                                          ccATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGAt                            >  2:200105/1‑150 (MQ=255)
                                                                                                                           cATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATc                           >  1:1027855/1‑150 (MQ=255)
                                                                                                                            aTCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑g                         <  2:1712479/150‑1 (MQ=255)
                                                                                                                                ccgctctGTGGT‑CCCGGTAATCATGACGGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAgg                     <  1:1099959/147‑1 (MQ=255)
                                                                                                                                ccTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAgg                     <  2:1115907/150‑1 (MQ=255)
                                                                                                                                 ctctctGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGa                    >  1:474491/1‑150 (MQ=255)
                                                                                                                                 ctctctGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGa                    >  2:1155818/1‑150 (MQ=255)
                                                                                                                                   ctctGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGC‑GATCAGCAGGAtc                  <  2:202516/150‑2 (MQ=255)
                                                                                                                                      tgtgGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTgg               <  1:36674/150‑2 (MQ=255)
                                                                                                                                       gtgGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGcc              <  2:1374602/150‑1 (MQ=255)
                                                                                                                                         ggT‑CCCGGTCATGATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAt            >  2:857897/1‑150 (MQ=255)
                                                                                                                                         ggT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAt            >  2:1699649/1‑150 (MQ=255)
                                                                                                                                         ggT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAt            >  2:354478/1‑150 (MQ=255)
                                                                                                                                          gT‑CCCGGTCATGATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAtt           >  2:1314802/1‑150 (MQ=255)
                                                                                                                                          gT‑CCCGGTCATCATGAACGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACAAGATATTCATCACGGAGTGACTATGTTACAGATAGTCTGCGCGCTGATCCAGATGATC‑TCACGATTTGCCAtt           >  1:781414/1‑150 (MQ=255)
                                                                                                                                           t‑CCCGGTCATGATGACCGGGACCACAGAGAGGCG‑‑ATGGTTCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAttct         >  1:119922/1‑150 (MQ=37)
                                                                                                                                             cccGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAttct         >  1:1527009/1‑150 (MQ=255)
                                                                                                                                              ccGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCAttctt        >  1:1518720/1‑150 (MQ=255)
                                                                                                                                                   catcatGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCt   >  2:1394586/1‑150 (MQ=255)
                                                                                                                                                    atcatGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCtt  >  2:877248/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGT‑CCCGGTCATCATGACCGGGACC‑CGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATC‑GCAGGATTTGCCATTCTTCGCCTT  >  NC_000913/2074542‑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: ‑