Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 962,939 C→A N315K (AAC→AAA rpsA → 30S ribosomal subunit protein S1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913962,9390CA98.9% 344.5 / ‑4.4 87N315K (AAC→AAArpsA30S ribosomal subunit protein S1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (41/45);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGG  >  NC_000913/962794‑963085
                                                                                                                                                 |                                                                                                                                                  
gggTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATcc                                                                                                                                                >  1:124390/1‑150 (MQ=255)
 ggTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCAc                                                                                                                                              >  2:227594/1‑151 (MQ=255)
 ggTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCGGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCAc                                                                                                                                              <  2:336874/151‑1 (MQ=255)
    aGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCAccc                                                                                                                                            <  2:741635/150‑1 (MQ=255)
    aGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTCCCGAAATGGACTGGACCAACAAAAAAATCCACCCg                                                                                                                                           >  2:440680/1‑151 (MQ=255)
        aTCGCTAAACGTTATCCTGAAGGTACCAAACTGACTGGTGGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCAccc                                                                                                                                            <  1:1027918/146‑1 (MQ=255)
        aTCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTcc                                                                                                                                        <  1:997794/150‑1 (MQ=255)
         tCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCa                                                                                                                                       >  1:1141931/1‑150 (MQ=255)
          cGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCaa                                                                                                                                      >  1:1197103/1‑150 (MQ=255)
           gCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCaa                                                                                                                                      >  2:165890/1‑149 (MQ=255)
           gCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAg                                                                                                                                    <  2:469540/151‑1 (MQ=255)
            cTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAgt                                                                                                                                   >  2:143305/1‑151 (MQ=255)
             tAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAgtt                                                                                                                                  >  1:45217/1‑151 (MQ=255)
              aaaCGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAgttg                                                                                                                                 <  2:877719/151‑1 (MQ=255)
               aaCGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAgttgt                                                                                                                                <  1:920073/151‑1 (MQ=255)
               aaCGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAgttgt                                                                                                                                >  2:375918/1‑151 (MQ=255)
                aCGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAgttgtt                                                                                                                               >  2:422904/1‑151 (MQ=255)
                  gTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTaa                                                                                                                             <  1:476641/151‑1 (MQ=255)
                  gTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTa                                                                                                                              >  2:933265/1‑150 (MQ=255)
                  gTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTa                                                                                                                              >  1:1049124/1‑150 (MQ=255)
                     aTCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGt                                                                                                                          >  2:616167/1‑151 (MQ=255)
                       ccGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTg                                                                                                                        >  1:509010/1‑151 (MQ=255)
                          gAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCg                                                                                                                     <  2:99129/151‑1 (MQ=255)
                            aGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGAt                                                                                                                   >  1:1034970/1‑151 (MQ=255)
                            aGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGAt                                                                                                                   <  2:846098/151‑1 (MQ=255)
                              gTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGt                                                                                                                 <  2:1051447/151‑1 (MQ=255)
                                  cAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTg                                                                                                             <  1:1061447/151‑1 (MQ=255)
                                  cAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTg                                                                                                             <  2:49715/151‑1 (MQ=255)
                                   aaaCTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTgg                                                                                                            <  1:522062/151‑1 (MQ=255)
                                   aaaCTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTgg                                                                                                            >  2:1239157/1‑151 (MQ=255)
                                   aaaCTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTgg                                                                                                            <  2:1192784/151‑1 (MQ=255)
                                     actgactgGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGaa                                                                                                          <  1:339606/151‑1 (MQ=255)
                                        gactgGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAg                                                                                                         <  2:19848/149‑1 (MQ=255)
                                            ggTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATg                                                                                                    <  2:1108713/150‑1 (MQ=255)
                                              tCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGtt                                                                                                 >  1:622949/1‑151 (MQ=255)
                                               cgcgTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTc                                                                                                >  1:912998/1‑151 (MQ=255)
                                               cgcgTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTc                                                                                                <  1:784684/151‑1 (MQ=255)
                                                  gTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTgg                                                                                             <  2:22264/151‑1 (MQ=255)
                                                  gTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTgg                                                                                             <  2:433164/151‑1 (MQ=255)
                                                    gACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGat                                                                                           >  1:730369/1‑151 (MQ=255)
                                                    gACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGat                                                                                           >  1:48752/1‑151 (MQ=255)
                                                    gACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGa                                                                                            <  2:428643/150‑1 (MQ=255)
                                                        aaCCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATcg                                                                                       <  2:296922/151‑1 (MQ=255)
                                                         aCCTGACCGACTACGGCGGCGTCGTTGAAATCGAAGAAGGCGGTGAAGGCCTGGTACGCGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTGAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATcg                                                                                       <  2:641824/150‑1 (MQ=255)
                                                              aCCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACgaag                                                                                 <  1:674941/151‑1 (MQ=255)
                                                                cGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACgaaga                                                                                >  1:422827/1‑150 (MQ=255)
                                                                 gACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAc                                                                              >  2:2543/1‑151 (MQ=255)
                                                                  aCTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcg                                                                             <  1:1068471/151‑1 (MQ=255)
                                                                  aCTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcg                                                                             <  1:354352/151‑1 (MQ=255)
                                                                    tACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtc                                                                           >  1:1132479/1‑151 (MQ=255)
                                                                     actgctgctTCGTTGAAATCGAAGAAGCCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATTTCGACGAAGAAcgt                                                                            <  2:640723/146‑1 (MQ=255)
                                                                     aCGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtc                                                                           <  2:508995/150‑1 (MQ=255)
                                                                       ggCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtcgtc                                                                        >  2:479150/1‑151 (MQ=255)
                                                                       ggCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtcgtc                                                                        <  1:670169/151‑1 (MQ=255)
                                                                         ctgctTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtcgtcgt                                                                      >  1:56695/1‑151 (MQ=255)
                                                                                gTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTccc                                                               >  2:1082060/1‑151 (MQ=255)
                                                                                  tGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCt                                                              <  2:360936/150‑1 (MQ=255)
                                                                                  tGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCt                                                              <  1:774297/150‑1 (MQ=255)
                                                                                   gAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTgg                                                            >  1:623322/1‑151 (MQ=255)
                                                                                   gAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTg                                                             <  2:678920/150‑1 (MQ=255)
                                                                                    aaaTCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTggg                                                           <  1:631190/151‑1 (MQ=255)
                                                                                    aaaTCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTggg                                                           >  2:460747/1‑151 (MQ=255)
                                                                                    aaaTCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTAGTCGTATCTCCCTggg                                                           >  1:638340/1‑151 (MQ=255)
                                                                                        cGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTg                                                       <  2:6018/151‑1 (MQ=255)
                                                                                         gaagaaGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTg                                                       >  2:1168922/1‑150 (MQ=255)
                                                                                         gaagaaGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGa                                                      <  2:1181001/151‑1 (MQ=255)
                                                                                         gaagaaGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGa                                                      <  1:152744/151‑1 (MQ=255)
                                                                                              agGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAg                                                 <  1:616184/151‑1 (MQ=255)
                                                                                                 cGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGc                                              >  2:536571/1‑151 (MQ=255)
                                                                                                      aaGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGc                                         >  1:589155/1‑151 (MQ=255)
                                                                                                      aaGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGc                                         >  1:912110/1‑151 (MQ=255)
                                                                                                         gCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTa                                       >  2:1224816/1‑150 (MQ=255)
                                                                                                                  acGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGgca                             <  2:405173/151‑1 (MQ=255)
                                                                                                                   cGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGgc                              >  2:806551/1‑149 (MQ=255)
                                                                                                                   cGTTTCCGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGgc                              >  1:1021768/1‑149 (MQ=255)
                                                                                                                        ccGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGtt                       <  2:1141878/151‑1 (MQ=255)
                                                                                                                         cGAAATGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTc                      <  1:242886/151‑1 (MQ=255)
                                                                                                                           aaaTGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTcgc                    <  2:873735/151‑1 (MQ=255)
                                                                                                                           aaaTGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTcgc                    <  2:124386/151‑1 (MQ=255)
                                                                                                                             aTGGACTGGACCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCgg                  <  1:949506/151‑1 (MQ=255)
                                                                                                                               ggactggacCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGa                 >  1:512644/1‑150 (MQ=255)
                                                                                                                                   tggacCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGTGGAAAccc            >  1:194448/1‑151 (MQ=255)
                                                                                                                                    ggacCAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAAccc            >  2:239911/1‑150 (MQ=255)
                                                                                                                                       ccAACAAAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCacaa        <  2:1000689/151‑1 (MQ=255)
                                                                                                                                        caacaaAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCGGTTCGCGGAAACCCacaac       >  1:929832/1‑151 (MQ=255)
                                                                                                                                           caaAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCacaaca      >  1:244005/1‑149 (MQ=255)
                                                                                                                                           caaAAAAATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCacaaca      >  2:807503/1‑149 (MQ=255)
                                                                                                                                             aaaaaaTCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAggg  <  1:1137726/151‑1 (MQ=255)
                                                                                                                                             aaaaaaTCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAg    >  2:1143389/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGG  >  NC_000913/962794‑963085

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