Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,888,375 A→G L516P (CTT→CCT)  cysJ ← sulfite reductase (NADPH)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,888,3750AG86.9% 671.9 / 90.7 274L516P (CTT→CCT) cysJsulfite reductase (NADPH)
Reads supporting (aligned to +/- strand):  ref base A (18/18);  new base G (127/111);  total (145/129)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

Only 100 of 276 total aligned reads displayed.

CCATCATTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAG  >  NC_000913/2888285‑2888464
                                                                                          |                                                                                         
ccATCATTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGa                                                                                         >  2:859999/1‑93 (MQ=255)
ccATCATTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGa                                                                                         >  2:4466207/1‑93 (MQ=255)
ccATCATTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGa                                                                                         >  2:1105646/1‑93 (MQ=255)
ccATCATTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGa                                                                                         <  2:4291683/93‑1 (MQ=255)
 catcatTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGAt                                                                                        <  2:702758/93‑1 (MQ=255)
   tcatTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCg                                                                                      >  2:1624821/1‑93 (MQ=255)
    catTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGa                                                                                     >  2:5584427/1‑93 (MQ=255)
    catTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGa                                                                                     <  1:774035/93‑1 (MQ=255)
      ttGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATa                                                                                   <  2:2151388/93‑1 (MQ=255)
       tGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATAc                                                                                  <  1:301555/93‑1 (MQ=255)
         aTCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACgt                                                                                >  1:6915344/1‑93 (MQ=255)
         aTCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACgt                                                                                >  1:1700184/1‑93 (MQ=255)
          tCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACgtg                                                                               >  2:2705177/1‑93 (MQ=255)
          tCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACgtg                                                                               <  1:7507808/93‑1 (MQ=255)
           ccAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACgtgt                                                                              <  2:5751256/93‑1 (MQ=255)
           ccAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACgtgt                                                                              <  2:3444337/93‑1 (MQ=255)
            cAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTc                                                                             <  1:1619958/93‑1 (MQ=255)
            cAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTc                                                                             >  1:4648493/1‑93 (MQ=255)
             aGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCa                                                                            <  2:1891130/93‑1 (MQ=255)
             aGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCa                                                                            <  1:6418152/93‑1 (MQ=255)
              gcgcCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAg                                                                           >  1:866739/1‑93 (MQ=255)
              gcgcCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAg                                                                           >  1:1023030/1‑93 (MQ=255)
                gcCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCa                                                                         <  2:4679453/93‑1 (MQ=255)
                gcCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCa                                                                         >  2:758721/1‑93 (MQ=255)
                   acaGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGACACGTGTCAGCACGc                                                                      >  2:3549451/1‑93 (MQ=255)
                    caGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGcc                                                                     >  1:3484800/1‑93 (MQ=255)
                    caGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGcc                                                                     <  1:5759696/93‑1 (MQ=255)
                     aGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCa                                                                    <  2:2088942/93‑1 (MQ=255)
                     aGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCa                                                                    <  2:1326393/93‑1 (MQ=255)
                     aGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCa                                                                    >  1:2817974/1‑93 (MQ=255)
                      gCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCAt                                                                   >  1:5032643/1‑93 (MQ=255)
                      gCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCAt                                                                   <  2:5176824/93‑1 (MQ=255)
                       cTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATc                                                                  <  1:3567247/93‑1 (MQ=255)
                        tCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCt                                                                 >  2:297298/1‑93 (MQ=255)
                        tCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCt                                                                 <  1:3073801/93‑1 (MQ=255)
                         ccGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCtt                                                                <  2:6139217/93‑1 (MQ=255)
                         ccGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCtt                                                                <  2:3222125/93‑1 (MQ=255)
                          cgcgCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCttt                                                               <  2:586113/93‑1 (MQ=255)
                            cgcCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGa                                                             <  1:6826213/93‑1 (MQ=255)
                             gcCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGAc                                                            <  1:6937555/93‑1 (MQ=255)
                              cccTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACg                                                           >  1:4983362/1‑93 (MQ=255)
                               ccTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGt                                                          >  1:7400798/1‑93 (MQ=255)
                                cTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTa                                                         >  1:532551/1‑93 (MQ=255)
                                  gTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAgc                                                       >  1:3370074/1‑93 (MQ=255)
                                  gTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAgc                                                       >  1:6089768/1‑93 (MQ=255)
                                          aGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCAc                                               >  2:149250/1‑93 (MQ=255)
                                          aGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCAc                                               >  2:439469/1‑93 (MQ=255)
                                          aGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCAc                                               >  1:4272/1‑93 (MQ=255)
                                          aGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCAc                                               <  2:1885220/93‑1 (MQ=255)
                                           gTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACt                                              >  1:6464617/1‑93 (MQ=255)
                                            tttGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTc                                             >  1:5670076/1‑93 (MQ=255)
                                              tGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCa                                           <  2:5433585/93‑1 (MQ=255)
                                               gTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCAc                                          >  1:947544/1‑93 (MQ=255)
                                               gTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCAc                                          <  1:2041853/93‑1 (MQ=255)
                                                  ttGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTg                                       <  1:7152111/93‑1 (MQ=255)
                                                  ttGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTg                                       >  1:2260878/1‑93 (MQ=255)
                                                  ttGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTg                                       <  1:5703640/93‑1 (MQ=255)
                                                   tGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTgg                                      <  1:4554040/93‑1 (MQ=255)
                                                     tACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTa                                    <  1:3724848/93‑1 (MQ=255)
                                                      aCGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTAc                                   >  1:2634096/1‑93 (MQ=255)
                                                        gTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAg                                 >  2:4540561/1‑93 (MQ=255)
                                                        gTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAg                                 >  1:7176930/1‑93 (MQ=255)
                                                         tAAACTTTTTCTTTTTGATCGCGCGACCGGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAgg                                >  2:472761/1‑93 (MQ=255)
                                                         tAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAgg                                <  2:2415937/93‑1 (MQ=255)
                                                          aaaCTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGa                               <  2:6951617/93‑1 (MQ=255)
                                                           aaCTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGaa                              >  2:6691023/1‑93 (MQ=255)
                                                           aaCTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGaa                              <  1:1386072/93‑1 (MQ=255)
                                                           aaCTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGaa                              <  1:2601586/93‑1 (MQ=255)
                                                            aCTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAg                             <  1:4823809/93‑1 (MQ=255)
                                                             ctttttctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAgc                            >  2:7278563/1‑92 (MQ=255)
                                                              tttttctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTc                           >  1:3846565/1‑93 (MQ=255)
                                                                tttctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCtt                         <  2:598159/93‑1 (MQ=255)
                                                                 ttctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTc                        >  1:1142986/1‑93 (MQ=255)
                                                                 ttctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTc                        <  1:1553099/93‑1 (MQ=255)
                                                                 ttctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTc                        >  2:238534/1‑93 (MQ=255)
                                                                 ttctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTc                        <  1:4396037/93‑1 (MQ=255)
                                                                 ttctttttGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTc                        <  1:3819174/93‑1 (MQ=255)
                                                                  tctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTcc                       >  2:7449276/1‑93 (MQ=255)
                                                                    tttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGt                     >  1:2110010/1‑93 (MQ=255)
                                                                    tttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGt                     <  1:1944264/93‑1 (MQ=255)
                                                                    tttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTCACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGt                     >  1:236567/1‑93 (MQ=255)
                                                                    tttttGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGt                     >  2:4983462/1‑93 (MQ=255)
                                                                      tttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTaa                   >  1:7454406/1‑93 (MQ=255)
                                                                        tGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAg                 >  1:6609219/1‑93 (MQ=255)
                                                                        tGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAg                 <  2:7100010/93‑1 (MQ=255)
                                                                         gATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGt                <  1:407808/93‑1 (MQ=255)
                                                                          aTCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTg               <  2:374336/93‑1 (MQ=255)
                                                                           tCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGc              <  2:2623894/93‑1 (MQ=255)
                                                                            cgcgcgACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCg             <  2:3898690/93‑1 (MQ=255)
                                                                              cgcgACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGa           <  2:6734233/93‑1 (MQ=255)
                                                                               gcgACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGAt          >  2:152649/1‑93 (MQ=255)
                                                                               gcgACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGCAAAGTGCGGAt          >  1:1429573/1‑93 (MQ=255)
                                                                                 gACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTa        >  2:5309233/1‑93 (MQ=255)
                                                                                 gACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTa        >  1:5997940/1‑93 (MQ=255)
                                                                                  aCCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTAc       <  2:1439567/93‑1 (MQ=255)
                                                                                    caggcaggATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCa     >  2:4106717/1‑93 (MQ=255)
                                                                                     aGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCaa    >  2:179228/1‑93 (MQ=255)
                                                                                      ggCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCaaa   >  2:6525884/1‑93 (MQ=255)
                                                                                      ggCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCaaa   >  2:6366851/1‑93 (MQ=255)
                                                                                      ggCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCaaa   <  2:5248876/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCATCATTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAG  >  NC_000913/2888285‑2888464

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