Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,019,735 G→T 11.5% P107T (CCA→ACA)  ompA ← outer membrane protein A (3a;II*;G;d)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,019,7350GT11.5% 253.1 / 3.5 91P107T (CCA→ACA) ompAouter membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base G (35/44);  new base T (10/2);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 1.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCAC  >  NC_000913/1019587‑1019878
                                                                                                                                                    |                                                                                                                                               
aTTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTgg                                                                                                                                                 >  1:370349/1‑149 (MQ=255)
    cAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCTGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTa                                                                                                                                              <  8:182989/148‑1 (MQ=255)
                cgTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCg                                                                                                                                 >  1:449141/1‑149 (MQ=255)
                     tCAACACCGCCAGCGAAGACCGGAGAATCGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGtt                                                                                                                                      <  1:131368/139‑1 (MQ=255)
                     tCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTTGGTAACCAAATTTAGCGGTCAg                                                                                                                            >  6:41981/1‑149 (MQ=255)
                     tCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGTGTAACCCAGTTTAGCGGTCAg                                                                                                                            >  4:353222/1‑149 (MQ=255)
                     tCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAg                                                                                                                            >  3:54749/1‑149 (MQ=255)
                       aaCACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGtt                                                                                                                          >  2:333812/1‑149 (MQ=255)
                       aaCACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGtt                                                                                                                          >  5:305095/1‑149 (MQ=255)
                       aaCACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGtt                                                                                                                          <  1:450552/149‑1 (MQ=255)
                       aaCACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGAATTGTTAAACAAGTTTAGCGGTCAGtt                                                                                                                          >  4:33942/1‑149 (MQ=255)
                         cacCGCCAGCGAAGACCTGAGAAACGCCCGTGTCCTGGTTTTTACCATAAACGTTGGATTTAGAGTCTGCACGCCATACAATGCCACCCAGACTAGTGTAGATCTCCAGGTAGTCAGTGATTTTATAACCAAGTTTAGCGGTCAGTTg                                                                                                                         >  6:359410/1‑148 (MQ=255)
                         cacCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCATGTCCTCAGTGAATTTGAAAACAAATTTAGCGGTCAGTTg                                                                                                                         >  6:3337/1‑148 (MQ=255)
                           ccGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAAc                                                                                                                      <  4:54749/149‑1 (MQ=255)
                            cGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGGTTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATATTGTAACCCAGTTTAGCGGTCAGTTGAACg                                                                                                                     >  4:365906/1‑149 (MQ=255)
                               cAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGccc                                                                                                                  >  6:70566/1‑149 (MQ=255)
                                aGCGATGACGGGAGAAACTCCGGTGTCGGGGTTTTTACCATAAACGATGGATTTAGTTGCTGATCGGCATACCATGCCACCCAGACGAGTGTAGATGTCAAGGTCGGCAGTTATTGGGTAACCCAGTTTAGCGGTCACTTCAACGCCCt                                                                                                                 <  5:401305/149‑1 (MQ=255)
                                  cGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTGGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCt                                                                                                                 <  5:434533/147‑1 (MQ=255)
                                     aGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGTGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCt                                                                                                            <  5:118369/149‑1 (MQ=255)
                                       aCCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCttt                                                                                                          <  1:333812/149‑1 (MQ=255)
                                       aCCGGAGAAACGCCGGTGTCGTGGTTGTTACCATAAACGTTGGATGTAGTGGCTGCCCGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGGCAGTTGAACGCCCTGAGCttt                                                                                                          <  5:372688/149‑1 (MQ=255)
                                         cGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCAAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGATAACCCAGTTTAGAGGTCAGTTGAACGCCCTGAGCTTTg                                                                                                         >  4:88691/1‑148 (MQ=255)
                                          ggAGAAACGCCCGTGTCGTGGTTTGTACCATAACCGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTa                                                                                                       <  6:296836/149‑1 (MQ=255)
                                               aaCGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACACAGACGAGTGTAGATGTCCATGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGGCAGTTGAACGCCCTGAGCTTTGTATGCa                                                                                                   <  6:277461/148‑1 (MQ=255)
                                                aCGCCGGTGTCGTGGTTTTTACCATATACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGaa                                                                                                                       <  5:218655/127‑1 (MQ=255)
                                                aCGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGTGTAAACCAGTTAAGCGGTCAGTTGaa                                                                                                                       >  6:218655/1‑127 (MQ=255)
                                                aCGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATCTCCATGTACTCAGTAAAATTATAAAAAATTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCAcc                                                                                                 >  6:41222/1‑149 (MQ=255)
                                                    caccggcGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGttt                                                                                             <  5:123394/143‑1 (MQ=255)
                                                    cGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGttt                                                                                             <  2:370349/149‑1 (MQ=255)
                                                         tCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAAc                                                                                        <  7:232710/149‑1 (MQ=255)
                                                           gTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGc                                                                                      >  8:227228/1‑149 (MQ=255)
                                                             ggTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTg                                                                                    >  6:99450/1‑149 (MQ=255)
                                                                  ttACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCtt                                                                                <  5:41981/148‑1 (MQ=255)
                                                                    aCCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTg                                                                              >  8:192639/1‑148 (MQ=255)
                                                                    aCCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTg                                                                              <  7:227228/148‑1 (MQ=255)
                                                                     ccATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGt                                                                             <  5:41222/148‑1 (MQ=255)
                                                                     ccATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTa                                                                            >  1:60280/1‑149 (MQ=255)
                                                                     ccATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTa                                                                            <  5:70566/149‑1 (MQ=255)
                                                                     ccATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTa                                                                            <  3:353222/149‑1 (MQ=255)
                                                                     ccATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTa                                                                            <  1:423234/149‑1 (MQ=255)
                                                                       aTAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACg                                                                          <  5:99450/149‑1 (MQ=255)
                                                                       aTAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACg                                                                          >  3:399218/1‑149 (MQ=255)
                                                                       aTAAACGTTGGATTTAGGGGCTGCACGCCATACCATGCCACCAAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGCACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACg                                                                          <  3:88691/149‑1 (MQ=255)
                                                                          aaCGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGc                                                                        <  7:192639/148‑1 (MQ=255)
                                                                           aCGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATCTCCATGTACTCCGTAAAATTAAAAAAAAAATTAGCTGTCAGTTGAACGCACTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCa                                                                       >  4:385963/1‑148 (MQ=255)
                                                                             gTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCtt                                                                                >  8:410433/1‑137 (MQ=255)
                                                                             gTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCtt                                                                                <  7:410433/137‑1 (MQ=255)
                                                                              ttGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACg                                                                   >  1:113621/1‑149 (MQ=255)
                                                                              ttGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTAGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCCTTTTCAACGCTGCCTTTGTACGGCATAc                                                                    >  1:80150/1‑148 (MQ=255)
                                                                                  aTTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCt                                                               <  2:449141/149‑1 (MQ=255)
                                                                                   tttAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCt                                                               >  5:72511/1‑148 (MQ=255)
                                                                                     tAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTaa                                                             <  6:72511/148‑1 (MQ=255)
                                                                                          tCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTc                                                       <  5:359410/149‑1 (MQ=255)
                                                                                          tCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTc                                                       <  1:359735/149‑1 (MQ=255)
                                                                                            tGCACGCCATACCATGCCACCCAGACGAGTGTAGTTGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCTCCGTTTTCAACGCTGCCTTTGTACGGCATACGACCt                                                               <  5:334971/139‑1 (MQ=255)
                                                                                            tGCACGCCATACCATGCCACCCAGACGAGTGTAGCTGTCCAGGTCGGCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCGGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCg                                                      >  1:270675/1‑148 (MQ=255)
                                                                                            tGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTa                                                                                                       <  8:446624/99‑1 (MQ=255)
                                                                                            tGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTa                                                                                                       >  7:440537/1‑99 (MQ=255)
                                                                                            tGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTa                                                                                                       >  7:446624/1‑99 (MQ=255)
                                                                                            tGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTa                                                                                                       <  8:440537/99‑1 (MQ=255)
                                                                                            tGCACGCCATACCATGCCACCCAGACGAGTGTAGATCTCAATGTACTCAGTAAAATTAAAAAACAATTTATCGGTCAGTTGAACGCCCTGAGATTTGTATGCACCGTTTTCAACGCAGCCTTAGTACGGCATACGACCt                                                               >  6:334971/1‑139 (MQ=255)
                                                                                              cACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTACTCAGTAAAATTAAAAAAAAATTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTaa                                                   >  4:64620/1‑149 (MQ=255)
                                                                                                    aTACCATGCCACCCAGACGAGTGTAGATCTCCATGTCCTCAGTTAATTTTTAACAAAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCAtt                                             >  2:149574/1‑149 (MQ=255)
                                                                                                     tACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCAttt                                            <  1:429675/149‑1 (MQ=255)
                                                                                                             cACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCa                                    <  2:60280/149‑1 (MQ=255)
                                                                                                              aCCCAGACGAGTGTAGATGTCCAGGGCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCaa                                   <  3:33942/149‑1 (MQ=255)
                                                                                                                   gACGAGTGTAGATGGCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGGCAGTTGACAGCCCAGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATAc                              >  1:148082/1‑149 (MQ=255)
                                                                                                                      gAGTGTAGATGTCCAGGGCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACggg                           <  6:305095/149‑1 (MQ=255)
                                                                                                                           tAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCaa                                                                                         >  6:248472/1‑82 (MQ=255)
                                                                                                                           tAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCaa                                                                                         <  5:248472/82‑1 (MQ=255)
                                                                                                                            aGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAAcc                     >  7:418068/1‑149 (MQ=255)
                                                                                                                             gATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCaaa                                        <  4:36840/129‑1 (MQ=255)
                                                                                                                             gATGTCAAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCaaa                                        >  3:36840/1‑129 (MQ=255)
                                                                                                                              aTGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTc                                           >  3:24080/1‑125 (MQ=255)
                                                                                                                              aTGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTc                                           <  4:24080/125‑1 (MQ=255)
                                                                                                                              aTGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGAGGTCAGTTGAACGCCCTGATCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGACGTAACCCATTTCAAAGCCAACATACGGGTAAACCTg                   >  6:303535/1‑149 (MQ=255)
                                                                                                                                gTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTgg                  >  2:410567/1‑148 (MQ=255)
                                                                                                                                   cAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGcc                                     <  3:307350/126‑1 (MQ=255)
                                                                                                                                   cAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGcc                                     >  4:307350/1‑126 (MQ=255)
                                                                                                                                   cAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAac              >  6:444963/1‑149 (MQ=255)
                                                                                                                                    aGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCaaa                                        >  7:7557/1‑122 (MQ=255)
                                                                                                                                    aGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCaaa                                        <  8:7557/122‑1 (MQ=255)
                                                                                                                                     ggTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAacca            <  1:410567/149‑1 (MQ=255)
                                                                                                                                     ggTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAacca            >  7:134588/1‑149 (MQ=255)
                                                                                                                                      gtcgtcTGTGTTTGTGTAACCCAGTTTAGCTGTCAGTTGAACGCCCTGAGCTTTGTATGGACCTTTTTCAACGCTGCCTTTGTACGGCATACGACCTAAc                                                            <  2:18210/100‑1 (MQ=255)
                                                                                                                                      gtcgtcAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAAc                                                            >  1:18210/1‑100 (MQ=255)
                                                                                                                                      gtcgtcAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAaccac           >  1:89211/1‑149 (MQ=255)
                                                                                                                                            aGTTATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCaaaa      <  3:365906/148‑1 (MQ=255)
                                                                                                                                            aGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAg     >  6:304391/1‑149 (MQ=255)
                                                                                                                                                aTTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCTTTTTCATCGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCAc  <  2:80150/148‑1 (MQ=255)
                                                                                                                                                aTTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCAc  <  4:399218/148‑1 (MQ=255)
                                                                                                                                                 ttGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCAc  >  5:404249/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCAC  >  NC_000913/1019587‑1019878

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: