Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,945,351 +C coding (1/600 nt)
coding (240/240 nt)
yifDA ←
yifDA ←
transcription regulator
transcription regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,945,3501.C99.6% 743.2 / ‑4.2 233coding (2/600 nt)yifDAtranscription regulator
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (122/110);  total (122/111)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 252 total aligned reads displayed.

AAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCT  >  NC_000913/3945260‑3945437
                                                                                            |                                                                                       
aaGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACctcc                                                                                       >  2:1790386/1‑92 (MQ=255)
aaGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACctcc                                                                                       >  2:2268273/1‑92 (MQ=255)
  gCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACctcctt                                                                                     >  2:3969978/1‑90 (MQ=255)
     tCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGGATGCGCCACCTCCTTACg                                                                                   >  1:3218273/1‑93 (MQ=255)
      cccAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACgg                                                                                  <  1:2475366/93‑1 (MQ=255)
      cccAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACgg                                                                                  >  2:1512873/1‑93 (MQ=255)
       ccAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACggg                                                                                 >  1:1819268/1‑93 (MQ=255)
       ccAGCCACTGATTAAGCATACATTCACACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGCTATGCGCCACCTCCTTACggc                                                                                 >  2:4357622/1‑92 (MQ=255)
        cAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGggc                                                                                >  1:2387194/1‑93 (MQ=255)
         aGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGggcg                                                                               >  2:2919992/1‑93 (MQ=255)
         aGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGggcg                                                                               >  2:2004069/1‑93 (MQ=255)
          gCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGggcgg                                                                              >  2:752908/1‑93 (MQ=255)
             aCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCt                                                                           >  1:3595407/1‑93 (MQ=255)
             aCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCt                                                                           >  2:3846807/1‑93 (MQ=255)
             aCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCt                                                                           >  1:1499139/1‑93 (MQ=255)
               tGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGc                                                                         >  2:1989987/1‑93 (MQ=255)
               tGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGc                                                                         >  2:1178494/1‑93 (MQ=255)
               tGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGc                                                                         >  1:2620604/1‑93 (MQ=255)
               tGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGc                                                                         >  2:4275447/1‑93 (MQ=255)
                gATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGcc                                                                        >  2:803028/1‑93 (MQ=255)
                gATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGcc                                                                        >  1:4401533/1‑93 (MQ=255)
                  ttAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAg                                                                      <  1:1414145/93‑1 (MQ=255)
                    aaGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAggg                                                                    >  2:3882203/1‑92 (MQ=255)
                        aTACATTCCCACAACGAGGCGCTGTCACCGATAGAAACCTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTc                                                                <  2:349010/93‑1 (MQ=255)
                        aTACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTc                                                                >  2:2120841/1‑93 (MQ=255)
                         tACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCa                                                               >  2:225561/1‑93 (MQ=255)
                          aCATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCAt                                                              <  2:266232/93‑1 (MQ=255)
                            aTTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGa                                                            >  2:2793439/1‑93 (MQ=255)
                            aTTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGa                                                            <  1:4579081/93‑1 (MQ=255)
                             ttCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAg                                                           >  2:4352915/1‑93 (MQ=255)
                             ttCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAg                                                           <  2:3748497/93‑1 (MQ=255)
                                ccACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGt                                                        >  2:4093167/1‑93 (MQ=255)
                                ccACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGt                                                        >  1:3731757/1‑93 (MQ=255)
                                ccACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGt                                                        >  2:2757003/1‑93 (MQ=255)
                                ccACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGt                                                        >  2:2538731/1‑93 (MQ=255)
                                 cacaACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGtt                                                       >  1:1274509/1‑93 (MQ=255)
                                 cacaACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGtt                                                       >  2:727719/1‑93 (MQ=255)
                                 cacaACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGtt                                                       <  2:3225998/93‑1 (MQ=255)
                                   caACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTc                                                     <  1:1743365/93‑1 (MQ=255)
                                   caACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTc                                                     <  2:2438007/93‑1 (MQ=255)
                                   caACGAGGAGCTGGCACCGATAGAAAACTCGTTATGTCGTGTGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTc                                                     >  2:4547939/1‑93 (MQ=255)
                                    aaCGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCt                                                    <  2:1242365/93‑1 (MQ=255)
                                    aaCGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCt                                                    <  2:3550157/93‑1 (MQ=255)
                                     aCGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTg                                                   >  1:4463886/1‑93 (MQ=255)
                                     aCGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTg                                                   >  1:1058019/1‑93 (MQ=255)
                                      cGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGc                                                  <  2:2994378/93‑1 (MQ=255)
                                        aGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCAt                                                >  1:1597721/1‑93 (MQ=255)
                                          gcgcTGGCACCAATAGAAACCTCGTTATTTCGTGAGTTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCTGAGCGTTTCTGCATAAg                                             <  1:3556067/93‑1 (MQ=255)
                                            gcTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAgg                                            >  1:3339238/1‑93 (MQ=255)
                                             cTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGc                                           >  2:4487519/1‑93 (MQ=255)
                                              tGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCa                                          <  1:3288337/93‑1 (MQ=255)
                                               ggCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCag                                         <  1:2036777/93‑1 (MQ=255)
                                               ggCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCag                                         >  1:2368217/1‑93 (MQ=255)
                                               ggCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCag                                         <  1:273356/93‑1 (MQ=255)
                                                 cACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagta                                       <  2:3034295/93‑1 (MQ=255)
                                                 cACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagta                                       <  2:1303491/93‑1 (MQ=255)
                                                 cACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagta                                       <  1:659778/93‑1 (MQ=255)
                                                 cACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagta                                       <  1:3373271/93‑1 (MQ=255)
                                                  aCCGATAGAAAACTCGTTATGTCGTGAGGGATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagtag                                      >  1:3004939/1‑93 (MQ=255)
                                                   ccGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagtagt                                     <  2:1363430/93‑1 (MQ=255)
                                                   ccGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagtagt                                     >  1:2756562/1‑93 (MQ=255)
                                                   ccGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagtagt                                     >  2:3306824/1‑93 (MQ=255)
                                                    cGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGtt                                    <  1:544636/93‑1 (MQ=255)
                                                     gATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGttt                                   >  1:758189/1‑93 (MQ=255)
                                                      aTAGAAAACTGGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGtttt                                  >  2:2868924/1‑93 (MQ=255)
                                                      aTAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGtttt                                  <  2:1963904/93‑1 (MQ=255)
                                                      aTAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGtttt                                  >  2:1563052/1‑93 (MQ=255)
                                                        aGAAAACTCGTTATGTCGTGAGGTATGCGCCACCGCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTc                                <  1:1915289/93‑1 (MQ=255)
                                                        aGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTc                                >  2:4265310/1‑93 (MQ=255)
                                                        aGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTc                                <  1:4256990/93‑1 (MQ=255)
                                                        aGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTGCTGCATAAGGCAGTAGTTTTTc                                >  1:973111/1‑93 (MQ=255)
                                                         gAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCa                               <  1:3627080/93‑1 (MQ=255)
                                                          aaaaCTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCAc                              >  1:4051393/1‑93 (MQ=255)
                                                            aaCTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCg                            <  1:3571795/93‑1 (MQ=255)
                                                            aaCTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCg                            >  1:2503060/1‑93 (MQ=255)
                                                             aCTCGTTATGTCGTGAGGTATGCGCCGCCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCgg                           <  1:3811489/93‑1 (MQ=255)
                                                              cTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGc                          >  1:1687071/1‑93 (MQ=255)
                                                               tCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCa                         >  1:4668285/1‑93 (MQ=255)
                                                                  ttATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCg                      <  2:4660758/93‑1 (MQ=255)
                                                                   tATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCgg                     >  1:1927747/1‑93 (MQ=255)
                                                                   tATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCgg                     >  1:2106137/1‑93 (MQ=255)
                                                                    aTGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGt                    <  1:3678365/93‑1 (MQ=255)
                                                                    aTGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGt                    <  1:2311883/93‑1 (MQ=255)
                                                                     tGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGtt                   <  1:3451376/93‑1 (MQ=255)
                                                                        cGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTaaa                <  1:519434/93‑1 (MQ=255)
                                                                          tGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACg              >  2:777450/1‑93 (MQ=255)
                                                                           gAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACgg             >  1:2868537/1‑93 (MQ=255)
                                                                            aGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTCGTTTTTCACCGGCAGCGGTTAAACGGa            <  1:1868794/93‑1 (MQ=255)
                                                                            aGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACTGa            >  2:707057/1‑93 (MQ=255)
                                                                            aGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGa            <  1:873606/93‑1 (MQ=255)
                                                                             ggTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGat           <  1:3779513/93‑1 (MQ=255)
                                                                             ggTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGat           <  1:385988/93‑1 (MQ=255)
                                                                             ggTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGat           >  1:4297535/1‑93 (MQ=255)
                                                                             ggTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGat           >  2:1657064/1‑93 (MQ=255)
                                                                               tATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGatat         >  2:3460543/1‑93 (MQ=255)
                                                                                  gcgcCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGt      <  1:1148340/93‑1 (MQ=255)
                                                                                    gcCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGttt    <  1:2155519/93‑1 (MQ=255)
                                                                                     ccACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTc   >  2:4414073/1‑93 (MQ=255)
                                                                                     ccACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTc   <  2:17508/93‑1 (MQ=255)
                                                                                      cACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCt  >  1:4094139/1‑93 (MQ=255)
                                                                                            |                                                                                       
AAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCT  >  NC_000913/3945260‑3945437

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