Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,494,883:1 (A)7→8 100% 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.8% 306.4 / ‑0.3 91coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (42/47);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 111 total aligned reads displayed.

CTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494780‑4494970
                                                                                                 |                                                                                                
cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaa                                                                                              >  2:2440973/1‑100 (MQ=255)
cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaa                                                                                              >  1:1521593/1‑100 (MQ=255)
cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaa                                                                                              >  2:659070/1‑100 (MQ=255)
  gatgaAGATCGACTTCCCTGGGTGGAACGCCTAAATTATGCTTCATACAGTCGTTATAAGAAGAATGAAACAGTATTTCTAGTCAGTTCATGATT‑‑caaa                                                                                             >  1:2508890/1‑99 (MQ=255)
  gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                            <  1:2531685/100‑1 (MQ=255)
  gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                            <  2:347027/100‑1 (MQ=255)
  gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                            <  1:309179/100‑1 (MQ=255)
  gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                            >  1:79320/1‑100 (MQ=255)
   atgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑a                                                                                          >  2:806355/1‑100 (MQ=255)
    tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑aa                                                                                         >  2:2113047/1‑100 (MQ=255)
    tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑aa                                                                                         >  2:1223690/1‑100 (MQ=255)
    tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑a                                                                                          >  2:1204026/1‑99 (MQ=255)
    tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGCTCATCATT‑‑aaaaa‑aa                                                                                         >  1:885319/1‑100 (MQ=255)
    agaAGATCGACTTCCCTGGCTGGAACGCTTTAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑aa                                                                                         >  1:1381982/2‑100 (MQ=255)
     gaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑aaa                                                                                        >  1:1816291/1‑99 (MQ=255)
      aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑aaac                                                                                       <  1:2283012/100‑3 (MQ=255)
      aaGATCGACTTCACTGGCTAGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑aaac                                                                                       <  1:177984/100‑3 (MQ=255)
       aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑aaaca                                                                                      <  2:2682279/100‑4 (MQ=255)
        gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑aaacag                                                                                     >  1:512750/1‑96 (MQ=255)
         aTCAACTTCCCTGGCTGTAACGCTTTGATAATCCTTTATAAAGACGTTAGAAAAAAAAAGAAACAGGATGTATAGCAGGTTCATCATT‑‑aaaaa‑aagcag                                                                                     <  1:729021/99‑5 (MQ=255)
          tCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTa                                                                                    <  2:1841985/100‑1 (MQ=255)
              cTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGt                                                                                <  2:636546/100‑1 (MQ=255)
              cTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATAATTA‑AAAAA‑AACAGTATCGt                                                                                <  1:571318/100‑1 (MQ=255)
                  ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGata                                                                            >  1:2007520/1‑100 (MQ=255)
                  ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGa                                                                              >  1:1057494/1‑98 (MQ=255)
                  ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAGAACAGTATCGTGata                                                                            >  2:2872239/1‑100 (MQ=255)
                   ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGatat                                                                           >  2:1268814/1‑100 (MQ=255)
                   ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGatat                                                                           >  2:1141410/1‑100 (MQ=255)
                    tggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATAtt                                                                          >  2:1594277/1‑100 (MQ=255)
                     ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATAttt                                                                         >  1:2716988/1‑100 (MQ=255)
                     ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATAttt                                                                         >  1:892116/1‑100 (MQ=255)
                       ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTa                                                                       >  2:1994396/1‑100 (MQ=255)
                       ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTa                                                                       <  2:575529/100‑1 (MQ=255)
                       ctggAACGCTTAAATGATGCTTCATACAGTATTCATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCATGATATTTTa                                                                       >  1:1914004/1‑100 (MQ=255)
                         ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACg                                                                     <  2:2046518/100‑1 (MQ=255)
                         ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACg                                                                     <  2:1052549/100‑1 (MQ=255)
                          gAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGc                                                                    <  2:1712396/100‑1 (MQ=255)
                               cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAggg                                                               >  1:1631646/1‑100 (MQ=255)
                               cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATCAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAgg                                                                <  2:160520/99‑1 (MQ=255)
                               cTTAAATCATGCTTCATACAGTCTTTATAAAAAGAATAAAACAGGATTTATTGTCTGTTCATAATTA‑AAAAA‑AACAGTATCGTGCTATTTTACGCAAg                                                                 >  1:900344/1‑98 (MQ=255)
                               cTTAAAGGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAg                                                                 <  2:2014361/98‑1 (MQ=255)
                                ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGt                                                              <  2:2259404/100‑1 (MQ=255)
                                taaaaTGATACTTAATACACTCTTAATAAAAATAATGAAAAAGGATTTATTTTCTGTTCATAATAA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAggg                                                               <  1:1138540/97‑1 (MQ=255)
                                  aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAg                                                            >  2:1790405/1‑100 (MQ=255)
                                   aaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGc                                                           >  2:1080054/1‑100 (MQ=255)
                                    atgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGcc                                                          <  2:1414309/100‑1 (MQ=255)
                                     tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                         <  1:1750324/100‑1 (MQ=255)
                                     tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                         <  1:2300462/100‑1 (MQ=255)
                                      gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                        <  1:20383/100‑1 (MQ=255)
                                      gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGcc                                                          >  1:2693175/1‑98 (MQ=255)
                                        tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGAGTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGccccc                                                       >  1:1561120/1‑99 (MQ=255)
                                          cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                    >  2:657158/1‑100 (MQ=255)
                                          cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                    <  2:2280609/100‑1 (MQ=255)
                                            tCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                  <  2:512753/100‑1 (MQ=255)
                                            tCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                  <  2:2240839/100‑1 (MQ=255)
                                             cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                                 <  2:2150654/100‑1 (MQ=255)
                                              aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                <  1:1524837/100‑1 (MQ=255)
                                               tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                               >  2:1396077/1‑100 (MQ=255)
                                               tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                               >  1:12636/1‑100 (MQ=255)
                                                 cAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAttt                                             >  2:2511688/1‑100 (MQ=255)
                                                  aatCTTTATAAAAACAATGAAACAGGATTTATTGTCTGTTAAAAAAAA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                            <  1:1986200/98‑1 (MQ=255)
                                                   gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                           <  1:2617729/100‑1 (MQ=255)
                                                   gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                           <  2:69458/100‑1 (MQ=255)
                                                    tgtttATAAAAAGAATGAAACAGGATCTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCct                                          <  1:275725/98‑1 (MQ=255)
                                                    tCTTTATAACAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATCTTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                           <  1:1917414/99‑1 (MQ=255)
                                                     cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                           <  1:607456/98‑1 (MQ=255)
                                                     cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         <  2:1632835/100‑1 (MQ=255)
                                                      tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         <  2:2041974/99‑1 (MQ=255)
                                                       ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                       <  1:2037484/100‑1 (MQ=255)
                                                        tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  2:850961/100‑1 (MQ=255)
                                                        tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  2:905897/100‑1 (MQ=255)
                                                         ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGTGTAGCCCCCATGTTCATTTCCTCTGGt                                     >  1:1298126/1‑100 (MQ=255)
                                                          taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCACCCATGTTCATTTCCTCTGGtt                                    >  1:2602028/1‑100 (MQ=255)
                                                          taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATAGTGATATTTTACGCAAGGGTTGCCCCCATGTTCATTTCCTCTGGt                                     >  1:1176622/1‑99 (MQ=255)
                                                           aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                   >  2:2361097/1‑100 (MQ=255)
                                                            aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  >  2:764991/1‑100 (MQ=255)
                                                              aaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                                <  1:1193007/100‑1 (MQ=255)
                                                               aGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATg                               <  2:441466/100‑1 (MQ=255)
                                                               aGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATg                               >  1:25828/1‑100 (MQ=255)
                                                                gAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTTTGGTTAGAt                                >  1:2502035/1‑98 (MQ=255)
                                                                 aaTGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                             >  2:1539850/1‑100 (MQ=255)
                                                                 aaTGAAACAGGATTTATTGTCTGTTCAACAGTAAAAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATgg                              <  2:405773/100‑1 (MQ=255)
                                                                  aTGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCAATGTTAATTTCCTCTGGGTAGATGGCg                            >  1:1131620/1‑100 (MQ=255)
                                                                     aaaCAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                         >  2:2789748/1‑100 (MQ=255)
                                                                     aaaCAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                         >  1:1547351/1‑100 (MQ=255)
                                                                       aCAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                       <  2:1882851/100‑1 (MQ=255)
                                                                         aGGATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCCACTATGa                     <  2:1529595/100‑1 (MQ=255)
                                                                           gATTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCAACAATATTCATTTACTCTGGTTCGATGGCGACTATGaaa                   >  1:1630958/1‑100 (MQ=255)
                                                                            aTTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                  <  2:2506925/100‑1 (MQ=255)
                                                                            aTTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                  <  1:174688/100‑1 (MQ=255)
                                                                            aTTTATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                  <  2:2717006/100‑1 (MQ=255)
                                                                            aTTTATTATCTGTTCAACATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                  <  1:606734/100‑1 (MQ=255)
                                                                              ttATTGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTa                <  1:1574747/100‑1 (MQ=255)
                                                                                  tGTCTGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct            <  2:1684211/100‑1 (MQ=255)
                                                                                  tGTCTGTTCATCATTA‑AAAAA‑AACAGTATCATGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct            >  2:807802/1‑100 (MQ=255)
                                                                                      tGTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCACCCATGTTAATTTTCTCTGGTTACATGGCGACTATGAAACTATTCTctcg        >  1:1915059/1‑100 (MQ=255)
                                                                                       gTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc       >  2:2451186/1‑100 (MQ=255)
                                                                                       gTTCATCATTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc       >  2:2410030/1‑100 (MQ=255)
                                                                                          catcatTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa    >  1:830916/1‑100 (MQ=255)
                                                                                            tcatTA‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg  <  1:724287/100‑1 (MQ=255)
                                                                                                 |                                                                                                
CTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAA‑AACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494780‑4494970

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