Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,790,919 A→C L143V (TTA→GTA)  tdh ← L‑threonine 3‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,790,9190AC98.3% 602.9 / 0.3 177L143V (TTA→GTA) tdhL‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base C (92/82);  total (92/85)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

Only 100 of 177 total aligned reads displayed.

GGGCCTGCACCAGAAACCAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGCCTAT  >  NC_000913/3790825‑3791015
                                                                                              |                                                                                                
gggCCTGCACCAGAAACCAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCAt                                                                                             >  2:2434219/1‑100 (MQ=255)
  gCCTGCACCAGAAACCAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCg                                                                                           >  1:1066483/1‑100 (MQ=255)
     tGCACCAGAAACCGGCACATCTTCGCCCACCAGATCAAACGACCGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGTaa                                                                                        >  2:799021/1‑100 (MQ=255)
     tGCACCAGAAACCAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGaa                                                                                        >  2:3181091/1‑100 (MQ=255)
     tGCACCAGAAACCAGCACATCTTCGCCCACCAGATCAAACGACAGCGCCGTCTGCACGTCGTTACCGAAGGGATCACAAATCGCTGCTACGTCATCgtac                                                                                        >  2:1658959/1‑97 (MQ=255)
      gCACCAGAAACCAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGaaa                                                                                       <  1:237071/100‑1 (MQ=255)
           aGAAACCAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTg                                                                                  >  2:3450844/1‑100 (MQ=255)
               aCCAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCgg                                                                              >  2:1372620/1‑100 (MQ=255)
                ccAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCggg                                                                             >  2:1081698/1‑100 (MQ=255)
                ccAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCggg                                                                             <  1:3436732/100‑1 (MQ=255)
                ccAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCggg                                                                             >  1:2227984/1‑100 (MQ=255)
                ccAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCggg                                                                             >  1:2668151/1‑100 (MQ=255)
                 cAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCgggg                                                                            <  2:1050593/100‑1 (MQ=255)
                 cAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCgggg                                                                            >  2:1467168/1‑100 (MQ=255)
                  aGCACATCTTCGCGCACCAGATCATACGACAGCGCGGTATGCGCGGCGTTACCAAATTGAACAAAACTCGCGGCTACGTCAGCGGAAACACTGTCgggta                                                                           >  2:2771367/1‑98 (MQ=255)
                      caTCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGAtttt                                                                       <  2:409225/100‑1 (MQ=255)
                      caTCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGAtttt                                                                       >  1:920374/1‑100 (MQ=255)
                      caTCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGActtt                                                                       >  2:3313869/1‑100 (MQ=255)
                      caTCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGAtttt                                                                       <  2:3004237/100‑1 (MQ=255)
                       aTCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTg                                                                      <  2:528253/100‑1 (MQ=255)
                        tCTTCGCTCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGCTCAAAAATTGCGGCTACGTCATCGGAAATATTGTTGGGGAttttca                                                                     <  2:3710805/100‑3 (MQ=255)
                        tCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGa                                                                     >  1:1187030/1‑100 (MQ=255)
                         cTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGACGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGaa                                                                    >  2:2546250/1‑100 (MQ=255)
                          ttCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAg                                                                   <  1:2392405/100‑1 (MQ=255)
                          ttCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAg                                                                   <  1:856255/100‑1 (MQ=255)
                            cGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGc                                                                 <  2:2228181/100‑1 (MQ=255)
                             gCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAGTTGCGGCTACGCCATCGGAAATATTGCCGGGGATTTCGAAGGCg                                                                >  1:3749063/1‑100 (MQ=255)
                             gCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCg                                                                <  1:1892456/100‑1 (MQ=255)
                              cccACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGt                                                               >  1:3996635/1‑100 (MQ=255)
                               ccaccaGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGtt                                                              >  1:1881635/1‑100 (MQ=255)
                                caccaGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTg                                                             <  1:1947796/100‑1 (MQ=255)
                                 accaGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGa                                                            >  1:4041357/1‑100 (MQ=255)
                                  ccaGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGCCGGGGATTTTGAAGGCGTTGaa                                                           >  2:2142946/1‑100 (MQ=255)
                                  ccaGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGCCTACGTCATCGGACATATTGTCGGGGATTTTGAAGGCGTTGaa                                                           >  2:2844797/1‑100 (MQ=255)
                                  ccaGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGCCTACGTCATCGGAAATATTGTCGGGGCTTTTGCAGGCGTTGaa                                                           >  2:2143311/1‑100 (MQ=255)
                                     gATCAAACGACAGCGCGGTATGCACGGCGTTACCGCAGGGATCAAAAATTGCGGCTCCGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGc                                                        >  2:459657/1‑100 (MQ=255)
                                     gATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGTGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGc                                                        <  1:2576708/100‑1 (MQ=255)
                                      aTCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGcc                                                       >  1:3216435/1‑100 (MQ=255)
                                        cAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCgg                                                     <  1:3296673/100‑1 (MQ=255)
                                        cAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCgg                                                     >  1:2276524/1‑100 (MQ=255)
                                         aaaCGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCggg                                                    <  2:2833925/100‑1 (MQ=255)
                                          aaCGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGa                                                   <  1:2092845/100‑1 (MQ=255)
                                          aaCGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGa                                                   <  2:2850684/100‑1 (MQ=255)
                                            cGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATc                                                 <  1:1727947/100‑1 (MQ=255)
                                            cGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATc                                                 <  1:633053/100‑1 (MQ=255)
                                             gACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCa                                                <  1:408459/100‑1 (MQ=255)
                                             gACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCa                                                <  2:3118646/100‑1 (MQ=255)
                                              aCAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCAc                                               <  2:2645119/100‑1 (MQ=255)
                                              aCAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCAc                                               <  1:1030151/100‑1 (MQ=255)
                                               cAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCAcc                                              >  1:2137900/1‑100 (MQ=255)
                                                aGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCa                                             <  1:912010/100‑1 (MQ=255)
                                                 gcgcGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAg                                            >  1:1806246/1‑100 (MQ=255)
                                                 gcgcGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAg                                            <  1:3385494/100‑1 (MQ=255)
                                                 gcgcGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAg                                            <  1:1276071/100‑1 (MQ=255)
                                                  cgcgGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGa                                           >  1:487511/1‑100 (MQ=255)
                                                  cgcgGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGa                                           >  1:308504/1‑100 (MQ=255)
                                                  cgcgGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGa                                           <  2:1376904/100‑1 (MQ=255)
                                                        aTGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCg                                     >  1:2431172/1‑100 (MQ=255)
                                                          gCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGc                                   >  2:1731290/1‑100 (MQ=255)
                                                          gCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGc                                   <  2:3512593/100‑1 (MQ=255)
                                                              ggCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAag                               <  1:3843805/100‑1 (MQ=255)
                                                               gCGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAagn                              <  1:3255634/100‑2 (MQ=255)
                                                                cGTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAagca                             >  1:586400/1‑100 (MQ=255)
                                                                 ggtACCTAAGGGATCAAAAATTGCGGCTACCTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAagcag                            <  2:3342338/98‑1 (MQ=255)
                                                                 gTTACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATCCCGGGATCACCAGATATTCGGCAAagcag                            >  2:3649940/1‑100 (MQ=255)
                                                                   tACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGcc                          <  1:1096062/100‑1 (MQ=255)
                                                                   tACCGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGcc                          >  2:1737512/1‑100 (MQ=255)
                                                                     ccGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCg                        >  2:945447/1‑100 (MQ=255)
                                                                      cGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCgg                       >  1:3707864/1‑100 (MQ=255)
                                                                      cGAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCgg                       >  1:332970/1‑100 (MQ=255)
                                                                       gAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCggg                      <  1:2109678/100‑1 (MQ=255)
                                                                       gAAGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCggg                      >  1:3568528/1‑100 (MQ=255)
                                                                        aaGGGATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATACCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGc                     >  1:1365215/1‑100 (MQ=255)
                                                                          gggATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGa                   >  2:1921825/1‑100 (MQ=255)
                                                                           ggATCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGAt                  <  1:122095/100‑1 (MQ=255)
                                                                           gaataaaaaaTTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGAt                  <  1:589763/95‑1 (MQ=255)
                                                                             ttCAAGAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTa                <  2:1794979/99‑1 (MQ=255)
                                                                             aTCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAACGCAGCCCGGGCGATTa                >  2:1519348/1‑100 (MQ=255)
                                                                             aTCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTa                <  1:1017941/100‑1 (MQ=255)
                                                                             aTCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTa                <  1:3183470/100‑1 (MQ=255)
                                                                             aTCAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTa                >  1:665276/1‑100 (MQ=255)
                                                                               cAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAac              <  2:716536/100‑1 (MQ=255)
                                                                               cAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAac              >  2:619410/1‑100 (MQ=255)
                                                                               cAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAac              <  2:3266169/100‑1 (MQ=255)
                                                                               cAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAac              >  2:2321194/1‑100 (MQ=255)
                                                                               cAAAAATTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAac              <  1:1502977/100‑1 (MQ=255)
                                                                                aaaaaTTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAaca             >  2:3032738/1‑100 (MQ=255)
                                                                                aaaaaTTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAaca             >  1:775428/1‑100 (MQ=255)
                                                                                 aaaaTTGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAacac            >  1:3340978/1‑100 (MQ=255)
                                                                                 aaaaTTGCGGCTACGTCATCGGAAATATTGTAGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAacac            <  1:386492/100‑1 (MQ=255)
                                                                                 aaaaTTGCGGCTACGTCATCGGAAAGATTGTCGGGGGTTTTGAAGGCGTTGAATGCCGGGATCGGCAGATATTCGGCGAAGCAGCCAGGGCGAGTAacac            >  1:3672903/1‑100 (MQ=255)
                                                                                 aaaaTTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAacac            <  1:1180598/100‑1 (MQ=255)
                                                                                     ttGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCAGGGCGATTAACACCAAc        >  2:3502973/1‑100 (MQ=255)
                                                                                      tGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACg       <  1:2107884/100‑1 (MQ=255)
                                                                                      tGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACg       <  1:155397/100‑1 (MQ=255)
                                                                                      tGCGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACg       <  1:2482940/100‑1 (MQ=255)
                                                                                        cGGCTGCGTCCTCGGAAATATGGTCGGGGATTTGGAAGGCGTTGAATGCCGGGATCACCAGATATTCTGCAAAGCAGCCCGGGCGCTTACCACCAAAGcc     >  2:3745917/1‑100 (MQ=255)
                                                                                        cGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGcc     <  1:2243505/100‑1 (MQ=255)
                                                                                        cGGCTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGcc     >  1:3523194/1‑100 (MQ=255)
                                                                                           cTACGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGCCTAt  >  1:653378/1‑100 (MQ=255)
                                                                                              |                                                                                                
GGGCCTGCACCAGAAACCAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGCCTAT  >  NC_000913/3790825‑3791015

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: ‑