Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 962,923 A→G D310G (GAC→GGC)  rpsA → 30S ribosomal subunit protein S1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913962,9230AG96.5% 315.1 / ‑1.5 86D310G (GAC→GGC) rpsA30S ribosomal subunit protein S1
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base G (41/42);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGG  >  NC_000913/962776‑963069
                                                                                                                                                   |                                                                                                                                                  
aGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAAtgggctg                                                                                                                                                 >  1:1320497/1‑151 (MQ=255)
 gCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAAtgggctgg                                                                                                                                                >  2:636242/1‑151 (MQ=255)
    gggCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGAcc                                                                                                                                             >  2:632534/1‑151 (MQ=255)
      gCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACcaa                                                                                                                                           <  2:1383461/151‑1 (MQ=255)
         aaGATCCGTGGGTAGGTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACcaaca                                                                                                                                         <  2:2799/150‑1 (MQ=255)
               cGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCGGGACCAACAAAAACa                                                                                                                                   <  1:598407/150‑1 (MQ=255)
                gTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATc                                                                                                                                 <  2:61050/151‑1 (MQ=255)
                    gTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCAcc                                                                                                                             >  2:732929/1‑151 (MQ=255)
                      aGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCg                                                                                                                           >  1:1033312/1‑151 (MQ=255)
                       gCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGt                                                                                                                          <  1:715613/151‑1 (MQ=255)
                       gCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGt                                                                                                                          <  1:657928/151‑1 (MQ=255)
                            cGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCaa                                                                                                                      >  2:800428/1‑150 (MQ=255)
                                         ccGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGtt                                                                                                         <  2:866960/150‑1 (MQ=255)
                                         ccGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTg                                                                                                        <  1:250714/151‑1 (MQ=255)
                                          cGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTgg                                                                                                       <  1:625462/151‑1 (MQ=255)
                                            gAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCg                                                                                                     <  1:271397/151‑1 (MQ=255)
                                              atgTACCAAAGTGACTGGTCGCGTGACAAAGCCGACCGACTACGGCTGCTTCGTTAAAATAGAAGCAGGCGTTGCAGGCCTGGTACATGATTCCGCAATGCGCTGGGCCAACACAAACATACACCCGTCCACAGCTGTTAAc                                                                                                            >  1:1304119/3‑142 (MQ=255)
                                              aGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGa                                                                                                    >  1:891835/1‑150 (MQ=255)
                                               ggTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATg                                                                                                  >  1:444392/1‑151 (MQ=255)
                                                  aCCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTa                                                                                                <  2:157010/150‑1 (MQ=255)
                                                    cAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTg                                                                                             <  2:675267/151‑1 (MQ=255)
                                                       actgactgGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGaa                                                                                          <  1:1288984/151‑1 (MQ=255)
                                                           actgGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGACATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAgct                                                                                       >  1:510243/1‑148 (MQ=255)
                                                             tgGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTAt                                                                                     >  1:785280/1‑150 (MQ=255)
                                                             tgGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATg                                                                                    >  2:1196194/1‑151 (MQ=255)
                                                              ggTCGCGTTACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATg                                                                                    <  1:1055863/150‑1 (MQ=255)
                                                              ggTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTAt                                                                                     <  2:530757/149‑1 (MQ=255)
                                                              ggTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATgg                                                                                   >  2:1222439/1‑151 (MQ=255)
                                                              ggTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATgg                                                                                   >  1:1068316/1‑151 (MQ=255)
                                                              ggTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATgg                                                                                   <  2:195055/151‑1 (MQ=255)
                                                              ggTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATg                                                                                    <  2:508991/150‑1 (MQ=255)
                                                                tCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGtt                                                                                 >  2:1293046/1‑151 (MQ=255)
                                                                 cgcgTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGtt                                                                                 <  1:250851/150‑1 (MQ=255)
                                                                 cgcgTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTc                                                                                >  1:625270/1‑151 (MQ=255)
                                                                 cgcgTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTc                                                                                <  1:641660/151‑1 (MQ=255)
                                                                   cgTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTg                                                                              <  1:423726/151‑1 (MQ=255)
                                                                    gTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTg                                                                              <  2:447572/150‑1 (MQ=255)
                                                                       aCCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGata                                                                          <  2:1329985/151‑1 (MQ=255)
                                                                          aaCCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATcg                                                                       <  2:815485/151‑1 (MQ=255)
                                                                          aaCCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATcg                                                                       <  2:1297215/151‑1 (MQ=255)
                                                                           aCCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATcg                                                                       >  2:479484/1‑150 (MQ=255)
                                                                            ccTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATcga                                                                      >  1:815178/1‑150 (MQ=255)
                                                                               gACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATcgacga                                                                   >  2:626222/1‑150 (MQ=255)
                                                                               gACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATcgacga                                                                   <  1:512224/150‑1 (MQ=255)
                                                                               gACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATcgacg                                                                    <  1:989061/149‑1 (MQ=255)
                                                                                 ccGACTACGGCTGCTTCGTTGAAATCGAAGCAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACgaaga                                                                <  2:784372/151‑1 (MQ=255)
                                                                                 ccGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACgaag                                                                 <  1:366536/150‑1 (MQ=255)
                                                                                    aCTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcg                                                             <  2:1186281/151‑1 (MQ=255)
                                                                                       aCGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtc                                                           <  2:135912/150‑1 (MQ=255)
                                                                                         ggCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtcgtc                                                        <  2:643366/151‑1 (MQ=255)
                                                                                          gctgctTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtcgtcg                                                       >  2:140/1‑151 (MQ=255)
                                                                                           ctgctTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtcgtcgt                                                      >  1:319516/1‑151 (MQ=255)
                                                                                             gctTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAAcgtcgtcg                                                       >  2:168231/1‑148 (MQ=255)
                                                                                              ctTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTAtc                                                   >  2:798206/1‑151 (MQ=255)
                                                                                              ctTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTAtc                                                   >  2:319231/1‑151 (MQ=255)
                                                                                                 cGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTAtctc                                                 <  1:1318483/150‑1 (MQ=255)
                                                                                                  gtggAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCGAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTccc                                               <  1:125914/148‑1 (MQ=255)
                                                                                                  gTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTccc                                               >  2:405615/1‑151 (MQ=255)
                                                                                                    tGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCt                                              <  1:94777/150‑1 (MQ=255)
                                                                                                     gAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTgg                                            <  2:746306/151‑1 (MQ=255)
                                                                                                       aaTCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGt                                          >  1:93752/1‑151 (MQ=255)
                                                                                                        aTCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTc                                         >  1:1067223/1‑151 (MQ=255)
                                                                                                              gaagGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACTTCGTCGTATCTCCCTGGGTCTGaaa                                    <  1:462363/150‑1 (MQ=255)
                                                                                                                   cGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGc                              <  2:37641/151‑1 (MQ=255)
                                                                                                                   cGTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGc                              >  1:696494/1‑151 (MQ=255)
                                                                                                                    gTTGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGc                              >  2:246625/1‑150 (MQ=255)
                                                                                                                      tGAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCaa                            >  2:65436/1‑150 (MQ=255)
                                                                                                                       gAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCaa                            >  1:979910/1‑149 (MQ=255)
                                                                                                                       gAAGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAg                          <  2:511353/151‑1 (MQ=255)
                                                                                                                        aaGGCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGc                         >  1:1032653/1‑151 (MQ=255)
                                                                                                                           gCCTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTa                       >  1:1374125/1‑150 (MQ=255)
                                                                                                                             cTGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAAcc                    >  2:1221328/1‑151 (MQ=255)
                                                                                                                              tGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAAccc                   >  2:662410/1‑151 (MQ=255)
                                                                                                                              tGGTACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAAccc                   >  2:770162/1‑151 (MQ=255)
                                                                                                                                ggacgcGTTTCGGAAATTGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAGCGTGGGCGATGTAGTGGAAAGTATGGTGCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCGAACACGt                 <  1:798140/146‑1 (MQ=255)
                                                                                                                                 tACACGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTg                >  1:529120/1‑151 (MQ=255)
                                                                                                                                  acacGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTgg               <  2:220275/151‑1 (MQ=255)
                                                                                                                                  acacGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTg                <  1:814148/150‑1 (MQ=255)
                                                                                                                                   cacGTTTCTGAAATGGACTGGACCAAAAAAAACATACCCCCGTGCAAAGTTGTGAACGTTGGCGAGGTAGTGGAAGTTATGGTTTTGGATAGAGACAAAGAACGTCGTCGTATCTCCCTAGGTCTGAAACGGTGCAAAGCAAACTCGTgc               <  1:57626/150‑2 (MQ=255)
                                                                                                                                    acGTTTCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGgc              >  1:1018704/1‑150 (MQ=255)
                                                                                                                                       tttCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGgcagc           >  1:511962/1‑150 (MQ=255)
                                                                                                                                        ttCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGgcagca          <  1:658614/150‑1 (MQ=255)
                                                                                                                                         tCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGCAACAGTGCAAAGCTAACCCGTGGCAGCAGt        >  2:790067/1‑151 (MQ=255)
                                                                                                                                         tCCGAAATGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGcagcag         <  1:349348/150‑1 (MQ=255)
                                                                                                                                             aaaTGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGtt       >  1:152371/1‑148 (MQ=255)
                                                                                                                                             aaaTGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTcgc    >  1:1380176/1‑151 (MQ=255)
                                                                                                                                             aaaTGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTcgc    <  2:407591/151‑1 (MQ=255)
                                                                                                                                               aTGGGCTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCgg  >  2:924018/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTGACTGGTCGCGTGACCAACCTGACCGACTACGGCTGCTTCGTTGAAATCGAAGAAGGCGTTGAAGGCCTGGTACACGTTTCCGAAATGGACTGGACCAACAAAAACATCCACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGG  >  NC_000913/962776‑963069

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