Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,372,620 A→C L532F (TTA→TTC)  groL → Cpn60 chaperonin GroEL, large subunit of GroESL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,372,6200AC98.5% 251.8 / ‑5.0 65L532F (TTA→TTC) groLCpn60 chaperonin GroEL, large subunit of GroESL
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (36/28);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTCATCTTTCTAGTATAAGATTCAGACACGGAC  >  NC_000913/4372472‑4372768
                                                                                                                                                    |                                                                                                                                                    
aagaaTACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCgg                                                                                                                                                    >  1:444380/1‑151 (MQ=255)
 agaaTACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGgc                                                                                                                                                   <  2:106115/151‑1 (MQ=255)
        ggCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGctgctg                                                                                                                                            >  2:169930/1‑151 (MQ=255)
                gATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTAt                                                                                                                                     >  1:212995/1‑150 (MQ=255)
                    gACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGCATGggc                                                                                                                                 >  2:23862/1‑150 (MQ=255)
                    gACATGGGAATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTAAGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACTTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATggg                                                                                                                                  >  2:1198445/1‑149 (MQ=255)
                       aTGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCa                                                                                                                             <  2:967868/151‑1 (MQ=255)
                          ggTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATgg                                                                                                                          >  2:505599/1‑151 (MQ=255)
                               ccTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTgg                                                                                                                      <  1:1019965/150‑1 (MQ=255)
                                 tGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCAt                                                                                                                   >  1:767759/1‑151 (MQ=255)
                                 tGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCAt                                                                                                                   >  2:284426/1‑151 (MQ=255)
                                  ggATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCa                                                                                                                    >  2:59948/1‑149 (MQ=255)
                                           caaAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCatg                                                                                                         >  1:949020/1‑151 (MQ=255)
                                              aGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCatgatg                                                                                                      >  1:367646/1‑151 (MQ=255)
                                               gTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGt                                                                                                     <  2:620101/151‑1 (MQ=255)
                                                tAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGGGGCgtg                                                                                                         >  1:631180/1‑146 (MQ=255)
                                                  aCTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTaa                                                                                                   >  2:387185/1‑150 (MQ=255)
                                                      gTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGcc                                                                                              >  1:631089/1‑151 (MQ=255)
                                                       ttCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGccc                                                                                             >  1:1090784/1‑151 (MQ=255)
                                                       ttCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGccc                                                                                             >  2:92515/1‑151 (MQ=255)
                                                        tctgctctgcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGccc                                                                                             <  1:636858/150‑1 (MQ=255)
                                                           gctctgcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGc                                                                                          <  1:1135387/150‑1 (MQ=255)
                                                             tctgcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCAcc                                                                                       >  2:105619/1‑151 (MQ=255)
                                                              ctgcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCAcc                                                                                       >  1:749475/1‑150 (MQ=255)
                                                                  aGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTGCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCg                                                                                    <  1:1018979/149‑1 (MQ=255)
                                                                     aCGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGa                                                                                >  1:455764/1‑150 (MQ=255)
                                                                       gCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAAt                                                                             >  1:319228/1‑151 (MQ=255)
                                                                         aGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATa                                                                            >  2:233778/1‑150 (MQ=255)
                                                                                gtgGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATaaacaa                                                                       <  1:813559/148‑1 (MQ=255)
                                                                                gtgGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATaaacaa                                                                       <  2:450170/148‑1 (MQ=255)
                                                                                gtgGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAAcc                                                                    >  1:891110/1‑151 (MQ=255)
                                                                                      ggcCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCgg                                                               <  1:1038942/150‑1 (MQ=255)
                                                                                          tgatgaTCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGc                                                             >  1:100540/1‑148 (MQ=255)
                                                                                           gatgatCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGaa                                                         <  2:477489/151‑1 (MQ=255)
                                                                                           gatgatCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGa                                                          >  1:925393/1‑150 (MQ=255)
                                                                                             tgatCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAAt                                                       >  2:445342/1‑151 (MQ=255)
                                                                                              gatCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGGAATAAACAAACCCCCGGGCAGAAAt                                                       <  1:803085/150‑1 (MQ=255)
                                                                                                 caccacCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCt                                                   <  1:268171/151‑1 (MQ=255)
                                                                                                    caccGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTggg                                                >  1:823711/1‑151 (MQ=255)
                                                                                                      ccGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTggggg                                              <  1:585710/151‑1 (MQ=255)
                                                                                                      ccGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTggggg                                              <  2:926000/151‑1 (MQ=255)
                                                                                                      ccGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTgggg                                               <  2:506621/150‑1 (MQ=255)
                                                                                                         aaTGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGttt                                           <  2:777615/151‑1 (MQ=255)
                                                                                                         aaTGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGttt                                           >  2:698582/1‑151 (MQ=255)
                                                                                                            gCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGttttt                                         >  2:200722/1‑150 (MQ=255)
                                                                                                            gCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGttttt                                         >  1:552048/1‑150 (MQ=255)
                                                                                                            gCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCCGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGttttt                                         >  2:611900/1‑150 (MQ=255)
                                                                                                              aTGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCtt                                      <  2:1090712/151‑1 (MQ=255)
                                                                                                               tGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGGTTTTCttt                                     >  2:942114/1‑151 (MQ=255)
                                                                                                                ggTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCtttt                                    >  1:1142458/1‑151 (MQ=255)
                                                                                                                 gTTACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTg                                   <  1:610735/151‑1 (MQ=255)
                                                                                                                   tACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTgg                                  >  2:9413/1‑150 (MQ=255)
                                                                                                                   tACCGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTgg                                  >  2:963397/1‑150 (MQ=255)
                                                                                                                     ccGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTc                                >  2:484110/1‑150 (MQ=255)
                                                                                                                      cGACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTCAt                              <  1:752799/151‑1 (MQ=255)
                                                                                                                       gACCTGCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTCa                               <  1:613330/149‑1 (MQ=255)
                                                                                                                            gCCGAAAAACGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTCATCttt                          <  1:272678/149‑1 (MQ=255)
                                                                                                                                aaaaaCGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTCATCTTTCTAGta                    <  1:959801/151‑1 (MQ=255)
                                                                                                                                   aaCGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTCATCTTTCTAGTATaa                 <  1:1101011/151‑1 (MQ=255)
                                                                                                                                   aaCGATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTCATCTTTCTAGTATaa                 <  2:212986/151‑1 (MQ=255)
                                                                                                                                      gATGCAGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTCATCTTTCTAGTATAAGAt              <  2:556173/151‑1 (MQ=255)
                                                                                                                                           aGCTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTCATCTTTCTAGTATAAGATTCAGa         <  2:1101960/151‑1 (MQ=255)
                                                                                                                                             cTGACTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTCATCTTTCTAGTATAAGATTCAGac        >  2:569533/1‑150 (MQ=255)
                                                                                                                                                aCTTCGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTCATCTTTCTAGTATAAGATTCAGACACgg    <  2:869102/151‑1 (MQ=255)
                                                                                                                                                  ttAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCCGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTCATCTTTCTAGTATAAGATTCAGACACGgac  >  2:169434/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAATAAACAAACCCCCGGGCAGAAATGTCTGGGGGTTTTTCTTTTGGTCATCTTTCTAGTATAAGATTCAGACACGGAC  >  NC_000913/4372472‑4372768

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: ‑