Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,378,977 T→A 6.9% intergenic (‑33/+30) ampC ← / ← frdD penicillin‑binding protein; beta‑lactamase, intrinsically weak/fumarate reductase (anaerobic), membrane anchor subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,9770TA 6.9% 331.4 / 8.3 118intergenic (‑33/+30)ampC/frdDpenicillin‑binding protein; beta‑lactamase, intrinsically weak/fumarate reductase (anaerobic), membrane anchor subunit
Reads supporting (aligned to +/- strand):  ref base T (54/55);  new base A (5/4);  total (59/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 = 1.00e+00

Only 100 of 118 total aligned reads displayed.

TGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGG  >  NC_000913/4378829‑4379124
                                                                                                                                                    |                                                                                                                                                   
tgttgCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTGGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATa                                                                                                                                                    <  2:1205897/150‑1 (MQ=255)
      tctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCggg                                                                                                                                              >  1:874132/1‑150 (MQ=255)
        ttttaGCGGGGTAATTGTGCGATGCAGAATATCGTTGATTTGTTTAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCCCCTATGGCGGGcc                                                                                                                                            <  1:831616/146‑1 (MQ=255)
        tataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGcc                                                                                                                                            >  2:928430/1‑150 (MQ=255)
        tataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGAAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGcc                                                                                                                                            <  2:742371/150‑1 (MQ=255)
         ataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGcag                                                                                                                                           >  1:323090/1‑148 (MQ=255)
         ataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCg                                                                                                                                           <  2:1724963/150‑1 (MQ=255)
          tcaGCGGGGTAATTGTGCGATGCACAATATAGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGCCgt                                                                                                                                          <  2:492924/148‑2 (MQ=255)
          taAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGCCgt                                                                                                                                          <  2:264306/150‑2 (MQ=255)
           aaGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCggt                                                                                                                                         >  1:837486/1‑149 (MQ=255)
            aGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTACATACAAAACTGCCCGCCATAGTCTGTCCGGAt                                                                                                                                        >  2:1558120/1‑150 (MQ=255)
             gCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                                                                                                                                       >  1:1415913/1‑150 (MQ=255)
              gggggTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGCGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTCAACATCGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAttt                                                                                                                                      <  2:2151170/149‑1 (MQ=255)
              cGGGGTATTGGTGCGATCCACAATAGCGTTGTTTTGTGGAGGGGCAGCAAATGTGGAGCATCAGGCGGTAATTAATAAGGCGCAGAGCGGCGTTTTGAACATAGGGGCTGGTATCCATACAAAACGGCCCGCCATAGGCGGGCCGGAttt                                                                                                                                      <  1:1942042/150‑1 (MQ=255)
               aacgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTa                                                                                                                                     >  1:1441244/7‑150 (MQ=255)
                 tttAATTGTTCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGTTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACa                                                                                                                                   <  2:487905/148‑1 (MQ=255)
                    aaTTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTg                                                                                                                                <  2:815419/150‑1 (MQ=255)
                     aTTGTGCGCTGCACAATATCGTTTATTTGTGGAGGGGCAGCAAATGTGGAGCAAGAGGCGGAAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTgg                                                                                                                               <  1:558027/150‑1 (MQ=255)
                         tgCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATTAGGCGCTGAGCGTCGTTTTGAACCTAGGGTCTGGTTTCCATACAAAACGGCCCGCCCTAGGCTGGCCGGATTTACATTCTCGAt                                                                                                                           >  1:350910/1‑150 (MQ=255)
                         tgCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATACGGTCTGGTTTCCATACAAAACGTCCCGCCATAGGCTGGCCGGATTTACATTGGCGAt                                                                                                                           >  2:413491/1‑150 (MQ=255)
                            gATGCACAATATCGTTGTTTTGTTGTGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGGCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCg                                                                                                                        <  2:1053465/150‑1 (MQ=255)
                              tGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGtt                                                                                                                      <  2:475033/150‑1 (MQ=255)
                              tGCACAATATCGTTGATTTGTGGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGtt                                                                                                                      <  2:210087/150‑1 (MQ=255)
                               gCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTa                                                                                                                     >  1:1770337/1‑150 (MQ=255)
                                cacaATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGACCGCCATAGGCGGGACGGATTTACATTGGCGATGCGTTAg                                                                                                                    >  1:984822/1‑150 (MQ=255)
                                  aaaTATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGAt                                                                                                                  <  2:909701/149‑1 (MQ=255)
                                  aaaTATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGGTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGAt                                                                                                                  <  1:475707/149‑1 (MQ=255)
                                   aaTATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGAtt                                                                                                                 <  1:1429441/150‑1 (MQ=255)
                                      atCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTa                                                                                                              <  2:1192725/150‑1 (MQ=255)
                                         gTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACg                                                                                                           <  2:1042933/150‑1 (MQ=255)
                                           tGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGac                                                                                                         <  1:2174085/150‑1 (MQ=255)
                                              tgtgttgAGGGGCAGCAAATGGGCCGCAAGAGGCGGTAATTAATAAGTCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACAcc                                                                                                      <  2:125310/148‑1 (MQ=255)
                                               ttgttgAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCa                                                                                                     >  1:1790929/1‑150 (MQ=255)
                                                 gttgAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGAGGCGTTAGATTGTAACGACACCAAt                                                                                                   >  2:173174/1‑150 (MQ=255)
                                                  ttgAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATc                                                                                                  <  2:531676/150‑1 (MQ=255)
                                                  ttgAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATc                                                                                                  <  1:303615/150‑1 (MQ=255)
                                                    gAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAg                                                                                                <  1:1311115/150‑1 (MQ=255)
                                                     aGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGc                                                                                               >  1:2027657/1‑150 (MQ=255)
                                                      ggggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCg                                                                                              <  2:616534/150‑1 (MQ=255)
                                                       gggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGCCACCAATCAGCGt                                                                                             >  2:1965781/1‑150 (MQ=255)
                                                        ggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGGCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTg                                                                                            <  1:1238961/150‑1 (MQ=255)
                                                        ggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTg                                                                                            <  1:742799/150‑1 (MQ=255)
                                                        ggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGACCGCCATAGGCGGGGCGGATTTACATTGGCGATGCGTTAGATTGTAACGCCACCAATCAGCGTg                                                                                            >  2:343096/1‑150 (MQ=255)
                                                         gcagcaAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGa                                                                                           <  1:1085932/150‑1 (MQ=255)
                                                           agcaAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACa                                                                                         <  1:590898/150‑1 (MQ=255)
                                                             caAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAAc                                                                                       <  2:637317/150‑1 (MQ=255)
                                                             caAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAAc                                                                                       >  2:1939073/1‑150 (MQ=255)
                                                                   tgGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGCCACCAATCAGCGTGACAACTGTccg                                                                                 >  1:1115925/1‑148 (MQ=255)
                                                                      aGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGAt                                                                              >  2:510759/1‑150 (MQ=255)
                                                                      aGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGAt                                                                              <  2:510684/150‑1 (MQ=255)
                                                                           gcggCGGTAATTAATAAGGCGCAGCGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTTAGGCGGGCCGGATTTACTTTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATagcag                                                                         <  1:1047529/148‑1 (MQ=255)
                                                                           gagGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATagcag                                                                         <  1:2059075/150‑1 (MQ=255)
                                                                             ggCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGcc                                                                       >  2:1727750/1‑150 (MQ=255)
                                                                                    attaatAAGGCGCAGAGCGTCGATTTGAACATAGGGTCAGGTATCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGCTTGTAACGACACCAATCCGCGTGACAACTGTCAGGATAGCAGCCAGACCGt                                                                >  1:1079851/1‑150 (MQ=255)
                                                                                     ttaatAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTa                                                               >  2:1810479/1‑150 (MQ=255)
                                                                                      taataaGGCGCAGAGCGTCGTGTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAg                                                              <  1:2118877/150‑1 (MQ=255)
                                                                                        ataaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGAAGCGTTAGATTGTAACGCCACCACTCCGCGTGACCACTGTCAGGATAGCAGCCAGACCGTATaa                                                            >  1:422543/1‑150 (MQ=255)
                                                                                            ggCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAAcc                                                        >  1:357803/1‑150 (MQ=255)
                                                                                              cgcAGAGCGTCGTTTTGAACATAGGGTCTGTTTTCCATACAAAACGGCCCGCCAAATGCGGGCAGAAATTACATTGGCGATGCGTTAGATTGTAACGACAGCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACAGTATAAACCCCa                                                      >  2:1333369/1‑150 (MQ=255)
                                                                                                 agagCGTCGTTTTGAACATAGGGTCTGTTTTCAATACAAAACGGCCCTCCATAGGCGGGCAGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCAttt                                                   >  2:869784/1‑150 (MQ=255)
                                                                                                  gagCGTCGTTTTGAAAAAAAGTTATTATATCAATACAAAACGGACCGCCATAGAACAGAAGAAAATAAATTGGCTATTCGTTAGATTGTAACGACAACAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTg                                                  >  2:79650/1‑150 (MQ=255)
                                                                                                    gCGTCGTTTTGAACATAGGGTCTGGTATCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGcc                                                <  2:93672/150‑1 (MQ=255)
                                                                                                           tttGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGt                                         >  2:1430108/1‑150 (MQ=255)
                                                                                                            ttGAACATAGGGTCTGGTTTCAATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTa                                        <  2:1662274/150‑1 (MQ=255)
                                                                                                             tGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTAc                                       >  2:2020610/1‑150 (MQ=255)
                                                                                                             tGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTAc                                       >  2:1023573/1‑150 (MQ=255)
                                                                                                              gAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACg                                      >  2:2064129/1‑150 (MQ=255)
                                                                                                                aCATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTg                                    <  2:932674/150‑1 (MQ=255)
                                                                                                                 cATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGGAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGACAGCCCGTAGAACACCCATTTGCCCGCAGGTACGTgg                                   >  1:918815/1‑150 (MQ=255)
                                                                                                                   tAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGAt                                 >  2:445432/1‑150 (MQ=255)
                                                                                                                     gggTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGAttt                               <  1:323668/150‑1 (MQ=255)
                                                                                                                      ggTCTGGTTTCCATACAAAACGGCCCGCCATAGGCTGGCCGGATTTACATTGGCGATGCGATAGATTATAACGACACCAAACAGCGTGACAAATTTCAGGATAGAAGCCAGCCCGTAGAAAACCCATTTGCCCGCAGCTAGGAGGTtttt                              >  1:1330761/1‑150 (MQ=255)
                                                                                                                       gTCTTGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTc                             <  2:235676/150‑1 (MQ=255)
                                                                                                                        tCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGGCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCa                            <  1:627220/150‑1 (MQ=255)
                                                                                                                         cTGGTTTCCATACAAAACGGCCCGCCAAAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAg                           >  1:462838/1‑150 (MQ=255)
                                                                                                                         cTGGTATCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGATAGATTGTAAAGACACCAATCAGCGTGACATCTGTCAGGATAGCAGACAGACCGTAGAAAACCCAATTGCCAGAAGGTACGTGGATTTACaa                           >  1:406290/1‑149 (MQ=255)
                                                                                                                         cTGGATTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAg                           <  1:969379/150‑1 (MQ=255)
                                                                                                                           ggTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGAt                         >  1:673803/1‑150 (MQ=255)
                                                                                                                            gTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATc                        <  2:1806224/150‑1 (MQ=255)
                                                                                                                               tCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTg                     <  2:468679/150‑1 (MQ=255)
                                                                                                                               tCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTg                     <  1:1558120/150‑1 (MQ=255)
                                                                                                                               tCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTg                     <  1:413491/150‑1 (MQ=255)
                                                                                                                                 cATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCa                   >  1:787519/1‑150 (MQ=255)
                                                                                                                                 cATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCa                   >  1:1296891/1‑150 (MQ=255)
                                                                                                                                  aTACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCAt                  >  2:38727/1‑150 (MQ=255)
                                                                                                                                   tACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATc                 <  2:894922/150‑1 (MQ=255)
                                                                                                                                    aCAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcg                >  2:408657/1‑150 (MQ=255)
                                                                                                                                    aCAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcg                <  1:206363/150‑1 (MQ=255)
                                                                                                                                     cAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTACATTGTAACCACACCAATCAGCGTGCCAACGGTCAGGATAGCAGCCAGACCGTAGAAAAACCATTTACCCGCAGGTACTTGGATTATCAGATAGTGCATcgc               >  2:476840/1‑150 (MQ=255)
                                                                                                                                       aaaCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAACCCCATTTGCCCGCAGGTACGTGGAGTTTCAGATCGTGCATCGCgt             >  1:1517279/1‑150 (MQ=255)
                                                                                                                                       aaaCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgt             >  2:924425/1‑150 (MQ=255)
                                                                                                                                         aCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtgg           <  2:1812653/150‑1 (MQ=255)
                                                                                                                                         aCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtgg           <  1:1134120/150‑1 (MQ=255)
                                                                                                                                         aCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtgg           <  1:584225/150‑1 (MQ=255)
                                                                                                                                          cGGCCCGCCTATGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggg          >  2:496208/1‑149 (MQ=255)
                                                                                                                                           ggCCCGCCTATGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggtg         >  1:1840690/1‑150 (MQ=255)
                                                                                                                                            gCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGc        <  2:1632588/150‑1 (MQ=255)
                                                                                                                                               cGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATa     <  2:837486/150‑1 (MQ=255)
                                                                                                                                               cGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATa     >  2:108864/1‑150 (MQ=255)
                                                                                                                                                  cATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACgg  >  1:29768/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGG  >  NC_000913/4378829‑4379124

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: