Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,472,951 C→T 5.0% pseudogene (2435/3497 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,472,9510CT 5.0% 257.4 / 2.3 79pseudogene (2435/3497 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (34/41);  new base T (2/2);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

TACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGTCGGCACGACTGCTGAACTGAACAGTGCTCTGAAAC  >  NC_000913/1472807‑1473096
                                                                                                                                                |                                                                                                                                                 
tACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCa                                                                                                                                                  <  3:125111/146‑1 (MQ=255)
       ggTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCaaaa                                                                                                                                         <  2:467731/148‑1 (MQ=255)
             gCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCt                                                                                                                                   <  2:220347/148‑1 (MQ=255)
                 acacCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCaaaa                                                                                                                                         <  4:228215/138‑1 (MQ=255)
                 acacCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCaaaa                                                                                                                                         >  3:228215/1‑138 (MQ=255)
                      acTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCACCGGTAACGTTGACGTCACCATGACCAAAAATGCCTACACCGAtg                                                                                                                          >  3:108747/1‑148 (MQ=255)
                       cTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGAtgt                                                                                                                         <  1:298043/148‑1 (MQ=255)
                       cTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGAtgt                                                                                                                         <  7:40536/148‑1 (MQ=255)
                       atgtgAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTTTGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGAtgt                                                                                                                         <  5:340492/147‑1 (MQ=255)
                         gtgAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACAcc                                                                                                                              <  8:431883/141‑1 (MQ=255)
                         gtgAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACAcc                                                                                                                              >  7:431883/1‑141 (MQ=255)
                         gtgAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTgg                                                                                                                       >  1:79114/1‑148 (MQ=255)
                         gtgAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTgg                                                                                                                       >  3:227392/1‑148 (MQ=255)
                           gAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACACCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGc                                                                                                                      >  3:112705/1‑147 (MQ=255)
                             gTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAAc                                                                                                                   <  1:275814/148‑1 (MQ=255)
                              ttttGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCGCCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCaa                                                                                                                    >  6:261120/1‑146 (MQ=255)
                              ttttGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCaa                                                                                                                    >  3:293683/1‑146 (MQ=255)
                                ttGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACa                                                                                                                  >  1:369454/1‑146 (MQ=255)
                                 tGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGa                                                                                                                <  5:51359/147‑1 (MQ=255)
                                    tAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGcc                                                                                                            <  4:227392/148‑1 (MQ=255)
                                    tAACGTAGTGGAAGGTAGCAACCTGACAGACGCTGACTCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGc                                                                                                             >  5:157197/1‑147 (MQ=255)
                                            tGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGa                                                                                                     <  2:79114/147‑1 (MQ=255)
                                               aaGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGa                                                                                                 <  4:293683/148‑1 (MQ=255)
                                                 ggTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACa                                                                                               >  6:25326/1‑148 (MQ=255)
                                                   tAGCAACCTGGCCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCa                                                                                                                                            <  2:146054/101‑1 (MQ=255)
                                                   tAGCAACCTGGCCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCa                                                                                                                                            >  1:146054/1‑101 (MQ=255)
                                                     gCACCCTGACCGCCGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAACTGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATcg                                                                                            <  5:129321/147‑1 (MQ=255)
                                                     gCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATcg                                                                                            <  4:29011/147‑1 (MQ=255)
                                                             aCCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCa                                                                                    <  3:47968/147‑1 (MQ=255)
                                                             aCCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCAc                                                                                   >  4:190885/1‑148 (MQ=255)
                                                              ccGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGTCTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACt                                                                                  <  5:24653/148‑1 (MQ=255)
                                                                gACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTgg                                                                                <  4:430993/148‑1 (MQ=255)
                                                                gACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTgg                                                                                <  6:398685/148‑1 (MQ=255)
                                                                 aCGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACTATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGa                                                                               >  2:242250/1‑148 (MQ=255)
                                                                      gACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGc                                                                            >  4:224335/1‑146 (MQ=255)
                                                                       aCGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCGCCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCgg                                                                          <  5:261120/147‑1 (MQ=255)
                                                                       aCGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCg                                                                           <  3:224335/146‑1 (MQ=255)
                                                                           cATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTac                                                                     <  2:395645/148‑1 (MQ=255)
                                                                           cATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTac                                                                     <  2:49834/148‑1 (MQ=255)
                                                                           cATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTa                                                                      >  5:22632/1‑147 (MQ=255)
                                                                            aTCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTaca                                                                    <  4:112705/148‑1 (MQ=255)
                                                                              cACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTacac                                                                   >  8:328202/1‑147 (MQ=255)
                                                                                 cTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTCACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGCTGCCTCGGTGAATGCCATCGCGAAAGCACTGGATGACGGTTACCCcaa                                                                >  3:47961/1‑147 (MQ=255)
                                                                                 cTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCTAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACcaa                                                                >  3:54047/1‑147 (MQ=255)
                                                                                  tCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACcaac                                                               >  5:435694/1‑147 (MQ=255)
                                                                                   cAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATg                                                                                                                                      <  2:125050/75‑1 (MQ=255)
                                                                                   cAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATg                                                                                                                                      >  1:125050/1‑75 (MQ=255)
                                                                                      cGTCCTTGGTATTGAATGACAAATGCAGCACCGATGCCAGCGTTAAAGTTGACGTCACCATTAGCAAAAATGCCTACACCGATTTGGCAACAGATTCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGa                                                          <  2:318620/148‑1 (MQ=255)
                                                                                            tGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACTATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGtt                                                    <  1:242250/148‑1 (MQ=255)
                                                                                              gTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGttt                                                   <  5:25326/147‑1 (MQ=255)
                                                                                              gTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTa                                                  <  7:124810/148‑1 (MQ=255)
                                                                                                  ggACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCaaaa                                                                                                                                         <  1:439085/57‑1 (MQ=255)
                                                                                                  ggACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCaaaa                                                                                                                                         >  2:439085/1‑57 (MQ=255)
                                                                                                  ggACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGttt                                                   >  5:179324/1‑143 (MQ=255)
                                                                                                  ggACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGttt                                                   <  6:179324/143‑1 (MQ=255)
                                                                                                    aCTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATg                                                                                                              >  7:300992/1‑82 (MQ=255)
                                                                                                    aCTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATg                                                                                                              <  8:300992/82‑1 (MQ=255)
                                                                                                    aCTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGcc                                            >  6:374375/1‑148 (MQ=255)
                                                                                                    aCTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGcc                                            >  8:5808/1‑148 (MQ=255)
                                                                                                     cTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCt                                           >  2:220511/1‑148 (MQ=255)
                                                                                                        ccAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGa                                         <  3:190885/147‑1 (MQ=255)
                                                                                                         cAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGAt                                                                                                                           >  6:196270/1‑64 (MQ=255)
                                                                                                         cAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGAt                                                                                                                           <  5:196270/64‑1 (MQ=255)
                                                                                                         cAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAAc                                       >  1:218042/1‑148 (MQ=255)
                                                                                                         cAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAAc                                       >  4:433739/1‑148 (MQ=255)
                                                                                                           aaGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACTATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGCATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGt                                     >  3:360532/1‑148 (MQ=255)
                                                                                                           aaGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGt                                     >  3:159221/1‑148 (MQ=255)
                                                                                                                 gcaCCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGCATGCCATCGCGACGCACTGGATGCGGGGTACACCAACACCGAACTGTTTACCAGCCTGAACGTCGGCACg                              >  8:166431/1‑148 (MQ=255)
                                                                                                                  caCCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTAGACCAACAACGAACTGTTTACCAGCCTGAACGTCGGCACg                              <  2:466239/148‑1 (MQ=255)
                                                                                                                     cGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGTCGGCACGAct                           >  1:64159/1‑148 (MQ=255)
                                                                                                                      gATGCCAGCGGTAACGTTGACGTCACCATTAGCAAAAATGCATACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGTCGGCACGActg                          <  4:37174/148‑1 (MQ=255)
                                                                                                                           cAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGTCGGCACGACTGctgaa                     <  7:418019/148‑1 (MQ=255)
                                                                                                                                gTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGTCGGCACGACTGctgaactgaa                >  8:165092/1‑148 (MQ=255)
                                                                                                                                gTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGTCGGCACGACTGctgaactgaa                <  2:218042/148‑1 (MQ=255)
                                                                                                                                  aaCGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGTCGGCACGACTGctgaactgaa                >  2:160225/1‑146 (MQ=255)
                                                                                                                                          cGTCACTATGAGCTAAAATGCCTTCACCGATGTTGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGTCTGCACGACTGCTGAACTGAACAGTGCTCTg      <  4:360532/148‑1 (MQ=255)
                                                                                                                                              aCCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGTCGGCACGACTGCTGAACTGAACAGTGCTCTGAAAc  <  2:64159/148‑1 (MQ=255)
                                                                                                                                              aCCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGTCGGCACGACTGCTGAACTGAACAGTGCTCTGAAAc  <  5:206496/148‑1 (MQ=255)
                                                                                                                                              aCCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTGTACCAGCCTGAACGTCTGCACGACTGCTGAACTGAACAGTGCTCTGAAAc  <  6:157197/148‑1 (MQ=255)
                                                                                                                                               ccATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGTCGGCACGACTGCTGAACTGAACAGTGCTCTGaaa   >  8:304498/1‑146 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGGGTTACACCAACAACGAACTGTTTACCAGCCTGAACGTCGGCACGACTGCTGAACTGAACAGTGCTCTGAAAC  >  NC_000913/1472807‑1473096

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: