Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 194,128 C→T 100% P203L (CCG→CTG)  dxr → 1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913194,1280CT94.3% 249.7 / 4.8 87P203L (CCG→CTG) dxr1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (37/45);  total (39/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

TGACCTTGAGCAAAATGGCGTGGTGTCCATTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCCGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAG  >  NC_000913/194033‑194223
                                                                                               |                                                                                               
tGACCTTGAGCAAAATGGCGTGGTGTCCATTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGAt                                                                                             <  1:2420686/100‑1 (MQ=255)
tGACCTTGAGCAAAATGGCGTGGTGTCCATTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGAt                                                                                             <  2:1265367/100‑1 (MQ=255)
tGACCTTGAGCAAAATGGCGTGGTGTCCATTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCCGGAt                                                                                             <  1:1228869/100‑1 (MQ=255)
tGACCTTGAGAAAAATGGCGTGGTGTCCATTCTACTCACCGGGTCTGGCGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGAt                                                                                             <  2:1527698/100‑1 (MQ=255)
    cTTGAGCAAAATGGCGTGGTGTCCATTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAg                                                                                         <  2:1021342/100‑1 (MQ=255)
      tGAGCAAAATGGCGTGGTGTCCATTTTACTTACCGTGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGcc                                                                                       <  2:599846/100‑1 (MQ=255)
          cAAAATGGCGTGGTGTCCATTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGcc                                                                                   >  2:260998/1‑100 (MQ=255)
          cAAAATGGCGTGGTGTCCATTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGcc                                                                                   >  2:1994088/1‑100 (MQ=255)
             aaTGGCGTGGTGTCCATTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTc                                                                                <  1:159048/100‑1 (MQ=255)
               tGGCGTGGTGTCCATTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCAt                                                                              <  1:527177/100‑1 (MQ=255)
                ggacaccacgccATTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATc                                                                             <  2:2360214/90‑1 (MQ=255)
                      gtgtCCATTTTACTTGCCGGGTCTGGTGGCCCTTTCCGTGAGACGCCAGGGCGCGATTTGGCAACAGTGACGCTGGATTAAGCCGGCCGTCATCCGAAc                                                                        <  2:1570851/99‑1 (MQ=255)
                      gtgtCCATTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACt                                                                       <  2:1412383/100‑1 (MQ=255)
                      gtgtCCACTTTACTTACGGGGTCTGGTGGCCCTTTCCGGGAGCCGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCGGCCGTCATCCGAACt                                                                       <  2:749821/100‑1 (MQ=255)
                        gtCCATTTTACTTACCGGGTCTGGTGTCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTgg                                                                     >  1:608147/1‑100 (MQ=255)
                        gtCCATTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTg                                                                      <  1:1694367/99‑1 (MQ=255)
                            aTTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGa                                                                 <  1:1273481/100‑1 (MQ=255)
                                 acttacCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATgggg                                                            <  2:1751610/100‑1 (MQ=255)
                                  cttacCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGc                                                           >  1:110390/1‑100 (MQ=255)
                                  cttacCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGc                                                           >  1:679947/1‑100 (MQ=255)
                                  cttacCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGc                                                           >  1:1239395/1‑100 (MQ=255)
                                   ttacCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCg                                                          <  1:1700860/100‑1 (MQ=255)
                                   ttacCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCg                                                          <  1:64151/100‑1 (MQ=255)
                                    tacCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCCACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGt                                                         >  2:685667/1‑100 (MQ=255)
                                     acCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTa                                                        >  1:409082/1‑100 (MQ=255)
                                       cGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTaaa                                                      >  1:106241/1‑100 (MQ=255)
                                       cGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTaaa                                                      <  1:2267371/100‑1 (MQ=255)
                                             tggtcgCCCCTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCt                                                <  1:511068/100‑1 (MQ=255)
                                              ggttgCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTg                                               <  1:1804384/100‑1 (MQ=255)
                                              ggtggCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTg                                               <  2:1992024/100‑1 (MQ=255)
                                              ggtggCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTg                                               >  2:1978114/1‑100 (MQ=255)
                                              ggtggCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTg                                               >  2:1796730/1‑100 (MQ=255)
                                              ggtggCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTg                                               >  2:1529445/1‑100 (MQ=255)
                                              ggtggCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTg                                               >  1:2276014/1‑100 (MQ=255)
                                               gtggCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTg                                               <  1:144016/99‑1 (MQ=255)
                                               gtggCCATTTCCGTGAGCCGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTg                                               <  2:1528017/99‑1 (MQ=255)
                                                tggCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTc                                             <  1:2247172/100‑1 (MQ=255)
                                                 gtcccTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAAGGACGCTGGAGCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTc                                             <  2:2357597/97‑1 (MQ=255)
                                                  gCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGa                                           >  1:1353686/1‑100 (MQ=255)
                                                   cccTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGAt                                          >  1:721296/1‑100 (MQ=255)
                                                   cccTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGAt                                          <  1:8831/100‑1 (MQ=255)
                                                   cccTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGAt                                          >  2:554603/1‑100 (MQ=255)
                                                   cccTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGAt                                          <  1:1764509/100‑1 (MQ=255)
                                                    ccTTTCCGTGCGACGCCATTGCGCTATTTGTCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGTTCGATGGGGCGTAAAATTTCTGTCg                                            <  2:2118025/97‑1 (MQ=255)
                                                    ccTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGAtt                                         >  2:1717783/1‑100 (MQ=255)
                                                    ccTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGAtt                                         >  2:294885/1‑100 (MQ=255)
                                                      tttCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCg                                       <  2:679940/100‑1 (MQ=255)
                                                       ttCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCgg                                      >  2:972891/1‑100 (MQ=255)
                                                       ttCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCgg                                      <  1:1978128/100‑1 (MQ=255)
                                                        tCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGc                                     <  1:1994078/100‑1 (MQ=255)
                                                         ccGTGAGACGCCATTGCGCGATTTGGCAACTATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTTGATGGGGCGTAAAATTTCTGTCGATTCGGc                                     <  2:1728190/99‑1 (MQ=255)
                                                           gTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGCACTGGTCGACGAGGCGAAAAATTAATGTAGATACGGCttc                                  >  2:2419939/1‑98 (MQ=255)
                                                           gTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTAc                                  >  2:750026/1‑100 (MQ=255)
                                                           gTGAGACGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGCCTAc                                  >  2:767510/1‑100 (MQ=255)
                                                           gTGAGACGCCATTGCGCGATTTGGCAACAATGACGCCGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTAc                                  >  1:2091234/1‑100 (MQ=255)
                                                             gagaCGCCATTGCGCGATTTGGCAACAATGACGCCGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCa                                <  1:2218590/100‑1 (MQ=255)
                                                               gaCGCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTTCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCa                                <  2:240215/98‑1 (MQ=255)
                                                                aCGCTATTGCGCTATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCatga                             <  1:1754826/100‑1 (MQ=255)
                                                                  gCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCatgatg                           <  1:266317/100‑1 (MQ=255)
                                                                  gCCATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCatgatg                           <  1:888303/100‑1 (MQ=255)
                                                                   ccATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCAtgatga                          >  2:1522850/1‑100 (MQ=255)
                                                                   ccATTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTAGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCatgatg                           >  2:1859449/1‑99 (MQ=255)
                                                                     aTTGCGCGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAACTTTCTGTCGATTCGGCTACCATCATGAAc                        >  2:1529296/1‑100 (MQ=255)
                                                                        gcgcGATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACaaa                     >  2:1766350/1‑100 (MQ=255)
                                                                         cgcgATTTGGCAACAAAGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCAGCTACCATGATGAACAAAg                    >  2:1545570/1‑100 (MQ=255)
                                                                          gcgATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAgg                   <  1:1796744/100‑1 (MQ=255)
                                                                           cgATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGt                  >  2:1450545/1‑100 (MQ=255)
                                                                           cgATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGt                  >  1:931637/1‑100 (MQ=255)
                                                                           cgATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCCACCATGATGAACAAAGGt                  >  1:1196516/1‑100 (MQ=255)
                                                                            gATTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTc                 >  1:2186921/1‑100 (MQ=255)
                                                                             aTTTGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTc                 >  2:975625/1‑99 (MQ=255)
                                                                                tGGGAACTAGGACGCTGGTTCACTCCTGCCCTGAGTCGAGCTGGTTGATGGGGCGTAAAGTTTCTGTCGATTCGGCTACCTTGATGAACAAAGGTCTGGa             <  2:182113/100‑1 (MQ=255)
                                                                                tGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGa             >  2:1540571/1‑100 (MQ=255)
                                                                                tGGCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGa             >  1:156964/1‑100 (MQ=255)
                                                                                 ggCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGaa            >  1:1473020/1‑100 (MQ=255)
                                                                                 ggCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGaa            <  1:1994102/100‑1 (MQ=255)
                                                                                 ggCAACAATGACGCCGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGaa            >  1:1466977/1‑100 (MQ=255)
                                                                                  gCAACAATGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAAt           >  2:1906642/1‑100 (MQ=255)
                                                                                      caaTGACGCCGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGTAATACAt       <  1:1342785/100‑1 (MQ=255)
                                                                                       aaTGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACAtt      <  2:197902/100‑1 (MQ=255)
                                                                                       aaTGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACAtt      <  2:721289/100‑1 (MQ=255)
                                                                                       aaTGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACAt       >  1:1951170/1‑99 (MQ=255)
                                                                                         tGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGa    <  1:145635/100‑1 (MQ=255)
                                                                                         tGACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGa    <  2:1473009/100‑1 (MQ=255)
                                                                                          gACGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTg     <  2:524744/98‑1 (MQ=255)
                                                                                           aCGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAg  <  2:520500/100‑1 (MQ=255)
                                                                                           aCGCTGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAg  <  1:1906656/100‑1 (MQ=255)
                                                                                               |                                                                                               
TGACCTTGAGCAAAATGGCGTGGTGTCCATTTTACTTACCGGGTCTGGTGGCCCTTTCCGTGAGACGCCATTGCGCGATTTGGCAACAATGACGCCGGATCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAG  >  NC_000913/194033‑194223

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: