Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,345 T→C 5.8% pseudogene (954/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,3450TC 5.8% 316.8 / 6.8 104pseudogene (954/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (36/60);  new base C (3/3);  total (41/63)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

ACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGG  >  NC_000913/1466196‑1466494
                                                                                                                                                     |                                                                                                                                                     
aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                                                                                                                       <  2:2751925/150‑1 (MQ=255)
aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGAGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCAACGGCACACATATAAAt                                                                                                                                                       >  1:804397/1‑150 (MQ=17)
aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACAGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACTAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                                                                                                                       >  2:1924996/1‑150 (MQ=255)
 ccGGTGACAGCGCGACGGTGGATAACTAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGTTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGTGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                                                                                                                      <  1:2071199/150‑1 (MQ=255)
 ccGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                                                                                                                      <  2:907098/150‑1 (MQ=255)
 ccGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                                                                                                                      <  1:2552558/150‑1 (MQ=255)
   ggTGACAGCGCGACGGTGGATAACGAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGt                                                                                                                                                    <  2:2465772/150‑1 (MQ=255)
   ggTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGt                                                                                                                                                    <  1:1680826/150‑1 (MQ=255)
     tGACAGCGCGACGGTGTATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                                                                                                                  <  2:1062631/150‑1 (MQ=255)
     tGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                                                                                                                  <  2:1062787/150‑1 (MQ=255)
         agagcgACGGTGTATAACAAGGTCACCATGACCGTCTCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGTCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGAc                                                                                                                                              <  2:2034879/147‑1 (MQ=18)
          gcgcGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                                                                                                             <  2:1087158/150‑1 (MQ=255)
               accgTGGATAAAAAGGGCACCCAGACCCTCACCGATCCGGAGTCCATGGGTATCCAGATAGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCCCCAACACCAACGGTGGCACCGGCACCCAGATTAAACGTGATAACGCCCCg                                                                                                                                        >  2:1796323/4‑150 (MQ=17)
               aCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                                                                                                                        >  1:1629949/1‑150 (MQ=255)
               aCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                                                                                                                        >  1:2607496/1‑150 (MQ=255)
                  gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGTTGATGACGCCACGGCa                                                                                                                                     >  1:1907672/1‑150 (MQ=255)
                      aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                                                                                                                 <  1:2135113/150‑1 (MQ=255)
                          cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaa                                                                                                                             <  1:2575992/150‑1 (MQ=255)
                          cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaa                                                                                                                             >  2:1641076/1‑150 (MQ=255)
                          cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaa                                                                                                                             >  2:248275/1‑150 (MQ=255)
                          cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaa                                                                                                                             >  2:258717/1‑150 (MQ=255)
                           aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaac                                                                                                                            <  2:13229/150‑1 (MQ=255)
                           aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaac                                                                                                                            <  2:13228/150‑1 (MQ=255)
                           aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaac                                                                                                                            <  2:13213/150‑1 (MQ=255)
                           aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGAAGCCACGGCAaacaacaac                                                                                                                            >  2:394247/1‑150 (MQ=255)
                            aGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACg                                                                                                                           <  2:15174/150‑1 (MQ=255)
                                   cATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATGAATGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                                    >  1:1611301/1‑150 (MQ=255)
                                      gACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGCCCATCGTCAATAACGAAGGCGAGAGCACCAACACCAACGGTGGCACCGGCACCCACATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAaccac                                                                                                                 <  1:1277896/150‑1 (MQ=255)
                                           tCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                                                                                                            <  1:861417/150‑1 (MQ=255)
                                            cACCGATCCGGAGTCCATGGGTATCCAGATCGCCGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGa                                                                                                           <  2:1414895/150‑1 (MQ=255)
                                            cACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGa                                                                                                           <  2:1414866/150‑1 (MQ=255)
                                            cACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGa                                                                                                           <  2:1414887/150‑1 (MQ=255)
                                            cACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGa                                                                                                           <  2:992497/150‑1 (MQ=255)
                                              ccGATCCGGAGTCCATGGGTATCCAGTTAGACGGTGACAAGGCCATGGTGAATAACGAAGGCGAGAGCACCATCACCAACGGTGGGACCGGCACGCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg                                                                                                         >  1:1188563/1‑150 (MQ=14)
                                              ccGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg                                                                                                         >  1:1188610/1‑150 (MQ=255)
                                               cGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                                                        <  1:1644079/150‑1 (MQ=255)
                                                gATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGc                                                                                                       >  2:982169/1‑150 (MQ=255)
                                                 aTCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCa                                                                                                      >  1:356649/1‑150 (MQ=255)
                                                   ccGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATAGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                                                                    <  2:1188610/150‑1 (MQ=255)
                                                    cGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                                     >  1:1285314/1‑148 (MQ=255)
                                                    cGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                                     <  2:1285314/148‑1 (MQ=255)
                                                      gAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt                                                                                                 >  2:779560/1‑150 (MQ=255)
                                                       cggccATGGGTATCCAGATCGCCGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTACTGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                                                                <  2:1167093/147‑1 (MQ=255)
                                                        gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAACGTTGACGGCAAGGATTc                                                                                               >  1:1798285/1‑150 (MQ=255)
                                                        gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGAGCTGACGCCACGGCAAACAACAACGGCACAACCACCGGTGACGGCAAGGATTc                                                                                               >  1:1753836/1‑150 (MQ=255)
                                                           cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                                                            <  2:1921098/150‑1 (MQ=255)
                                                           cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                                                            >  2:1406655/1‑150 (MQ=255)
                                                              gggAATCCAGATCGACGGTGACAAGGCCAACGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                                                         >  1:114720/1‑150 (MQ=255)
                                                                gTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                                                       <  1:679255/150‑1 (MQ=255)
                                                                atATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                                                       <  2:1213748/149‑1 (MQ=255)
                                                                  aTCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTCCg                                                                                     >  1:2132461/1‑150 (MQ=255)
                                                                   tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                                                    >  2:2336119/1‑150 (MQ=255)
                                                                    ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                                                   >  2:953750/1‑150 (MQ=255)
                                                                     cAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCACGGATTCCACCGGTACGGaa                                                                                  >  1:1023043/1‑150 (MQ=255)
                                                                     cAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaa                                                                                  <  2:2598873/150‑1 (MQ=255)
                                                                     cAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACAACCGTTGACGGCAAGGATTCCACCGGTACGGaa                                                                                  >  1:1022898/1‑150 (MQ=255)
                                                                       gATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAt                                                                                <  2:802031/150‑1 (MQ=255)
                                                                       gATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAt                                                                                <  2:1356033/150‑1 (MQ=255)
                                                                       gATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAt                                                                                <  1:66789/150‑1 (MQ=255)
                                                                       gATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAt                                                                                <  1:66528/150‑1 (MQ=255)
                                                                        aTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATc                                                                               >  1:665470/1‑150 (MQ=255)
                                                                           gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAttgct                                                                            <  1:1740921/150‑5 (MQ=255)
                                                                           gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAttgct                                                                            <  2:2624132/150‑5 (MQ=255)
                                                                           gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAttgct                                                                            >  1:2624132/1‑146 (MQ=255)
                                                                           gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAttgct                                                                            >  2:1740921/1‑146 (MQ=255)
                                                                              ggTGACTTGTCCTTCGTCTATAACGAAGGCGAGTTCACCATAACCTACGCTGTCACCGTCACCGAGTCTAATGGTTATTACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGt                                                                         <  2:927543/150‑1 (MQ=14)
                                                                                         cttCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCGATATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaa                                                              <  1:844166/148‑1 (MQ=255)
                                                                                          ttggTCTATCACGGAGGCGAGAGCTCGATCTCCAACGGGGGCTCCGGGACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                                                             <  2:492639/147‑1 (MQ=17)
                                                                                          aTCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                                                             <  1:1406655/150‑1 (MQ=255)
                                                                                               cAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                                                        <  1:671779/150‑1 (MQ=255)
                                                                                                aaTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAtt                                                       <  1:1557039/150‑1 (MQ=255)
                                                                                                aaTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAtt                                                       <  1:1484364/150‑1 (MQ=255)
                                                                                                aaTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAtt                                                       <  1:1484375/150‑1 (MQ=255)
                                                                                                 aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc                                                      <  2:472066/150‑1 (MQ=255)
                                                                                                    aCGAAGGCGAGATCTCGATCTACAACGGTGGCTCCGGCACCCAGATTAATGGTGATTACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                                                   <  2:329232/150‑1 (MQ=17)
                                                                                                        aGGCGAGAGCACCATCAACAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACgg                                               <  1:963382/150‑1 (MQ=255)
                                                                                                          gCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCg                                             >  1:1657333/1‑150 (MQ=255)
                                                                                                               agCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTg                                        <  1:248275/150‑1 (MQ=255)
                                                                                                               agCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTg                                        <  1:1641076/150‑1 (MQ=255)
                                                                                                                 cACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAaccacc                                                                                                                <  2:160770/76‑1 (MQ=255)
                                                                                                                 cACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAaccacc                                                                                                                >  1:160770/1‑76 (MQ=255)
                                                                                                                 cACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAGGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGa                                      >  1:2060592/1‑150 (MQ=255)
                                                                                                                   ccAACACCAACCCTCTCACCGCCAAACAGATTAATGGTGATCACGCCACTGCAAACAACAACGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                                    <  2:90633/150‑1 (MQ=17)
                                                                                                                      tCACCAACTGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCa                                 <  2:2607496/150‑1 (MQ=255)
                                                                                                                      tCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCa                                 <  2:1629949/150‑1 (MQ=255)
                                                                                                                        aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgc                               <  1:1068384/150‑1 (MQ=255)
                                                                                                                        aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgc                               >  2:830396/1‑150 (MQ=255)
                                                                                                                          cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcgg                             >  2:239142/1‑150 (MQ=255)
                                                                                                                          cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcgg                             >  1:1642180/1‑150 (MQ=255)
                                                                                                                            aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACg                                                                                                                           >  1:1083890/1‑54 (MQ=255)
                                                                                                                            aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACg                                                                                                                           <  2:1083890/54‑1 (MQ=255)
                                                                                                                              ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcggtgat                         <  2:147070/150‑5 (MQ=255)
                                                                                                                                     ccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg                  >  1:1868915/1‑150 (MQ=255)
                                                                                                                                     ccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg                  >  1:488123/1‑150 (MQ=255)
                                                                                                                                     ccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg                  <  1:2290009/150‑1 (MQ=255)
                                                                                                                                        tgattgtgATTAATGGTGATTACGCCACTGAAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTa               <  2:1739053/143‑1 (MQ=17)
                                                                                                                                               gATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGatat        <  2:1907672/150‑1 (MQ=255)
                                                                                                                                                 ttAATGGTGATTACGCCACTGCAAACAACAACGGCAAAACCACCGTTGACGGAAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCa      <  2:1873392/150‑1 (MQ=255)
                                                                                                                                                     tGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg  >  1:1597882/1‑150 (MQ=255)
                                                                                                                                                     tGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg  >  1:1597875/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGG  >  NC_000913/1466196‑1466494

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