Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,494,883 (A)7→8 coding (261/564 nt) idnK → D‑gluconate kinase, thermosensitive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,494,8761.A97.6% 458.3 / 2.5 127coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base A (57/66);  total (61/66)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAT  >  NC_000913/4494778‑4494969
                                                                                                   |                                                                                                
   tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aa                                                                                            >  1:1081468/1‑100 (MQ=255)
   tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aa                                                                                            >  2:2516695/1‑100 (MQ=255)
      tgcataccGACTTCCCTTGCTAGAACGCTTTAATGATGATTCATACAGTGTTTATACTAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑a                                                                                             <  1:2797811/89‑1 (MQ=255)
      tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaaa                                                                                         >  2:50323/1‑100 (MQ=255)
      tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaaa                                                                                         >  2:1297264/1‑100 (MQ=255)
       gaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaaaa                                                                                        >  1:792551/1‑99 (MQ=255)
        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGTATTTATTGTCTGTTCATCATT‑‑aa‑aaaaaac                                                                                       <  1:1857594/100‑3 (MQ=255)
        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaaaac                                                                                       <  2:2233937/100‑3 (MQ=255)
        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaaaac                                                                                       <  1:135351/100‑3 (MQ=255)
        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaaaac                                                                                       <  2:2027916/100‑3 (MQ=255)
          gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaaaacag                                                                                     <  2:1682996/100‑5 (MQ=255)
           aTCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaaaacagt                                                                                    >  2:864337/1‑95 (MQ=255)
            tCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTa                                                                                   <  2:428690/100‑1 (MQ=255)
            tCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTa                                                                                   >  1:1667715/1‑100 (MQ=255)
            tCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTa                                                                                   <  2:801002/100‑1 (MQ=255)
             cGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTAt                                                                                  >  2:600876/1‑100 (MQ=255)
               aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCg                                                                                >  2:2541397/1‑100 (MQ=255)
               aCTTCACTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGAAAATAATa‑‑aa‑aaaaaacag                                                                                     <  1:263093/95‑5 (MQ=255)
                 ttCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTg                                                                              <  1:1373021/100‑1 (MQ=255)
                 ttCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATAATTA‑AA‑AAAAAC‑AGTATCGt                                                                               <  1:1196131/99‑1 (MQ=255)
                  tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAACAGAATGAAACAGGATTTATTGTTTGTTCATCATT‑‑AACAAAAAC‑AGTAtagaga                                                                             >  2:2564870/1‑95 (MQ=255)
                   cccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGa                                                                             <  1:1304001/99‑1 (MQ=255)
                      tggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATAtt                                                                         >  1:1270290/1‑100 (MQ=255)
                       tgctggAACGATTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATAttt                                                                        <  1:1016222/99‑1 (MQ=255)
                       ggctggAACGCTTAAATGATGCTTCATACAGTCTTTTTAAAATGTCTCCAACAGTATTTATTGTCTATACATCATT‑‑AA‑AAAAACGAGTAACGAGATCttt                                                                        >  2:315076/1‑100 (MQ=255)
                       ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATAttt                                                                        >  1:2039951/1‑100 (MQ=255)
                        gctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATAtttt                                                                       >  1:2729659/1‑100 (MQ=255)
                         ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTa                                                                      <  2:2488732/100‑1 (MQ=255)
                         ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTa                                                                      <  1:1404808/100‑1 (MQ=255)
                         ctagAACACTTAAATGATGCTTCATACAGTCTTTATAAAAAAAATAAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATAtttt                                                                       <  1:2541377/96‑1 (MQ=255)
                          tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTAc                                                                     <  2:1351673/100‑1 (MQ=255)
                          aatctcGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTAc                                                                     <  2:1616176/95‑1 (MQ=255)
                             aaCGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCa                                                                  >  1:2443540/1‑100 (MQ=255)
                               cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTAcagatc                                                                >  1:1616164/1‑95 (MQ=255)
                               cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAg                                                                <  2:210339/100‑1 (MQ=255)
                               cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAg                                                                <  1:2618365/100‑1 (MQ=255)
                                 cTTCAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGCTTTCTTGTCTGTTCATCATTC‑AA‑CAAAAC‑AGTATTGTGATAtctttc                                                                     >  1:2283004/1‑91 (MQ=255)
                                 cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAggg                                                              >  2:1024789/1‑100 (MQ=255)
                                   tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAggata                                                            >  1:1775013/1‑97 (MQ=255)
                                   tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTa                                                            >  1:2714353/1‑100 (MQ=255)
                                   tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTa                                                            >  1:578353/1‑100 (MQ=255)
                                    aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAg                                                           >  1:1488863/1‑100 (MQ=255)
                                          tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                     <  1:600871/100‑1 (MQ=255)
                                              tCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                 <  1:998619/100‑1 (MQ=255)
                                               cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                                >  2:35661/1‑100 (MQ=255)
                                                aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                               <  2:1234390/100‑1 (MQ=255)
                                                aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                               <  2:403110/100‑1 (MQ=255)
                                                 tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                              >  2:1446304/1‑100 (MQ=255)
                                                 tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                              >  1:433778/1‑100 (MQ=255)
                                                 tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                              >  1:1897572/1‑100 (MQ=255)
                                                  aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                             >  2:566588/1‑100 (MQ=255)
                                                   cAGTCTTTATAAAAAGAATGAAACAGGACGTATTGTCTGTTGATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGACACCATGTTCCttt                                            >  1:165626/1‑100 (MQ=255)
                                                     gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                          >  1:2645559/1‑100 (MQ=255)
                                                     gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                           >  2:2071085/1‑99 (MQ=255)
                                                      ncTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCct                                         >  1:2560387/2‑100 (MQ=255)
                                                       cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                        <  2:1481384/100‑1 (MQ=255)
                                                       cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                        <  2:1972402/100‑1 (MQ=255)
                                                         ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                      <  2:2719120/100‑1 (MQ=255)
                                                         ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                      <  1:355736/100‑1 (MQ=255)
                                                         ttATAAAAAGAATGAAAAAGGATTTATTGTCTGTTAATAATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                      <  1:1592162/100‑1 (MQ=255)
                                                          tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                     <  2:1004451/100‑1 (MQ=255)
                                                           ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                    >  1:483280/1‑100 (MQ=255)
                                                           ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                    <  2:2073861/100‑1 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                   >  1:1428600/1‑100 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                   >  2:2241619/1‑100 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                   <  1:2082231/100‑1 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                   >  1:1436755/1‑100 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                   >  1:362542/1‑100 (MQ=255)
                                                             aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                  <  2:1589988/100‑1 (MQ=255)
                                                              gatctaaTGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                 <  1:2428615/95‑1 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                 >  1:1866606/1‑100 (MQ=255)
                                                                ctgAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                               <  1:2617794/98‑1 (MQ=255)
                                                                aaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                >  1:2823613/1‑100 (MQ=255)
                                                                aaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATTGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                               <  2:792557/100‑1 (MQ=255)
                                                                   aaTGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGagatc                            >  2:2428635/1‑96 (MQ=255)
                                                                   aaTGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                            <  2:91230/100‑1 (MQ=255)
                                                                   aaTGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                            <  1:2071070/100‑1 (MQ=255)
                                                                    aTGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                           <  2:351743/100‑1 (MQ=255)
                                                                    aTGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGATAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                           <  2:1775025/100‑1 (MQ=255)
                                                                     tGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                          >  1:806214/1‑100 (MQ=255)
                                                                      gAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                         <  2:1261290/100‑1 (MQ=255)
                                                                        aaCAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                       <  1:864330/100‑1 (MQ=255)
                                                                        aaCAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                       <  1:2220548/100‑1 (MQ=255)
                                                                        aaCAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                       <  1:786938/100‑1 (MQ=255)
                                                                         aCAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                      >  2:2284033/1‑100 (MQ=255)
                                                                         aCAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                      >  1:1194570/1‑100 (MQ=255)
                                                                         aCAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                      <  1:1525648/100‑1 (MQ=255)
                                                                         aCAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                      <  2:2641735/100‑1 (MQ=255)
                                                                              aTTTATTGTGTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                 <  1:2494663/100‑1 (MQ=255)
                                                                              aTTTATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                 <  2:2798934/100‑1 (MQ=255)
                                                                                ttATTGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGAGACTATGAAACTa               <  2:1914314/100‑1 (MQ=255)
                                                                                   ttGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtc            >  2:2286591/1‑100 (MQ=255)
                                                                                    tGTCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct           >  2:1828822/1‑100 (MQ=255)
                                                                                      tCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg         >  1:2772777/1‑100 (MQ=255)
                                                                                      tCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg         >  1:2772767/1‑100 (MQ=255)
                                                                                      tCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg         <  2:2729683/100‑1 (MQ=255)
                                                                                      tCTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg         <  1:1024782/100‑1 (MQ=255)
                                                                                       cTGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc        >  2:1457591/1‑100 (MQ=255)
                                                                                        tGTTCATCATTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCATCTGGTAAGATAGCGACTATGAAACTATTCTcgc        >  1:308417/1‑99 (MQ=255)
                                                                                            catcatTA‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGTTTAGATGGCGACTATGAAACTATTCTCGCGCGaa   >  2:2645176/1‑100 (MQ=255)
                                                                                                   |                                                                                                
ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AA‑AAAAAC‑AGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAT  >  NC_000913/4494778‑4494969

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