Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 502,019 T→C Q148Q (CAA→CAG ybaL ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913502,0190TC98.0% 461.7 / ‑2.7 153Q148Q (CAA→CAGybaLputative transport protein
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (67/83);  total (69/84)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 163 total aligned reads displayed.

GCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACAC  >  NC_000913/501929‑502108
                                                                                          |                                                                                         
ccccTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTg                                                                                         <  2:868564/92‑1 (MQ=255)
 cccTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGATCCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGc                                                                                        <  2:854499/93‑1 (MQ=255)
 cccTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGc                                                                                        >  2:4268308/1‑93 (MQ=255)
  ccTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACGATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGcc                                                                                       >  2:4620182/1‑93 (MQ=255)
  ccTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGcc                                                                                       >  2:3109570/1‑93 (MQ=255)
   cTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGccc                                                                                      <  2:248400/93‑1 (MQ=255)
    tGTTCCATCATCCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCa                                                                                     <  2:4087258/93‑1 (MQ=255)
      ttCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACg                                                                                   >  2:2520112/1‑93 (MQ=255)
       tCCATCAGTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGc                                                                                  <  2:1448181/93‑1 (MQ=255)
         catcatTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTg                                                                                >  1:2009660/1‑93 (MQ=255)
         catcatTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTg                                                                                >  2:1594937/1‑93 (MQ=255)
          atcatTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGa                                                                               >  1:3363630/1‑93 (MQ=255)
           tcatTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGAc                                                                              >  1:3555623/1‑93 (MQ=255)
           tcatTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGAc                                                                              >  1:2273423/1‑93 (MQ=255)
              tcccTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGt                                                                           <  1:2001179/91‑1 (MQ=255)
               tCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTc                                                                          <  1:2875144/93‑1 (MQ=255)
                ccTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCa                                                                         <  2:4655783/93‑1 (MQ=255)
                 cTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCaa                                                                        >  1:1912507/1‑93 (MQ=255)
                 cTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCaa                                                                        >  1:2045240/1‑93 (MQ=255)
                 cTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCaa                                                                        <  1:1239340/93‑1 (MQ=255)
                 cTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCaa                                                                        >  2:4562372/1‑93 (MQ=255)
                   gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCaatt                                                                      >  2:811195/1‑93 (MQ=255)
                         gCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACACTGACTGTCAATTAATTGc                                                                >  2:3081276/1‑93 (MQ=255)
                         gCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGCCGATCTGCCCACGCTGACTGTCAATTAATTGc                                                                <  1:1854286/93‑1 (MQ=255)
                            gtcagcaACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGt                                                             <  1:1029391/91‑1 (MQ=255)
                            ggCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGt                                                             <  1:2347501/93‑1 (MQ=255)
                              gcacaACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGttc                                                           <  2:4331792/90‑1 (MQ=255)
                               agcaACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGttct                                                          >  2:2471/1‑93 (MQ=255)
                               agcaACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGttct                                                          >  1:105653/1‑93 (MQ=255)
                               agcaACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGttct                                                          >  2:208378/1‑93 (MQ=255)
                               agcaACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGttct                                                          <  2:48246/93‑1 (MQ=255)
                                 caACACCAGCGTCAGAACCATTCCCAGGTCTGCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGttcttc                                                        <  2:4221779/93‑1 (MQ=255)
                                 caACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGttcttc                                                        <  1:516917/93‑1 (MQ=255)
                                  aaCACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCa                                                       >  2:3139495/1‑93 (MQ=255)
                                  aaCACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCa                                                       <  1:3147934/93‑1 (MQ=255)
                                  aaCACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCa                                                       <  1:1801013/93‑1 (MQ=255)
                                  aaCACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCa                                                       <  1:170587/93‑1 (MQ=255)
                                  aaCACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCa                                                       <  1:1499159/93‑1 (MQ=255)
                                   acacCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCaa                                                      <  2:4441583/93‑1 (MQ=255)
                                   acacCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCaa                                                      >  2:1773310/1‑93 (MQ=255)
                                   acacCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCaa                                                      <  2:1578733/93‑1 (MQ=255)
                                   acacCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCaa                                                      <  2:90392/93‑1 (MQ=255)
                                   acacCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCaa                                                      <  2:1154701/93‑1 (MQ=255)
                                     acCAGCGTCAGAGCCATTACCAGGTCTGCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGt                                                    <  2:499044/93‑1 (MQ=255)
                                     acCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGt                                                    >  1:720063/1‑93 (MQ=255)
                                     acCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGt                                                    >  1:596095/1‑93 (MQ=255)
                                     acCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGt                                                    <  2:3356538/93‑1 (MQ=255)
                                        aGCGTCAGAACCATTCCCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTgcg                                                 <  1:3568611/93‑1 (MQ=255)
                                         gCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTgcgc                                                >  1:3253042/1‑93 (MQ=255)
                                            tCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCa                                             >  1:3472656/1‑93 (MQ=255)
                                                aaCCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTaa                                         <  2:3990489/93‑1 (MQ=255)
                                                aaCCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTaa                                         <  2:3203085/93‑1 (MQ=255)
                                                      tACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAAcaccac                                   >  2:514669/1‑93 (MQ=255)
                                                      tACCAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAAcaccac                                   <  2:422103/93‑1 (MQ=255)
                                                        ccAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACgg                                 >  1:3877577/1‑93 (MQ=255)
                                                        ccAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACgg                                 <  1:760442/93‑1 (MQ=255)
                                                         cAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGt                                >  2:4285251/1‑93 (MQ=255)
                                                         cAGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGt                                <  1:717742/93‑1 (MQ=255)
                                                          aGGTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTa                               <  1:3824415/93‑1 (MQ=255)
                                                           ggTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTAc                              <  1:921664/93‑1 (MQ=255)
                                                           ggTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTAc                              >  2:2641474/1‑93 (MQ=255)
                                                            gTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACt                             <  2:1045645/93‑1 (MQ=255)
                                                            gTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACt                             >  1:1385348/1‑93 (MQ=255)
                                                            gTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACt                             >  2:1992908/1‑93 (MQ=255)
                                                            gTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACt                             >  2:2285995/1‑93 (MQ=255)
                                                            gTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACt                             <  1:2928317/93‑1 (MQ=255)
                                                            gTCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACt                             <  1:4008689/93‑1 (MQ=255)
                                                             tCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTg                            >  2:646611/1‑93 (MQ=255)
                                                             tCTTCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTg                            >  2:2495890/1‑93 (MQ=255)
                                                               ttCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGc                          <  1:1098683/93‑1 (MQ=255)
                                                               ttCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGc                          <  2:3682773/93‑1 (MQ=255)
                                                               ttCCACAATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGc                          >  1:1691455/1‑93 (MQ=255)
                                                                tCCACAATCAACCAACCGATGGCGATCTCCCCACGCTGACTGTCAATTAATTGCCGTTCTTCCAGCGCGCGCAGTAACACCACGGTACTGGCg                         >  2:612838/1‑93 (MQ=255)
                                                                  cacaATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGt                       <  1:2377394/93‑1 (MQ=255)
                                                                  cacaATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGt                       <  1:436836/93‑1 (MQ=255)
                                                                  cacaATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTGCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGt                       <  2:1045055/93‑1 (MQ=255)
                                                                   acaATCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTg                      <  1:4515934/93‑1 (MQ=255)
                                                                     aaTCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGa                    >  2:2814438/1‑93 (MQ=255)
                                                                       tCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGaaa                  <  1:2222781/93‑1 (MQ=255)
                                                                       tCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGaaa                  <  1:461514/93‑1 (MQ=255)
                                                                       tCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGaaa                  >  1:4579014/1‑93 (MQ=255)
                                                                       tCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGaaa                  <  1:4070785/93‑1 (MQ=255)
                                                                       tCAACCAACCGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGaaa                  <  1:4178515/93‑1 (MQ=255)
                                                                         aaccaaccGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGa                >  2:1896918/1‑93 (MQ=255)
                                                                         aaccaaccGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGa                <  2:215312/93‑1 (MQ=255)
                                                                          accaaccGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGAc               <  1:1117261/93‑1 (MQ=255)
                                                                           ccaaccGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACa              <  1:791278/93‑1 (MQ=255)
                                                                            caaccGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACAt             <  1:3771038/93‑1 (MQ=255)
                                                                            caaccGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTAGCGGTGGAAAGACAt             >  1:1400985/1‑93 (MQ=255)
                                                                             aaccGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATa            <  2:1188565/93‑1 (MQ=255)
                                                                             aaccGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATa            <  1:2479284/93‑1 (MQ=255)
                                                                             aaccGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGAAACACCACGGTACTGGCGGTGGAAAGACATa            <  2:4641334/93‑1 (MQ=255)
                                                                              accGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATaa           >  1:503764/1‑93 (MQ=255)
                                                                              accGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATaa           <  2:3971513/93‑1 (MQ=255)
                                                                                cGATGGCGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAAc         <  2:137462/93‑1 (MQ=255)
                                                                                      cGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAaca   <  2:3686723/93‑1 (MQ=255)
                                                                                      cGATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAaca   >  2:113269/1‑93 (MQ=255)
                                                                                       gATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAacac  <  1:4403996/93‑1 (MQ=255)
                                                                                       gATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAacac  <  2:3769763/93‑1 (MQ=255)
                                                                                       gATCTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAacac  >  2:1155176/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACAC  >  NC_000913/501929‑502108

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