Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,224 A→G I40V (ATC→GTC)  ilvE → branched‑chain amino‑acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,2240AG99.6% 695.8 / ‑4.8 224I40V (ATC→GTC) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (121/102);  total (121/103)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

Only 100 of 229 total aligned reads displayed.

TGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTATCGC  >  NC_000913/3950134‑3950313
                                                                                          |                                                                                         
tGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTc                                                                                         <  2:1135664/93‑1 (MQ=255)
tGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTc                                                                                         <  2:2168384/93‑1 (MQ=255)
    tCAATGGGGAGATGGTTCGCTGGGAATGCGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGtt                                                                                     <  1:1350268/93‑1 (MQ=255)
    tCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGtt                                                                                     >  1:3292360/1‑93 (MQ=255)
    tCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGtt                                                                                     >  2:1515549/1‑93 (MQ=255)
    tCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGtt                                                                                     >  2:568264/1‑93 (MQ=255)
     cAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTg                                                                                    >  1:2097538/1‑93 (MQ=255)
      aaTGGGGAGATGGTTCGCTGGGAAGACGCGGAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGc                                                                                   >  1:4001282/1‑93 (MQ=255)
      aaTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGc                                                                                   >  1:454023/1‑93 (MQ=255)
      aaTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGc                                                                                   >  1:2092782/1‑93 (MQ=255)
       aTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCt                                                                                  >  1:4134187/1‑93 (MQ=255)
       aTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCt                                                                                  >  2:67911/1‑93 (MQ=255)
        tGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTa                                                                                 >  2:4104462/1‑93 (MQ=255)
            gagaTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGAc                                                                             <  1:3405610/93‑1 (MQ=255)
            gagaTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGAc                                                                             <  2:1354018/93‑1 (MQ=255)
            gagaTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGAc                                                                             >  1:3953771/1‑93 (MQ=255)
              gaTGGTTCGCTGGGAAGACGCGAAGGCGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTCCGACTc                                                                           >  2:1883702/1‑93 (MQ=255)
                ttgTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGc                                                                         <  2:776034/91‑1 (MQ=255)
                tGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGc                                                                         <  2:1064970/93‑1 (MQ=255)
                tGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGc                                                                         <  1:2958008/93‑1 (MQ=255)
                tGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGc                                                                         >  1:2820945/1‑93 (MQ=255)
                 ggTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGca                                                                        <  2:4569196/93‑1 (MQ=255)
                  gTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGcac                                                                       >  2:3290572/1‑93 (MQ=255)
                  gTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGcac                                                                       <  1:1971797/93‑1 (MQ=255)
                  gTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTAGcac                                                                       >  2:3171158/1‑93 (MQ=255)
                    tCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACaa                                                                     <  1:4514231/93‑1 (MQ=255)
                        tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGa                                                                 >  2:2692722/1‑93 (MQ=255)
                         gggAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGAc                                                                <  2:1330096/93‑1 (MQ=255)
                             aGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGt                                                            <  2:1542369/93‑1 (MQ=255)
                             aGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGt                                                            <  2:2824830/93‑1 (MQ=255)
                             aGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGt                                                            >  1:343842/1‑93 (MQ=255)
                                 gcgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTa                                                        <  1:4654555/93‑1 (MQ=255)
                                 gcgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGGTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTa                                                        >  1:1146742/1‑93 (MQ=255)
                                 gcgAAGGTGCATGTGATGTCGCACGCGCCGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTa                                                        >  1:1905939/1‑93 (MQ=255)
                                   gAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTAtt                                                      <  2:2791237/93‑1 (MQ=255)
                                   gAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTAtt                                                      >  2:3082225/1‑93 (MQ=255)
                                    aaGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTc                                                     >  1:3752508/1‑93 (MQ=255)
                                    aaGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTc                                                     <  2:2625810/93‑1 (MQ=255)
                                    aaGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTc                                                     >  2:4315182/1‑93 (MQ=255)
                                     aGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTcc                                                    >  2:3243040/1‑93 (MQ=255)
                                     aGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTcc                                                    >  2:4294582/1‑93 (MQ=255)
                                      ggTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCg                                                   >  2:4601125/1‑93 (MQ=255)
                                      ggTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCg                                                   >  1:1295212/1‑93 (MQ=255)
                                        tGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGcc                                                 <  2:185065/93‑1 (MQ=255)
                                        tGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGcc                                                 >  1:3713876/1‑93 (MQ=255)
                                        tGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGcc                                                 <  1:821370/93‑1 (MQ=255)
                                           aTGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATc                                              >  2:1936556/1‑93 (MQ=255)
                                           aTGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATc                                              >  2:3587905/1‑93 (MQ=255)
                                            tgtgATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCg                                             >  2:3381988/1‑93 (MQ=255)
                                            tgtgATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCg                                             <  1:1205389/93‑1 (MQ=255)
                                            tgtgATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCg                                             >  1:3885853/1‑93 (MQ=255)
                                              tgATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTg                                           >  1:1111714/1‑93 (MQ=255)
                                               gATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGa                                          >  1:3013030/1‑93 (MQ=255)
                                               gATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGa                                          <  1:641929/93‑1 (MQ=255)
                                                aTGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTTCGCCATCGTGac                                         >  2:1070861/1‑92 (MQ=255)
                                                aTGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAg                                         <  2:4096263/93‑1 (MQ=255)
                                                aTGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAg                                         >  1:4285682/1‑93 (MQ=255)
                                                    cGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCata                                     >  2:1769198/1‑93 (MQ=255)
                                                    cGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCata                                     <  1:230228/93‑1 (MQ=255)
                                                      cACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATg                                   <  1:536911/93‑1 (MQ=255)
                                                      cACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATg                                   >  1:772345/1‑93 (MQ=255)
                                                       aCGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGc                                  >  2:1437631/1‑93 (MQ=255)
                                                       aCGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGc                                  <  1:4346782/93‑1 (MQ=255)
                                                       aCGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGc                                  <  1:1778390/93‑1 (MQ=255)
                                                        cgcgCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCa                                 >  1:885103/1‑93 (MQ=255)
                                                        cgcgCTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCa                                 <  2:2490096/93‑1 (MQ=255)
                                                         gcgcTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAg                                >  1:3084668/1‑93 (MQ=255)
                                                          cgcTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGc                               >  2:3123479/1‑93 (MQ=255)
                                                           gcTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg                              >  2:372575/1‑93 (MQ=255)
                                                            cTGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGt                             >  2:2360999/1‑93 (MQ=255)
                                                             tGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTc                            <  1:2608725/93‑1 (MQ=255)
                                                             tGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTc                            <  1:3840419/93‑1 (MQ=255)
                                                             tGCACTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTc                            <  2:1284579/93‑1 (MQ=255)
                                                                 cTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCa                        >  2:2265968/1‑93 (MQ=255)
                                                                 cTATGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGGTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCa                        >  2:1015733/1‑93 (MQ=255)
                                                                    tGGCACTTCTGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGCCCGGTTGTATTCCGCAATCGTGAGCATATGCAGCGTCTGCATGa                     <  2:3321374/93‑1 (MQ=255)
                                                                    tGGCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGa                     >  2:4369336/1‑93 (MQ=255)
                                                                     ggCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGAc                    <  2:390202/93‑1 (MQ=255)
                                                                     ggCACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGAc                    <  1:3360353/93‑1 (MQ=255)
                                                                       cACTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTc                  >  1:1767282/1‑93 (MQ=255)
                                                                        ccTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTcc                 <  2:3504353/92‑1 (MQ=255)
                                                                        aCTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTcc                 <  1:2380980/93‑1 (MQ=255)
                                                                        aCTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTcc                 <  1:947418/93‑1 (MQ=255)
                                                                        aCTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGCATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTcc                 >  2:2916772/1‑93 (MQ=255)
                                                                         cTTCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCg                >  1:1600515/1‑93 (MQ=255)
                                                                          ttCGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGc               >  1:1000162/1‑93 (MQ=255)
                                                                            cGGTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCa             >  1:4367944/1‑93 (MQ=255)
                                                                              gTTTTTGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCaca           >  2:453290/1‑91 (MQ=255)
                                                                               tttttGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCaaaa          >  1:3776091/1‑93 (MQ=255)
                                                                               tttttGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCaaaa          >  2:4316749/1‑93 (MQ=255)
                                                                                ttttGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAAt         >  2:3033061/1‑93 (MQ=255)
                                                                                ttttGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAAt         <  1:4232864/93‑1 (MQ=255)
                                                                                 tttGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATc        >  1:331697/1‑93 (MQ=255)
                                                                                 tttGAAGGCGTCCGTTGCTACGACTCGCACAAAGGAACGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATc        >  2:2437314/1‑93 (MQ=255)
                                                                                   tGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAACATCTa      >  1:722535/1‑93 (MQ=255)
                                                                                   tGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTa      <  1:3487539/93‑1 (MQ=255)
                                                                                   tGAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTa      >  1:1819514/1‑93 (MQ=255)
                                                                                    gAAGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTAt     <  1:3796402/93‑1 (MQ=255)
                                                                                     aaGGCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTATc    <  1:256888/93‑1 (MQ=255)
                                                                                       ggCGTCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTATCGc  <  1:4090016/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCATGACTCCGCCAAAATCTATCGC  >  NC_000913/3950134‑3950313

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