Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,537,076 C→G 17.2% L54L (CTG→CTC pdxK ← pyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,537,0760CG17.2% 313.9 / 58.5 122L54L (CTG→CTCpdxKpyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base C (43/58);  new base G (12/9);  total (55/67)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.32e-01

Only 100 of 127 total aligned reads displayed.

GCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGC  >  NC_000913/2536977‑2537175
                                                                                                   |                                                                                                   
gcgcagcgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTc                                                                                                     <  2:3410798/100‑1 (MQ=255)
gcgcagcgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTc                                                                                                     >  1:3098886/1‑100 (MQ=255)
   cagcgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGc                                                                                                  >  2:2803306/1‑100 (MQ=255)
   cagcgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGc                                                                                                  <  2:632089/100‑1 (MQ=255)
      cgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAAt                                                                                               <  1:135340/100‑1 (MQ=255)
         aTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATAcc                                                                                            <  2:902034/100‑1 (MQ=255)
          tCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgc                                                                                           <  1:3635060/100‑3 (MQ=255)
          tCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgc                                                                                           <  1:489263/100‑3 (MQ=255)
          tCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgc                                                                                           <  2:2291988/100‑3 (MQ=255)
          tCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgc                                                                                           <  2:2166816/100‑3 (MQ=255)
          tCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCg                                                                                           <  2:2741005/100‑1 (MQ=255)
          tCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTAATAATGCGGCGTATTGCTCAGCAATACCg                                                                                           >  1:315442/1‑100 (MQ=255)
           cACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgcc                                                                                          <  1:2810312/100‑4 (MQ=255)
           cACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgcc                                                                                          <  2:1792824/100‑4 (MQ=255)
           cACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgcc                                                                                          <  1:1076456/100‑4 (MQ=255)
            aCGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgccg                                                                                         <  1:3648186/100‑5 (MQ=255)
             cGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgccgc                                                                                        <  1:2565150/100‑6 (MQ=255)
             cGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcg                                                                                        <  1:3637568/100‑1 (MQ=255)
              gCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcgg                                                                                       <  1:2748508/100‑1 (MQ=255)
               ccccTGAAGCGCACGCAAATAGCCGCTAAACCATGCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggc                                                                                      <  1:807111/98‑1 (MQ=255)
                  cTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacg                                                                                   >  2:1025252/1‑100 (MQ=255)
                  cTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacg                                                                                   <  1:2087218/100‑1 (MQ=255)
                  cTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacg                                                                                   <  2:278163/100‑1 (MQ=255)
                   tGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacgg                                                                                  <  2:1632358/100‑1 (MQ=255)
                       gcgcACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCaaa                                                                              <  1:2086929/100‑1 (MQ=255)
                         gcacgcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAg                                                                             <  2:1699109/99‑1 (MQ=255)
                         gcacgcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGa                                                                            >  2:1212974/1‑100 (MQ=255)
                          cacgcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGTACGGCAAAGAc                                                                           >  1:915407/1‑100 (MQ=255)
                            cgcaCATAACCGCTAACCCATGTTTCCGCACTCGCACTACCGTATAAAGCGTCATAATGCGGCGTATTGCTCTGCAAGACCGTCGGCACGGCtaat                                                                             <  2:3383021/96‑2 (MQ=255)
                              caAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTc                                                                       <  2:471368/100‑1 (MQ=255)
                              caAATAGCCGCTAAACCATTCGTCCGGAATCACACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTc                                                                       >  1:947217/1‑100 (MQ=255)
                               aaaTAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCa                                                                      <  1:1017584/100‑1 (MQ=255)
                               aaaTAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCa                                                                      <  2:2811791/100‑1 (MQ=255)
                                  tAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGc                                                                   <  1:1499959/100‑1 (MQ=255)
                                    gCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCg                                                                 <  1:2621408/100‑1 (MQ=255)
                                    gCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCg                                                                 >  2:3738042/1‑100 (MQ=255)
                                     ccGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGCCATTCACGCCGt                                                                >  2:2638545/1‑100 (MQ=255)
                                     ccGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGCGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCCGGCCg                                                                 >  2:2638878/1‑99 (MQ=255)
                                      cGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGt                                                                >  1:3445507/1‑99 (MQ=255)
                                      cGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCACGCCg                                                                 >  2:762045/1‑98 (MQ=255)
                                       gCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTTCGGCCGTTc                                                              >  2:1543143/1‑100 (MQ=255)
                                       gCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTc                                                              <  1:1090548/100‑1 (MQ=255)
                                        cTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCt                                                             >  1:1764094/1‑100 (MQ=255)
                                           aaCCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGtt                                                          <  1:2111047/100‑1 (MQ=255)
                                           aaCCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGtt                                                          >  2:3715759/1‑100 (MQ=255)
                                           aaCCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGtt                                                          >  1:867053/1‑100 (MQ=255)
                                            aCCATTCGTCCGGAATCGCACCACCGTAGAAATTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGtt                                                          <  2:2356493/99‑1 (MQ=255)
                                            aCCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGttt                                                         <  2:2285276/100‑1 (MQ=255)
                                            aCCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGttt                                                         <  2:1070904/100‑1 (MQ=255)
                                            aCCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGttt                                                         <  1:126172/100‑1 (MQ=255)
                                             ccATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTg                                                        >  1:420855/1‑100 (MQ=255)
                                             ccATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTg                                                        >  1:3343365/1‑100 (MQ=255)
                                              cATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGa                                                       <  1:289043/100‑1 (MQ=255)
                                               aTTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGAt                                                      <  2:3027738/100‑1 (MQ=255)
                                               aTTCGACCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGACACCGCAAAGACATTCAGGCCGTTCTGTTTGAt                                                      <  2:1054375/100‑1 (MQ=255)
                                                 tCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAg                                                    <  1:1411793/100‑1 (MQ=255)
                                                   gTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCa                                                  <  2:1979078/100‑1 (MQ=255)
                                                     ccGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAgg                                                >  1:1246326/1‑100 (MQ=255)
                                                      cGGAATCGCACCACTGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGATCTGTTTGATAGCAg                                                 >  2:950009/1‑98 (MQ=255)
                                                        gAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATATCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCAc                                             >  1:1053846/1‑100 (MQ=255)
                                                           tCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGc                                          >  1:240540/1‑100 (MQ=255)
                                                            cGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCa                                         >  1:1784084/1‑100 (MQ=255)
                                                             gCACCACCGTAGAAAGTGTCATAATGCGTCGTATCGCTCAGAAATACCGTCGGCACGGCTTAGACATTCAGGCCGGTCTGTTTGATAGCAGGCACGGCa                                         >  2:1810562/1‑99 (MQ=255)
                                                                  accGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTg                                   >  1:2165395/1‑100 (MQ=255)
                                                                   ccGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGt                                  <  1:3090553/100‑1 (MQ=255)
                                                                        gAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGccc                             >  2:966764/1‑100 (MQ=255)
                                                                        gAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGccc                             >  1:2374466/1‑100 (MQ=255)
                                                                         aaaGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCa                            <  2:1295365/100‑1 (MQ=255)
                                                                             tgtCATAATGCGGCGTATTGCTCAGCGATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCt                        <  2:1064445/100‑1 (MQ=255)
                                                                             tgtCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCt                        >  2:831068/1‑100 (MQ=255)
                                                                              gtCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTg                       >  1:2060891/1‑100 (MQ=255)
                                                                               tCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGc                      >  1:3273205/1‑100 (MQ=255)
                                                                                cATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGcc                     >  1:1717940/1‑100 (MQ=255)
                                                                                cATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGcc                     <  2:1846122/100‑1 (MQ=255)
                                                                                 cgtgccgacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCg                    >  2:2088140/11‑100 (MQ=38)
                                                                                  tAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGt                   <  2:1467871/100‑1 (MQ=255)
                                                                                  tAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGt                   >  2:1210603/1‑100 (MQ=255)
                                                                                  tAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGt                   <  1:729870/100‑1 (MQ=255)
                                                                                   aaTGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTa                  >  2:833236/1‑100 (MQ=255)
                                                                                   aaTGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTa                  <  2:3170107/100‑1 (MQ=255)
                                                                                    gccgacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTaa                 >  1:609924/8‑100 (MQ=255)
                                                                                     tGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTaaa                <  1:795179/100‑1 (MQ=255)
                                                                                     ccgacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTaaa                >  1:2468065/7‑100 (MQ=255)
                                                                                      gcggcgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAac               <  2:3115079/100‑1 (MQ=255)
                                                                                      gcggcgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAac               >  2:3230820/1‑100 (MQ=255)
                                                                                       gacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAacc              >  1:2983749/5‑100 (MQ=255)
                                                                                       gacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAacc              >  1:290823/5‑100 (MQ=255)
                                                                                       gacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAacc              >  2:3484116/5‑100 (MQ=255)
                                                                                       cggcgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAacc              >  2:1275674/1‑100 (MQ=255)
                                                                                       cggcgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAacc              <  1:1351289/100‑1 (MQ=255)
                                                                                         gcgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAaccac            >  2:1792261/1‑100 (MQ=255)
                                                                                         cggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAaccac            >  2:2914692/3‑100 (MQ=255)
                                                                                         cggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTATCATTGATAGCAGGCACGGCAATGCTGATGCCCCCGCTACCGTCAaccac            >  2:106713/3‑100 (MQ=255)
                                                                                              ttGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCg       <  1:307724/100‑1 (MQ=255)
                                                                                              ttGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCg       >  1:1927494/1‑100 (MQ=255)
                                                                                                gCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCGATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGAc     >  1:1834391/1‑100 (MQ=255)
                                                                                                 cTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACt    >  1:521088/1‑100 (MQ=255)
                                                                                                  tCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTATTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTg   <  1:344371/100‑1 (MQ=255)
                                                                                                   cAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGc  >  2:1879319/1‑100 (MQ=255)
                                                                                                   cAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGATCTGTTTGATAGCAGGCACGCCAATGATGTCGTCCACGCTGCGCTAAACCACCTGCGACTGc  >  2:438671/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGC  >  NC_000913/2536977‑2537175

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: