Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,608,604 C→G 8.5% I140M (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,6040CG 8.5% 223.3 / 10.7 71I140M (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.10e-01

GAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCGCT  >  NC_000913/4608455‑4608753
                                                                                                                                                     |                                                                                                                                                     
gaACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatc                                                                                                                                                       <  1:122449/150‑1 (MQ=255)
 aaCTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatca                                                                                                                                                      <  1:296104/150‑1 (MQ=255)
  aCTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcat                                                                                                                                                     >  1:819142/1‑150 (MQ=255)
     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:162494/1‑150 (MQ=255)
                    ggCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAg                                                                                                                                   >  2:1250430/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                                                                                               <  2:580895/150‑1 (MQ=255)
                                                        tGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTgg                                                                                               <  1:1372115/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:424517/1‑150 (MQ=255)
                                                                      tGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGc                                                                                 <  1:1027932/150‑1 (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)
                                                                                                                                                     |                                                                                                                                                     
GAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCGCT  >  NC_000913/4608455‑4608753

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: