Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,540,613 A→G intergenic (‑1112/‑344) nanC ← / → fimB N‑acetylneuraminic acid outer membrane channel/regulator for fimA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,540,6130AG98.9% 314.9 / NA 89intergenic (‑1112/‑344)nanC/fimBN‑acetylneuraminic acid outer membrane channel/regulator for fimA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (38/50);  minor base T (1/0);  total (39/50)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATA  >  NC_000913/4540494‑4540742
                                                                                                                       |                                                                                                                                 
aGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAAAGTCATATGTGAGTCGGTCAGTTGAACGGATCATAAATTTTg                                                                                                                      >  2:285704/1‑133 (MQ=255)
        aaGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAAttt                                                                                                              >  1:111484/1‑133 (MQ=255)
         aGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAAtttt                                                                                                             >  1:386989/1‑133 (MQ=255)
         aGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAAtttt                                                                                                             >  1:563997/1‑133 (MQ=255)
         aGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAATTGAACGTATAGTAAATTTTGATGAAtttt                                                                                                             >  1:179568/1‑133 (MQ=255)
             ttAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGGCGGCTAGTGGAAGGGAAATTAAATTATGCTAAATTTTTTAt                                                                                                         >  2:167727/1‑133 (MQ=255)
             ttAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTAt                                                                                                         >  2:596582/1‑133 (MQ=255)
             ttAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTAt                                                                                                         >  2:259931/1‑133 (MQ=255)
             ttAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTAt                                                                                                         >  2:421709/1‑133 (MQ=255)
             ttAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTAt                                                                                                         <  2:310406/133‑1 (MQ=255)
             ttAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGAAATTAAATTTTGCTCAATTTTTTAt                                                                                                         >  2:212236/1‑133 (MQ=255)
                gCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGGTGAATTTTTTATGtt                                                                                                      >  1:273412/1‑133 (MQ=255)
                    gATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGAtt                                                                                                  <  1:596582/133‑1 (MQ=255)
                     ataataGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGAttt                                                                                                 >  1:339406/1‑133 (MQ=255)
                     ataataGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTAGCTGAATTTTTTATGTTGAttt                                                                                                 >  1:279469/1‑133 (MQ=255)
                      taataGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGGCTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGAtttt                                                                                                <  2:403916/133‑1 (MQ=255)
                         taGCCACTCAGAAATCACAGCGCTCCATGATATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACt                                                                                             <  1:285704/133‑1 (MQ=255)
                         taGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACt                                                                                             <  2:455904/133‑1 (MQ=255)
                         taGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACt                                                                                             <  2:623272/133‑1 (MQ=255)
                           gCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTg                                                                                           >  1:81513/1‑133 (MQ=255)
                              aCTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTa                                                                                        >  1:490984/1‑133 (MQ=255)
                              aCTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTa                                                                                        >  1:338023/1‑133 (MQ=255)
                              aCTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTa                                                                                        <  1:154722/133‑1 (MQ=255)
                              aCTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTa                                                                                        >  1:326591/1‑133 (MQ=255)
                               cTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTa                                                                                        <  2:583469/132‑1 (MQ=255)
                               cTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTa                                                                                        >  1:583469/1‑132 (MQ=255)
                               cTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTAc                                                                                       <  2:111484/133‑1 (MQ=255)
                               cTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTAc                                                                                       <  1:259931/133‑1 (MQ=255)
                               cTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGATGAATTTATTAGGTCGATTTTGCTTGTTAc                                                                                       >  1:184528/1‑133 (MQ=255)
                                tAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACa                                                                                      >  2:102848/1‑133 (MQ=255)
                                tAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACa                                                                                      <  2:563997/133‑1 (MQ=255)
                                tAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACa                                                                                      <  1:560711/133‑1 (MQ=255)
                                tAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACa                                                                                      <  2:386989/133‑1 (MQ=255)
                                tAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACa                                                                                      >  2:37824/1‑133 (MQ=255)
                                 aaGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAg                                                                                     <  1:578161/133‑1 (MQ=255)
                                  aGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGa                                                                                    >  1:2788/1‑133 (MQ=255)
                                   gAAATTACTGGGCTCGATGAAATAGCCAATTTTTGGCAAATAGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGttctgtaa                                                                                   <  2:179568/133‑7 (MQ=255)
                                       ttACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACata                                                                               <  2:339406/133‑1 (MQ=255)
                                         aCTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAGATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATc                                                                             <  1:102848/133‑1 (MQ=255)
                                         aCTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATc                                                                             <  1:501169/133‑1 (MQ=255)
                                         aCTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATc                                                                             <  2:490984/133‑1 (MQ=255)
                                         aCTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATc                                                                             <  2:563709/133‑1 (MQ=255)
                                           tGCGCTCCATTAAATCGCCATTTTCTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATcac                                                                           <  2:170426/133‑1 (MQ=255)
                                           tGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATcac                                                                           <  1:421709/133‑1 (MQ=255)
                                           tGCGCTCCATAAAATTGCCATTTTATTGAAAAAGGAGTTGAATAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATcac                                                                           <  2:265835/133‑1 (MQ=255)
                                            gggcTCCATTCAATATCCATTTTTAGGCTAATGGAGTTGACAAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATcaca                                                                          <  2:289310/131‑1 (MQ=255)
                                            gcgcTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATcaca                                                                          <  2:81513/133‑1 (MQ=255)
                                            gcgcTCCACGAAATAACCATTTCGTGGCAAGTGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATcaca                                                                          <  1:161194/133‑1 (MQ=255)
                                             cgcTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACAt                                                                         <  1:269090/133‑1 (MQ=255)
                                             cgcTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACAt                                                                         <  2:418630/133‑1 (MQ=255)
                                             cgcTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACAt                                                                         >  2:572213/1‑133 (MQ=255)
                                              gcTCCATGAACTAGCTATTTTGTGCCAACTGGCGTTGACTATTAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATg                                                                        <  2:184528/133‑1 (MQ=255)
                                              gcTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATcac                                                                           >  1:505740/1‑130 (MQ=255)
                                              gcTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATcac                                                                           <  2:505740/130‑1 (MQ=255)
                                              gcTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATg                                                                        >  2:522120/1‑133 (MQ=255)
                                              gcTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATg                                                                        <  1:572213/133‑1 (MQ=255)
                                              gcTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATg                                                                        <  2:21293/133‑1 (MQ=255)
                                              gcTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATg                                                                        <  1:144231/133‑1 (MQ=255)
                                              gcTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTGGATTTTAGTTGGTATAGAACAGATCACATg                                                                        >  1:180324/1‑133 (MQ=255)
                                              gcTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACAAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATg                                                                        >  2:21931/1‑133 (MQ=255)
                                               cTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATcaca                                                                          >  2:590108/1‑130 (MQ=255)
                                               cTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATcaca                                                                          <  1:590108/130‑1 (MQ=255)
                                                            ccATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGa                                                          >  2:24710/1‑133 (MQ=255)
                                                            ccATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGa                                                          >  1:12254/1‑133 (MQ=255)
                                                            ccATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGa                                                          >  1:382048/1‑133 (MQ=255)
                                                                tttGTGGCCAATGGACTTTACTCATAATGTGATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAg                                                      <  1:171852/133‑1 (MQ=255)
                                                                  tgtgGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACGTGTTACAGCACATATCGCATGATATGTAGAAAAGATTAGtt                                                    >  2:171852/1‑133 (MQ=255)
                                                                     ggCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCa                                                 <  2:12254/133‑1 (MQ=255)
                                                                     ggCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCa                                                 <  2:558781/133‑1 (MQ=255)
                                                                         aaTGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCAAATGATATATAGGTATGATTAGGTGCATTaa                                             >  1:179352/1‑133 (MQ=255)
                                                                             gAGTTGACTAATAATGTCGTATCTGAGACGACTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAatgat                                         <  1:185203/133‑1 (MQ=255)
                                                                             gAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAatgat                                         <  2:338023/133‑1 (MQ=255)
                                                                             gAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAatgat                                         <  2:326591/133‑1 (MQ=255)
                                                                                ttGACTAATAATGTCTTATGTGGGACGGCTCGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAg                                      <  1:212236/133‑1 (MQ=255)
                                                                                          atGTCATACGTGAGAAGGCTAGGTGAACGGCTATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGttattatt                            <  1:160681/133‑1 (MQ=255)
                                                                                             tCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGa                         >  2:396632/1‑133 (MQ=255)
                                                                                             tCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGa                         >  2:543705/1‑133 (MQ=255)
                                                                                             tCATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGAATTAATGATGAGGGTTATTATTAGa                         >  2:489663/1‑133 (MQ=255)
                                                                                              cATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATTGATAAGATTAGTAGCAGTAATGATGGGGGTTATTATTAGAt                        >  2:160681/1‑133 (MQ=255)
                                                                                              cATATGTGAGACGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGAt                        <  1:525795/133‑1 (MQ=255)
                                                                                              cATATGTGAGACGGCAAGTTGAACCGATTTTATATTTTGCTCAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGAt                        <  2:180324/133‑1 (MQ=255)
                                                                                                      agaCGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATcc                >  2:575627/1‑133 (MQ=255)
                                                                                                      agaCGGCTAGTTGAACGGATATAAAATTTTCCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATcc                <  1:167727/133‑1 (MQ=255)
                                                                                                        aCGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGa              <  1:24710/133‑1 (MQ=255)
                                                                                                        aCGGCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGa              <  1:522120/133‑1 (MQ=255)
                                                                                                           gCTAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTg           >  2:606380/1‑133 (MQ=255)
                                                                                                             tAGTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGAt         <  2:104163/133‑1 (MQ=255)
                                                                                                               gTTGAACGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATaa       <  1:56270/133‑1 (MQ=255)
                                                                                                                    aCGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAatata  <  1:575627/133‑1 (MQ=255)
                                                                                                                    aCGGATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAatata  <  1:606380/133‑1 (MQ=255)
                                                                                                                       |                                                                                                                                 
AGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATA  >  NC_000913/4540494‑4540742

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: