Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,113,475 G→A 7.3% L61L (TTG→TTA) ‡ uspD → universal stress protein D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,113,4750GA 7.3% 200.2 / 5.7 78L61L (TTG→TTA) ‡uspDuniversal stress protein D
Reads supporting (aligned to +/- strand):  ref base G (46/25);  new base A (3/3);  total (50/28)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

GAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCT  >  NC_000913/4113335‑4113617
                                                                                                                                            |                                                                                                                                               
gaagaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCTCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgtt                                                                                                                                               >  1:572455/1‑143 (MQ=255)
gaagaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCGCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgtt                                                                                                                                               <  1:793441/143‑1 (MQ=255)
  agaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACGTAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTga                                                                                                                                             >  2:675646/1‑143 (MQ=255)
       tcTTGCTCGTCAATAGAGCCCGGGAGCTCGCCGGACATAATGCCGCTCCCCTGACGTTAATTCATATTGATGATCGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAAt                                                                                                                                        <  2:606190/142‑1 (MQ=255)
        cTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                                                                       >  1:451920/1‑143 (MQ=255)
          tACGGGTGAAGAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatc                                                                                                                                     <  1:698940/143‑3 (MQ=255)
           aCTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGCTAATTCATATTGAAGATGGCTTAAGCGAGATCTATCCGGGCATCTACTTCCCTGCAAAAGACGATATTCTTCAATTGTTGAAGAATAAGt                                                                                                                                    >  2:340491/1‑143 (MQ=255)
            cgtgcgcaTAAAGCCCTTGACCTCGCCAGACATAATGACGCGCACCTGACGTTAATTCATATTGATGATGGCTTACGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTc                                                                                                                                   <  2:575833/136‑1 (MQ=255)
             tGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTGTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                                                                  <  1:745411/143‑6 (MQ=255)
             tGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGACGATATTCTTCAATTGTTGAGGAATAAGTCg                                                                                                                                  >  1:760315/1‑143 (MQ=255)
                tGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCCACAGAAGATATTCTTCAATTGTTGCAGAATCACTCggct                                                                                                                               >  1:144504/1‑141 (MQ=255)
                 gAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAGTAAGTGGGATa                                                                                                                              >  1:109385/1‑143 (MQ=255)
                 gAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAACACTAAGTCGGATa                                                                                                                              >  2:163212/1‑143 (MQ=255)
                 gAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTCAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGCTa                                                                                                                              >  1:331162/1‑143 (MQ=255)
                  aaTAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAACATATTCTCCAATTGTTGCAGAATAAGTCGGATaa                                                                                                                             >  2:608195/1‑143 (MQ=255)
                   agataggccTGGAGCTCGTCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAAc                                                                                                                            <  2:402517/136‑1 (MQ=255)
                       aGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTACTTCATATTGATGATGGCCTAAGAGCGTTGTATATGGGTATCTACTTCCATGCAACAGAAGATATTCTTCAATTATTGAAGAACAAGTCGGAAAACaacc                                                                                                                        >  2:353771/1‑141 (MQ=255)
                        gCCCTGGATCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCt                                                                                                                       <  1:69360/143‑1 (MQ=255)
                         cccTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTg                                                                                                                      >  1:94838/1‑143 (MQ=255)
                          ccTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCAGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGt                                                                                                                     >  2:792348/1‑143 (MQ=255)
                           cTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGta                                                                                                                    >  2:491968/1‑143 (MQ=255)
                            tGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGATGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATGGTTGAAGACTAAGTCGGATAACAAGATGtat                                                                                                                   >  2:840826/1‑143 (MQ=255)
                             ggAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                                                                  >  1:523151/1‑143 (MQ=255)
                                 cTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATGTACTTCCCTGCAACAGAAGATATTCTTCAATCGTTGAAGAATAAGTCGGATAAGAAGCTGTATaaaat                                                                                                              >  2:51953/1‑141 (MQ=255)
                                 cTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTGCCTGCCCCAGAAGATATTCTTCAATTGTTGCAGAATAAGTCGGTTAACACGCTGTATaaaat                                                                                                              >  2:58074/1‑141 (MQ=255)
                                  tCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTg                                                                                                             >  1:239496/1‑143 (MQ=255)
                                   cGCCAGACATAATGACGCTCACCTGACGTTAATTCATACTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGa                                                                                                            <  1:434193/143‑1 (MQ=255)
                                    gtcAGACATAATGACGCTCACCTGACGTTAATACAGATTGATGAGGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAc                                                                                                           <  1:5907/141‑1 (MQ=255)
                                    gCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTAGAAGAATAAGACGGATAACAAACTGTATAAAATGAc                                                                                                           >  2:465480/1‑143 (MQ=255)
                                      cAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGACTAAGTCGGATAACAAGCTGTATAAACTGACgc                                                                                                         >  1:275707/1‑142 (MQ=255)
                                         aCATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGcaaa                                                                                                      >  1:72964/1‑143 (MQ=255)
                                         aCATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGagaa                                                                                                      >  1:276435/1‑143 (MQ=255)
                                            tAATGACGCTCACCTGACGTTAATTCATATGTATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                                                                   <  1:304567/143‑1 (MQ=255)
                                             aaTGTCGCTCACCTTACGTTAATTCATATTGTTGATGGCTGAAGCGAGTTGTATCCGGGTATCTACTTCCCGGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAatat                                                                                                  <  2:187732/143‑1 (MQ=255)
                                             aaTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAatat                                                                                                  <  2:458163/143‑1 (MQ=255)
                                              aTGACGCTCACCTGACTTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATAtt                                                                                                 <  2:125516/143‑1 (MQ=255)
                                              aTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCCACAGAAGATATTCTTCACTTGTTGAAGAATAAGCCCGATAACAAGCTGTATAAACTGACGAAAAATAtt                                                                                                 >  2:623503/1‑143 (MQ=255)
                                              aTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGTAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATAtt                                                                                                 <  2:711206/143‑1 (MQ=255)
                                                   gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCACTg                                                                                            >  2:559915/1‑143 (MQ=255)
                                                    ctctCCTGACGTTAATTCCTATTGATCATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                                                                           <  1:709277/143‑1 (MQ=255)
                                                     ggaCCTGAGGGTAATTCATAGTGATGATGGCTTAAGAGAGGTGTATCCGGGGATCTACTTCCCTGCAACAGAAGATATGCTTCAATTGTTGAAGAATAAGTGGGAGAACAAGGTGTATAAACTGACGAAAAATATTCAATGGc                                                                                          <  1:870984/141‑1 (MQ=255)
                                                      cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGcc                                                                                         >  2:167262/1‑143 (MQ=255)
                                                      cACCGGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGcc                                                                                         <  1:536538/143‑1 (MQ=255)
                                                         cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACCTCCCTGCAACAGCAGATATTCTCCATTCGAAAGAAAATAAGTCGGATAACAAGCTGTAGAAACTGACGAAAAATATTCAATGGCCGaa                                                                                      <  2:628711/143‑1 (MQ=255)
                                                            aCGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGAc                                                                                   >  1:711185/1‑143 (MQ=255)
                                                             cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAATAGGAGATATGCTTCAATTGTTGAAGAATAAGTCGGATAACAACCTGTATAAACTGACGAAAAATATTCAATGGccacacaca                                                                                  >  2:820492/1‑136 (MQ=255)
                                                                tAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAACCCGACGAAAAATATTCAATGGCAGAAGACaaaa                                                                              >  2:141252/1‑143 (MQ=255)
                                                                    tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAGAAATATTAAAGGGCCGAAGACAAAAATga                                                                          >  2:705660/1‑142 (MQ=255)
                                                                    tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGc                                                                          >  2:841528/1‑143 (MQ=255)
                                                                       tatTGATGATGGCTTAAGCGAGTTGTGTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTa                                                                       >  1:403924/1‑143 (MQ=255)
                                                                       tatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAAAAGAAGATATTCTTCAATAGTTGAAGAATAAGTCGGATAAAAAGCTGTAGAAACTGACGAAAAAAATCCAAATGACGAAGGAAAAACTGCATa                                                                       >  2:588892/1‑143 (MQ=255)
                                                                         ttGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTCCTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGAATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAtt                                                                     >  1:549202/1‑143 (MQ=255)
                                                                         ttGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCACTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAtt                                                                     <  2:715149/143‑1 (MQ=255)
                                                                         ttGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAACAACCTTCAATGGCCGAAGACACAACTGCGCAtt                                                                     >  1:158189/1‑143 (MQ=255)
                                                                             tgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAc                                                                 >  1:511257/1‑143 (MQ=255)
                                                                              gatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGACcg                                                                <  1:559915/143‑1 (MQ=255)
                                                                                   cTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCTGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACAGCGTATAGAACGCAgaa                                                           >  2:11719/1‑142 (MQ=255)
                                                                                    ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGCTATTCTTCAATTTTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATGGAAGGCGgaca                                                          >  2:583318/1‑141 (MQ=255)
                                                                                         cGAGTTGTATCCGGGTATCTACTTCCCTGCGACAGAAGATATTCTTCAATTGTTGAAGAATGAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAAGCGGAGAAATGc                                                     >  2:128597/1‑143 (MQ=255)
                                                                                         cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGc                                                     <  1:482232/143‑1 (MQ=255)
                                                                                           cgctGTACCCGTGTATCTACCTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAATCTGTATAACCTGACGCACAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCg                                                   <  1:263760/140‑1 (MQ=255)
                                                                                           aGTTGTATCCGGGTATCCACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAACATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCg                                                   <  1:848015/143‑1 (MQ=255)
                                                                                             ttGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATAGTTGATGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACCAAACTGCGTATTGACAGCGGAAAAACGCCGGa                                                 >  2:883734/1‑143 (MQ=255)
                                                                                                     gggTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACCGCGTATTGAACGCGGAGAAATGCCGGAAACActgc                                         >  1:575913/1‑143 (MQ=255)
                                                                                                      ggTATGTACTTTCCTGCAACAGAAGAGATTCTTCAATTGTGGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACActgct                                        <  2:625407/143‑1 (MQ=255)
                                                                                                           cTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAACAACACCGGCGGTTTGCGCGCGAAGCACTGCCGGAAGCACTGCTGGaga                                   >  1:197589/1‑141 (MQ=255)
                                                                                                             aCTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAAtt                                 >  1:813867/1‑143 (MQ=255)
                                                                                                              cTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTa                                <  1:167262/143‑1 (MQ=255)
                                                                                                               ttCCCTGCAACAGAAGATACTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTAt                               <  1:17741/143‑1 (MQ=255)
                                                                                                                 ccttGCAACAGAAGATATTCTTCAATTGTTGAATAATTAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGc                             <  2:474017/140‑1 (MQ=255)
                                                                                                                  ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCa                            >  1:298002/1‑143 (MQ=255)
                                                                                                                   cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCaa                           >  1:861271/1‑143 (MQ=255)
                                                                                                                   cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCAGAAGCACTGCTGGGAATTATGcgt                           >  1:682244/1‑141 (MQ=255)
                                                                                                                          agaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAAAATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGCAAACACCGCTGGAAATTAAGCAAAAAGAGc                    >  2:683146/1‑143 (MQ=255)
                                                                                                                           gaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCa                   <  1:657173/143‑1 (MQ=255)
                                                                                                                             acttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGCCCGAAGACAAAACTGCGTATTGAACGCTGAGGAATGCCGGGAACACTGCTGGAGCTTATGCAGAAAGCGCAGt                 >  1:240651/4‑143 (MQ=255)
                                                                                                                               ttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCACAAAGAGCAGTgg               >  1:449107/2‑142 (MQ=255)
                                                                                                                                      ttaaTTGTTGAAGAATAAGTCTGATAACAAGCTGTATAAACGGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGAcctcc        <  1:822266/141‑1 (MQ=255)
                                                                                                                                        aaTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGAcctccct      >  2:320508/1‑141 (MQ=255)
                                                                                                                                        aaTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGAcctccat      >  2:170262/1‑141 (MQ=255)
                                                                                                                                          tggttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGt    <  2:173406/141‑1 (MQ=255)
                                                                                                                                           tgtggAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTc   <  1:768371/143‑1 (MQ=255)
                                                                                                                                            gttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACGGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCt  <  1:463744/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                               
GAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCT  >  NC_000913/4113335‑4113617

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: