Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,608,607 C→G 8.4% I141M (ATC→ATG rimI → ribosomal‑protein‑S18‑alanine N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,608,6070CG 8.4% 234.2 / 10.7 71I141M (ATC→ATGrimIribosomal‑protein‑S18‑alanine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (25/40);  new base G (3/3);  total (28/43)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

GGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCGCTGAA  >  NC_000913/4608460‑4608756
                                                                                                                                                   |                                                                                                                                                     
cgAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatcat                                                                                                                                                     >  2:793048/2‑150 (MQ=255)
 gAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATg                                                                                                                                                    <  2:747396/150‑1 (MQ=255)
   aaaaCGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGc                                                                                                                                                  <  1:268407/150‑1 (MQ=255)
    aaaCGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCg                                                                                                                                                 <  2:1093518/150‑1 (MQ=255)
     aaCGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGt                                                                                                                                                >  1:1396920/1‑150 (MQ=255)
       cgcgGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtct                                                                                                                                              <  1:55078/150‑3 (MQ=255)
          ggcgTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                                                                           <  1:891738/150‑6 (MQ=255)
           gcgTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcgc                                                                                                                                          <  1:584155/150‑7 (MQ=255)
             gTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATc                                                                                                                                        <  2:307268/150‑1 (MQ=255)
               ggCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAg                                                                                                                                      >  2:1250430/1‑150 (MQ=255)
               ggCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAg                                                                                                                                      >  2:162494/1‑150 (MQ=255)
                 cGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTa                                                                                                                                    <  2:1396920/150‑1 (MQ=255)
                      cTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTaa                                                                                                                               <  2:14540/150‑1 (MQ=255)
                          ggCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACa                                                                                                                           >  2:1027932/1‑150 (MQ=255)
                               gAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTg                                                                                                                      <  2:376629/150‑1 (MQ=255)
                                aaGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTgg                                                                                                                     <  2:307764/150‑1 (MQ=255)
                                aaGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTgg                                                                                                                     <  1:236826/150‑1 (MQ=255)
                                 aGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGa                                                                                                                    <  2:985804/150‑1 (MQ=255)
                                      gTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaataat                                                                                                               <  1:891105/150‑1 (MQ=255)
                                       tGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATg                                                                                                              <  1:1348937/150‑1 (MQ=255)
                                       tGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATg                                                                                                              <  2:937807/150‑1 (MQ=255)
                                         cTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGaa                                                                                                            <  1:295089/150‑1 (MQ=255)
                                           tCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGt                                                                                                          <  1:1430762/150‑1 (MQ=255)
                                           tCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGt                                                                                                          <  1:785276/150‑1 (MQ=255)
                                             aaaCGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTgg                                                                                                        <  1:363979/150‑1 (MQ=255)
                                                   tGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTgg                                                                                                  <  1:1372115/150‑1 (MQ=255)
                                                   tGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTgg                                                                                                  <  2:580895/150‑1 (MQ=255)
                                                    gccgccATTTCCCTGAACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACTGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTTGAATAATGAAGTGGGACTTGa                                                                                                 <  2:750234/150‑1 (MQ=255)
                                                    gccgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGa                                                                                                 >  2:1011695/1‑150 (MQ=255)
                                                     ccgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGAt                                                                                                <  2:390731/150‑1 (MQ=255)
                                                      cgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGAtt                                                                                               >  2:699022/1‑150 (MQ=255)
                                                       gccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGAttt                                                                                              >  2:1216706/1‑150 (MQ=255)
                                                           ttGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCtt                                                                                          <  2:976251/150‑1 (MQ=255)
                                                            tGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGGAATTACTACCCCACCACGGACGGTCGGGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCttt                                                                                         <  1:1425475/150‑1 (MQ=255)
                                                             gCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCACAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTg                                                                                        <  1:1321590/150‑1 (MQ=255)
                                                               ccTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGAt                                                                                      <  1:1423367/150‑1 (MQ=255)
                                                                 tGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGc                                                                                    <  1:1027932/150‑1 (MQ=255)
                                                                 tGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGc                                                                                    >  1:424517/1‑150 (MQ=255)
                                                                   tACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCg                                                                                  <  2:1214580/150‑1 (MQ=255)
                                                                     cGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGAt                                                                                >  2:707572/1‑150 (MQ=255)
                                                                         aGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGaaa                                                                            >  1:1310003/1‑150 (MQ=255)
                                                                          gTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAAc                                                                           >  2:1036413/1‑150 (MQ=255)
                                                                            ttAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGc                                                                         <  2:1309192/150‑1 (MQ=255)
                                                                              aGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTg                                                                       <  1:314594/150‑1 (MQ=255)
                                                                               ggCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGt                                                                      >  1:396189/1‑150 (MQ=255)
                                                                                gCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCAAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGtt                                                                     >  1:467977/1‑150 (MQ=255)
                                                                                   ttAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAc                                                                  <  2:254205/150‑1 (MQ=255)
                                                                                      aCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCtt                                                               <  2:424517/150‑1 (MQ=255)
                                                                                      aCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCtt                                                               >  1:607109/1‑150 (MQ=255)
                                                                                      aCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACCCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAAAGAACTGGCACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCtt                                                               >  2:629043/1‑150 (MQ=255)
                                                                                       cGAGGCGACGATACGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCttt                                                              >  2:484743/1‑150 (MQ=255)
                                                                                          ggCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGAc                                                           >  2:1212491/1‑150 (MQ=255)
                                                                                           gCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACt                                                          <  1:463837/150‑1 (MQ=255)
                                                                                           gCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACt                                                          <  1:1225038/150‑1 (MQ=255)
                                                                                             gacgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCa                                                        >  1:1222545/1‑150 (MQ=255)
                                                                                               cgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCAtt                                                      >  2:724602/1‑150 (MQ=255)
                                                                                                  ttCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCAc                                                   <  1:162494/150‑1 (MQ=255)
                                                                                                       cgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTATTTACCTTTGACTCATTCACCGGcc                                              >  1:1357131/1‑150 (MQ=255)
                                                                                                        gcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCt                                             <  1:1212491/150‑1 (MQ=255)
                                                                                                                  ccccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGATTACCTTTGACt                                                          <  2:898724/127‑1 (MQ=255)
                                                                                                                  ccccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGATTACCTTTGACt                                                          >  1:898724/1‑127 (MQ=255)
                                                                                                                   cccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGt                                  >  1:1065932/1‑150 (MQ=255)
                                                                                                                     caccacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGttt                                >  1:286830/1‑150 (MQ=255)
                                                                                                                      accacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTTTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTc                               <  2:1357131/150‑1 (MQ=255)
                                                                                                                      accacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTc                               <  2:1065932/150‑1 (MQ=255)
                                                                                                                        cacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCtt                             >  2:1286396/1‑150 (MQ=255)
                                                                                                                         acggacggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTg                            <  2:19435/150‑1 (MQ=255)
                                                                                                                                  cgcgAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCg                   >  2:918231/1‑150 (MQ=255)
                                                                                                                                   ttggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGt                  >  2:1384983/7‑150 (MQ=255)
                                                                                                                                     ggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGCCTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTcc                >  1:1318985/5‑149 (MQ=255)
                                                                                                                                       aGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCAcc              <  1:754563/150‑1 (MQ=255)
                                                                                                                                         aCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCtt            >  1:693800/1‑150 (MQ=255)
                                                                                                                                             catcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGTAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGTAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTAcc        <  1:55526/150‑1 (MQ=255)
                                                                                                                                             catcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTAcc        <  2:467977/150‑1 (MQ=255)
                                                                                                                                              atcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCg       <  1:707572/150‑1 (MQ=255)
                                                                                                                                                catcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCGCt     <  1:559015/150‑1 (MQ=255)
                                                                                                                                                  tcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGCCTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCGCTGa   >  2:408384/1‑150 (MQ=255)
                                                                                                                                                   catGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCGCTGaa  <  1:629043/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCGCTGAA  >  NC_000913/4608460‑4608756

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