Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,537,076 C→G 16.3% 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,0760CG16.3% 375.1 / 59.0 141L54L (CTG→CTCpdxKpyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base C (69/48);  new base G (14/9);  total (84/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 145 total aligned reads displayed.

GCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGC  >  NC_000913/2536977‑2537175
                                                                                                   |                                                                                                   
gcgcagcgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTc                                                                                                     >  2:688511/1‑100 (MQ=255)
  gcgtcgcaTTACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCa                                                                                                    <  2:1727873/95‑1 (MQ=255)
  gcagcgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAg                                                                                                   <  1:1995116/100‑1 (MQ=255)
  gcagcgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAg                                                                                                   <  2:2873656/100‑1 (MQ=255)
       gcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATa                                                                                              <  1:835057/100‑1 (MQ=255)
        caTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATAc                                                                                             <  2:2521750/100‑1 (MQ=255)
         aTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATAcc                                                                                            <  1:924086/100‑1 (MQ=255)
           cACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgcc                                                                                          <  2:1901854/100‑4 (MQ=255)
           cACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgcc                                                                                          <  1:1593862/100‑4 (MQ=255)
           cACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgcc                                                                                          <  1:402180/100‑4 (MQ=255)
           cACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGACCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgcc                                                                                          <  2:1071659/100‑4 (MQ=255)
            aCGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTc                                                                                         >  2:1226375/1‑100 (MQ=255)
            aCGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTc                                                                                         <  1:2572511/100‑1 (MQ=255)
            aCGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGAGTCATAATGCGGCGTATTGCTCAGCAATACCGTc                                                                                         >  2:334615/1‑100 (MQ=255)
             cGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcg                                                                                        <  2:2201670/100‑1 (MQ=255)
             cGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcg                                                                                        >  1:3187528/1‑100 (MQ=255)
                  cTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacg                                                                                   >  1:2988426/1‑100 (MQ=255)
                  cTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacg                                                                                   >  2:1372471/1‑100 (MQ=255)
                   tGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacgg                                                                                  >  2:2740759/1‑100 (MQ=255)
                    gAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacggn                                                                                 <  1:68013/100‑2 (MQ=255)
                    gAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacggc                                                                                 <  1:2998794/100‑1 (MQ=255)
                     aaGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATACCGTcggcacggc                                                                                 <  2:2218293/99‑1 (MQ=255)
                     aaGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacggca                                                                                <  1:1035803/100‑1 (MQ=255)
                      aGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAACCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacggc                                                                                 >  2:1946503/1‑98 (MQ=255)
                       gcgcACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCaaa                                                                              >  1:3369881/1‑100 (MQ=255)
                        cgcacgcaAATAGTCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACTGTCGGCACGGCAAAg                                                                             <  2:629446/100‑1 (MQ=255)
                          cacgcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGAc                                                                           >  2:3469710/1‑100 (MQ=255)
                            cgcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACAt                                                                         >  2:2439887/1‑100 (MQ=255)
                            cgcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACAt                                                                         >  1:2896653/1‑100 (MQ=255)
                            cgcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACAt                                                                         >  1:241345/1‑100 (MQ=255)
                              caAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTc                                                                       <  1:687338/100‑1 (MQ=255)
                              caAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTc                                                                       >  1:3437736/1‑100 (MQ=255)
                              caAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTc                                                                       >  1:435706/1‑100 (MQ=255)
                              caAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTc                                                                       >  2:2617034/1‑100 (MQ=255)
                              caAATAGCCGCTAAACCATTCGTCCGGAATCACACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTc                                                                       >  2:1865834/1‑100 (MQ=255)
                                aaTAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAg                                                                     <  1:3247885/100‑1 (MQ=255)
                                 aTAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAgg                                                                    <  1:1226383/100‑1 (MQ=255)
                                 aTAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAgg                                                                    <  1:382039/100‑1 (MQ=255)
                                  tAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGc                                                                   <  1:537975/100‑1 (MQ=255)
                                   aGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGcc                                                                  >  1:3316916/1‑100 (MQ=255)
                                    gCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCg                                                                 >  1:2123411/1‑100 (MQ=255)
                                     ccGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCg                                                                 <  2:2530164/99‑1 (MQ=255)
                                     ccGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGt                                                                >  2:3608977/1‑100 (MQ=255)
                                      cGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGtt                                                               >  1:2542537/1‑100 (MQ=255)
                                      cGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGtt                                                               >  2:3375185/1‑100 (MQ=255)
                                       gCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTc                                                              >  1:1177769/1‑100 (MQ=255)
                                       gCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTc                                                              <  1:1479924/100‑1 (MQ=255)
                                           aaCCATTCGTCCGGAATAGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGtt                                                          <  1:3275571/100‑1 (MQ=255)
                                            aCCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGttt                                                         <  1:2175109/100‑1 (MQ=255)
                                              cATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTg                                                        <  2:3015237/99‑1 (MQ=255)
                                               aTTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGAt                                                      >  2:1999739/1‑100 (MQ=255)
                                                   gTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCa                                                  >  1:3397525/1‑100 (MQ=255)
                                                   gTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACACTCAGGCCGTTCTGTTTGATAGc                                                   >  2:1447820/1‑99 (MQ=255)
                                                    tCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAg                                                 <  2:755331/100‑1 (MQ=255)
                                                     ccGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAgg                                                >  2:1090953/1‑100 (MQ=255)
                                                     ccGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAgg                                                >  1:244515/1‑100 (MQ=255)
                                                      cGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCacg                                                >  2:1678057/1‑97 (MQ=255)
                                                      cGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGc                                               <  2:496650/100‑1 (MQ=255)
                                                       ggAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCa                                              >  1:2067716/1‑100 (MQ=255)
                                                       ggAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCa                                              >  1:519898/1‑100 (MQ=255)
                                                          aTCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACgg                                           <  2:3391408/100‑1 (MQ=255)
                                                            cGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCACGCACGGc                                          >  2:2088309/1‑99 (MQ=255)
                                                              caccacCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAAt                                       >  2:484234/1‑100 (MQ=255)
                                                                ccaccGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCTATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGc                                     <  2:1646894/100‑1 (MQ=255)
                                                                ccaccGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGc                                     <  2:1540459/100‑1 (MQ=255)
                                                                     gTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTg                                >  1:375751/1‑100 (MQ=255)
                                                                        gAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGccc                             >  1:2697998/1‑100 (MQ=255)
                                                                        gAAAGTGTCATAATGCGGCGTATTGCGCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTTTTTGATAGCAGGCACGGCAATGCTGTTGcc                              <  1:3039323/99‑1 (MQ=255)
                                                                            gtgtCATAATGCGGCGTATTACTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACg                          >  2:2231678/1‑99 (MQ=255)
                                                                              gtCATAATGCGGCGTATTGCTCCGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCACGCACGGCAATGCTGTTGCCCACGCTg                       >  2:125532/1‑100 (MQ=255)
                                                                              gtCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTg                       >  2:2635944/1‑100 (MQ=255)
                                                                              gtCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTg                       >  2:1873460/1‑100 (MQ=255)
                                                                                cATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACTGCAAAGACCTTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGc                      >  2:2020139/1‑99 (MQ=255)
                                                                                   aaTGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTa                  >  2:3107977/1‑100 (MQ=255)
                                                                                    aTGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTaa                 <  1:569011/100‑1 (MQ=255)
                                                                                    aTGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTaa                 <  1:1149278/100‑1 (MQ=255)
                                                                                      gcggcgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAac               >  1:3363424/1‑100 (MQ=255)
                                                                                      cgacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAac               >  1:225131/6‑100 (MQ=255)
                                                                                      cgacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAac               >  2:1125424/6‑100 (MQ=255)
                                                                                      cgacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAac               >  1:865011/6‑100 (MQ=255)
                                                                                       gacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAacc              >  1:2798128/5‑100 (MQ=255)
                                                                                       gacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAacc              >  1:109190/5‑100 (MQ=255)
                                                                                       cggcgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAacc              >  2:2268930/1‑100 (MQ=255)
                                                                                        ggcgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAacca             >  1:1626094/1‑100 (MQ=255)
                                                                                        acggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAacca             >  1:471582/4‑100 (MQ=255)
                                                                                         gcgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAaccac            >  2:889708/1‑100 (MQ=255)
                                                                                         cggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAaccac            >  2:2312934/3‑100 (MQ=255)
                                                                                         cggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAaccac            >  2:1727358/3‑100 (MQ=255)
                                                                                         cggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAaccac            >  1:953687/3‑100 (MQ=255)
                                                                                          cgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAaccacc           <  2:1097312/100‑1 (MQ=255)
                                                                                            tATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTg         <  1:920827/100‑1 (MQ=255)
                                                                                             aTTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGc        >  1:2847804/1‑100 (MQ=255)
                                                                                             aTTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGc        >  1:1931832/1‑100 (MQ=255)
                                                                                               tGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGa      >  1:2738398/1‑100 (MQ=255)
                                                                                                gCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGAc     >  1:1212208/1‑100 (MQ=255)
                                                                                                 cTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTTCCCACGCTGCCGTAAACCACCTGCGACt    >  2:2176486/1‑100 (MQ=255)
                                                                                                 cTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACt    >  1:1615803/1‑100 (MQ=255)
                                                                                                  tCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTg   >  1:351204/1‑100 (MQ=255)
                                                                                                   cAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGc  >  1:1306425/1‑100 (MQ=255)
                                                                                                   cAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGc  >  1:1078164/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: