Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,956,720 T→G K722T (AAA→ACA)  ptrA ← protease III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,956,7200TG95.3% 208.5 / 0.4 64K722T (AAA→ACA) ptrAprotease III
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (34/27);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

ATTGTTCTTCGGTACGCAACTGATTGTAGAACCACGGCTGTACGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTTTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACTTGCGAGAGCATCTGCGCGGGCATAATCGCCTGCTCAAACGCTTTACCCTTTTCTGCGG  >  NC_000913/2956496‑2956945
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                 
aTTGTTCTTCGGTACGCAACTGATTGTAGAACCACGGCTGTACGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATcg                                                                                                                                                                                                                        <  2:563587/236‑1 (MQ=255)
      cttcGGTACGCAACTGATTGTAGAACCACGGCTGTACGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGcc                                                                                                                                                                                                                  >  1:164029/1‑236 (MQ=255)
         cGGTACGCAACTGATTGTAGAACCACGGCTGTACGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGcc                                                                                                                                                                                                                  >  2:83180/1‑233 (MQ=255)
           gTACGCAACTGATTGTAGAACCACGGCTGTACGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCgt                                                                                                                                                                                                             >  1:426081/1‑236 (MQ=255)
                      aTTGTAGAACCACGGCTGTACGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTggg                                                                                                                                                                                                  >  2:21541/1‑236 (MQ=255)
                                  cGGCTGTACGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATg                                                                                                                                                                                       >  1:131781/1‑235 (MQ=255)
                                  cGGCTGTACGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATg                                                                                                                                                                                       >  2:763189/1‑235 (MQ=255)
                                    gCTGTACGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCCCTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCAGTCAt                                                                                                                                                                                        >  1:691218/1‑232 (MQ=255)
                                         aCGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGa                                                                                                                                                                                <  1:579211/235‑1 (MQ=255)
                                          cGATCTGCCCCAACCGGAAGCTATAGGCTGAGCTGGTGTTTTCATCGTAGCCAGTCGGTACAAATTCCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATa                                                                                                                                                                              <  1:351503/236‑1 (MQ=255)
                                          cGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATa                                                                                                                                                                              <  1:83180/236‑1 (MQ=255)
                                             tCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATTGCGTTCCAGCGTTGTTGCCTGGGCCTCGGTCCTGTTTCCGATAAcc                                                                                                                                                                           >  2:761933/1‑236 (MQ=255)
                                             tCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTTCCGATAAcc                                                                                                                                                                           >  2:233694/1‑236 (MQ=255)
                                               tGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGAAACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTTCCTGGGCCTCGGTCATGTTGCCGATAATCa                                                                                                                                                                          >  2:690465/1‑235 (MQ=255)
                                                 ccccAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATa                                                                                                                                                                        >  1:179347/1‑235 (MQ=255)
                                                         gagagCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGt                                                                                                                                                               >  1:564806/1‑236 (MQ=255)
                                                                 taGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCCAGcccc                                                                                                                                                       >  1:131859/1‑236 (MQ=255)
                                                                  aGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGcccc                                                                                                                                                       <  2:131781/235‑1 (MQ=255)
                                                                    gCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTTTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGa                                                                                                                                                    >  2:650007/1‑236 (MQ=255)
                                                                    gCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGtt                                                                                                                                                                                     >  2:183144/1‑203 (MQ=255)
                                                                    gCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGtt                                                                                                                                                                                     <  1:183144/203‑1 (MQ=255)
                                                                        aGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGcccc                                                                                                                                                       >  2:79631/1‑229 (MQ=255)
                                                                        aGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGcccc                                                                                                                                                       <  1:79631/229‑1 (MQ=255)
                                                                                    cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGa                                                                                                                                    <  1:763189/236‑1 (MQ=255)
                                                                                    cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGa                                                                                                                                    <  2:426081/236‑1 (MQ=255)
                                                                                        gTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCTCGATAgg                                                                                                                                >  2:43754/1‑236 (MQ=255)
                                                                                          aGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGc                                                                                                                               >  1:385707/1‑235 (MQ=255)
                                                                                           gCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCa                                                                                                                             <  2:564806/236‑1 (MQ=255)
                                                                                              aGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCa                                                                                                                          <  1:21541/236‑1 (MQ=255)
                                                                                                     aCAGATACCGCTTCCAGTGCGGAGTCGTGACTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCAGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCttt                                                                                                                   <  2:691218/236‑1 (MQ=255)
                                                                                                     aCAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTTTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCttt                                                                                                                   <  1:650007/236‑1 (MQ=255)
                                                                                                        aaTACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCttt                                                                                                                   <  2:617484/233‑1 (MQ=255)
                                                                                                                                gTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAAttt                                                                                            >  2:563783/1‑232 (MQ=255)
                                                                                                                                 tGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGc                                                                                       >  2:466220/1‑236 (MQ=255)
                                                                                                                                        aCCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACGCAACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCAt                                                                                <  2:608379/236‑1 (MQ=255)
                                                                                                                                         ccGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCAtc                                                                               <  1:563783/236‑1 (MQ=255)
                                                                                                                                           ggCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGcccc                                                                                                                                                       >  2:558351/1‑162 (MQ=255)
                                                                                                                                           ggCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGcccc                                                                                                                                                       <  1:558351/162‑1 (MQ=255)
                                                                                                                                             cTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTcgc                                                                           >  2:614983/1‑236 (MQ=255)
                                                                                                                                                 tacAAAGATGACTGATTGGTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGACCCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCgag                                                                        <  2:373643/233‑1 (MQ=255)
                                                                                                                                                            cGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGa                                                                                                                                    <  2:132669/164‑1 (MQ=255)
                                                                                                                                                            cGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGa                                                                                                                                    >  1:132669/1‑164 (MQ=255)
                                                                                                                                                            cGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACt                                                            >  1:183846/1‑236 (MQ=255)
                                                                                                                                                                   ttattaaTCGACCACTACATCGTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACTTGCgag                                                      <  2:164029/229‑1 (MQ=255)
                                                                                                                                                                      tttATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACTTGCgaga                                                     >  1:488527/1‑233 (MQ=255)
                                                                                                                                                                            gACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACTTGCGAGAGCATCTgcg                                            <  1:466220/236‑1 (MQ=255)
                                                                                                                                                                             aCCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTTTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACtt                                                           <  2:27948/220‑1 (MQ=255)
                                                                                                                                                                             aCCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTcgc                                                                                                                                      <  2:280631/145‑1 (MQ=255)
                                                                                                                                                                             aCCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTcgc                                                                                                                                      >  1:280631/1‑145 (MQ=255)
                                                                                                                                                                             aCCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGa                                                                                                                                    <  1:339636/147‑1 (MQ=255)
                                                                                                                                                                             aCCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGa                                                                                                                                    >  2:339636/1‑147 (MQ=255)
                                                                                                                                                                             aCCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACtt                                                           >  1:27948/1‑220 (MQ=255)
                                                                                                                                                                                            tttCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCAtctc                                                                             >  2:100634/1‑187 (MQ=255)
                                                                                                                                                                                            tttCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCAtctc                                                                             <  1:100634/187‑1 (MQ=255)
                                                                                                                                                                                            tttCGACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACTTGCGAGAGCATCTGCGCGGGCATAATCGCCTg                            >  2:201059/1‑236 (MQ=255)
                                                                                                                                                                                                gACACCACTCTGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACTTGCGAGAGCATCTGCGCGGGCATAATCGCCTGCTCa                        >  2:471269/1‑236 (MQ=255)
                                                                                                                                                                                                        tctGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACTTGCGAGAGCATCTGCGCGGGCATAATCGCCTGCTCaa                       >  2:398386/1‑229 (MQ=255)
                                                                                                                                                                                                        tctGAACCATCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACTTGCGAGAGCATCTGCGCGGGCATAATCGCCTGCTCaa                       <  1:398386/229‑1 (MQ=255)
                                                                                                                                                                                                                 tCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCAc                                                                                                                         <  2:582349/122‑1 (MQ=255)
                                                                                                                                                                                                                 tCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCAc                                                                                                                         >  1:582349/1‑122 (MQ=255)
                                                                                                                                                                                                                 tCAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACTTGCGAGAGCATCTGCGCGGGCATAATCGCCTGCTCAAACGCTTTACCCTTTTc       >  2:190740/1‑236 (MQ=255)
                                                                                                                                                                                                                  cAGCGCCCAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACTTGCGAGAGCATCTGCGCGGGCATAATCGCCTGCTCAAACGCTTTACCCTTTTCt      >  2:83768/1‑236 (MQ=255)
                                                                                                                                                                                                                         cAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACTTGCGAGAGCATCTGCGCGGGCATAATCGCCTGCTCAAACGCTTTACCCTTTTCTGCgg  <  1:614983/233‑1 (MQ=255)
                                                                                                                                                                                                                         cAACTGTGTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACTTGCGAGAGCATCTGCGCGGGCATAATCGCCTGCTCAAACGCTTTACCCTTTTCTGCg   <  1:535306/232‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                 
ATTGTTCTTCGGTACGCAACTGATTGTAGAACCACGGCTGTACGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAACCATCAGCGCCCAACTGTTTTTGCACATCGCGTGCCAGCGTTGTTGCCTGGGCCTCGGTCATGTTGCCGATAACCATAAACTCTGGTCGAGCCCCTGATTTTAAGGCGTCGCGATAGGCCAGCACCTCTTTCAACGTAATGGAGGGCAAAATTTTACGCCGTTCATCTCGCGAGAAGTACGGCACTTGCGAGAGCATCTGCGCGGGCATAATCGCCTGCTCAAACGCTTTACCCTTTTCTGCGG  >  NC_000913/2956496‑2956945

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