Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,907,388 T→C E213E (GAA→GAG pstC ← PstC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,907,3880TC100.0% 774.9 / NA 243E213E (GAA→GAGpstCPstC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (116/127);  total (116/127)

Only 100 of 254 total aligned reads displayed.

CATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTA  >  NC_000913/3907298‑3907477
                                                                                          |                                                                                         
  tgttttCCCCGAGCGCGCGGCCCAGCCCCAGCATGAGGCCCCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTccc                                                                                       <  2:2701042/90‑1 (MQ=255)
  tGGTTTCACCGAGCGCGCGGCCCCGCCCCCGCATGATGCCGCCGATAACACCATTTTTTGTGGACGGAAGAAAGATACGCCAGATAAACTccc                                                                                       >  2:4031475/1‑93 (MQ=255)
  tGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTccc                                                                                       <  1:1863028/93‑1 (MQ=255)
   ggTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCa                                                                                      <  2:1797491/93‑1 (MQ=255)
    gTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAg                                                                                     <  1:2028827/93‑1 (MQ=255)
    gTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAg                                                                                     >  2:2927258/1‑93 (MQ=255)
    gTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAg                                                                                     <  1:1846268/93‑1 (MQ=255)
      ttCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAggt                                                                                   <  2:3213752/93‑1 (MQ=255)
       tCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAggtg                                                                                  <  1:1907454/93‑1 (MQ=255)
        cACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAggtgg                                                                                 <  1:93408/93‑1 (MQ=255)
         cccGAGCGCGCGGCCAAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAggtggt                                                                                <  2:4171001/92‑1 (MQ=255)
         aCCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAggtggg                                                                                >  2:3802557/1‑92 (MQ=255)
          ccGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCGCAGgttgtg                                                                               >  2:2285133/1‑93 (MQ=255)
          ccGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGgtggtg                                                                               <  1:147952/93‑1 (MQ=255)
            gAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCa                                                                             >  1:745386/1‑93 (MQ=255)
            gAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCa                                                                             <  1:3556156/93‑1 (MQ=255)
             aGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAg                                                                            <  1:3470207/93‑1 (MQ=255)
             aGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAg                                                                            >  2:2853106/1‑93 (MQ=255)
                 cgcgGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCaa                                                                        <  1:1107966/93‑1 (MQ=255)
                   cgGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAAtc                                                                      >  1:2302161/1‑92 (MQ=255)
                   cgGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATa                                                                      <  1:4537240/93‑1 (MQ=255)
                    ggCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATAc                                                                     <  1:4366683/93‑1 (MQ=255)
                    ggCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATAc                                                                     >  1:4492955/1‑93 (MQ=255)
                        cacccccagaATGAAGCCTCCGATAACACCATTGTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTa                                                                 <  1:383074/90‑1 (MQ=255)
                         agccccagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAg                                                                <  2:2786220/93‑1 (MQ=255)
                              cagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGa                                                           >  2:3016487/1‑93 (MQ=255)
                               agcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGAc                                                          <  1:1225535/93‑1 (MQ=255)
                               agcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGAc                                                          <  2:752363/93‑1 (MQ=255)
                                 cATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGActc                                                        >  1:2803169/1‑93 (MQ=255)
                                 cATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGActc                                                        <  1:2135718/93‑1 (MQ=255)
                                  atgatgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGActct                                                       >  1:582877/1‑93 (MQ=255)
                                  atgatgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGActct                                                       <  2:1100184/93‑1 (MQ=255)
                                   tgatgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCtt                                                      >  2:3148259/1‑93 (MQ=255)
                                   tgatgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCtt                                                      >  2:1607219/1‑93 (MQ=255)
                                    gatgCCGCTGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCttt                                                     <  2:1167033/93‑1 (MQ=255)
                                    gatgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCttt                                                     >  1:2257791/1‑93 (MQ=255)
                                     atgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtc                                                    <  1:1350965/93‑1 (MQ=255)
                                        ccgccgATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtcatc                                                 <  2:2729056/93‑1 (MQ=255)
                                        ccgccgATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtcatc                                                 <  1:2393288/93‑1 (MQ=255)
                                         cgccgATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtcatca                                                <  1:316785/93‑1 (MQ=255)
                                          gccgATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtcatcat                                               <  2:3062078/93‑1 (MQ=255)
                                          gccgATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtcatcat                                               <  2:4636614/93‑1 (MQ=255)
                                            cgATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTtcatcatca                                             >  1:2465703/1‑93 (MQ=255)
                                             gATAACACCATTTGTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATCCCGTAGGCCGACTCTTTCATCATCAc                                            <  1:1730085/93‑1 (MQ=255)
                                              aTAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCAcc                                           <  1:595638/93‑1 (MQ=255)
                                              aTAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCAcc                                           >  2:2865961/1‑93 (MQ=255)
                                               tAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCg                                          <  2:3934247/93‑1 (MQ=255)
                                               tAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCg                                          >  2:2101767/1‑93 (MQ=255)
                                               tAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCg                                          >  2:626237/1‑93 (MQ=255)
                                               tAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCg                                          >  2:660505/1‑93 (MQ=255)
                                                 acacCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCggg                                        <  2:1945003/93‑1 (MQ=255)
                                                 acacCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCggg                                        <  2:112630/93‑1 (MQ=255)
                                                 acacCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAAGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCggg                                        >  1:470680/1‑93 (MQ=255)
                                                  cacCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCgggg                                       >  2:1572998/1‑93 (MQ=255)
                                                   acCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGt                                      >  1:2963588/1‑93 (MQ=255)
                                                    ccATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGtt                                     >  1:763950/1‑93 (MQ=255)
                                                     cATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCCTCATCACCGGGGttt                                    >  2:4093190/1‑93 (MQ=255)
                                                     cATTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGttt                                    >  1:88907/1‑93 (MQ=255)
                                                      aTTTTTGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTg                                   >  2:4338248/1‑93 (MQ=255)
                                                        ttttGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGtt                                 <  1:2668698/93‑1 (MQ=255)
                                                        ttttGGTGAACGAAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGtt                                 <  2:2405010/93‑1 (MQ=255)
                                                          ttGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCg                               <  1:2625144/93‑1 (MQ=255)
                                                           tGGTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGCCTCTTTCATCATCACCGGGGTTTGTTCGa                              <  1:3289772/93‑1 (MQ=255)
                                                            ggTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGaa                             <  1:4393625/93‑1 (MQ=255)
                                                             gTGAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGGTTGTTCGAac                            >  2:3345486/1‑93 (MQ=255)
                                                             gTGAACGAAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAac                            <  2:80387/93‑1 (MQ=255)
                                                             gTGAACAGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAac                            <  2:2389810/93‑1 (MQ=255)
                                                              tGAACGGAACAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAaca                           <  1:1270342/93‑1 (MQ=255)
                                                               gAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAacac                          <  2:1362784/93‑1 (MQ=255)
                                                               gAACGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAacac                          <  2:3855768/93‑1 (MQ=255)
                                                                aaCGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAAcaca                         >  2:2465444/1‑93 (MQ=255)
                                                                 aCGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGATTGTTCGAACACAt                        >  1:4660475/1‑93 (MQ=255)
                                                                  cGGAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATc                       <  1:2548997/93‑1 (MQ=255)
                                                                   ggAAGAACGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCa                      >  2:3652379/1‑93 (MQ=255)
                                                                    gaagaaCGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCAc                     >  1:1312549/1‑93 (MQ=255)
                                                                    gaagaaCGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCAc                     >  2:1714463/1‑93 (MQ=255)
                                                                      agaaCGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGc                   <  2:1264868/93‑1 (MQ=255)
                                                                      agaaCGATACGCCAGATAACCTCCCAGGTGGGGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGc                   <  1:2557055/93‑1 (MQ=255)
                                                                      agaaCGAGACGCCAGATACCCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGc                   <  2:1754511/93‑1 (MQ=255)
                                                                       gaaCGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCa                  <  1:4422524/93‑1 (MQ=255)
                                                                       gaaCGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCa                  <  1:1912801/93‑1 (MQ=255)
                                                                          cGATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTa               >  2:1851532/1‑93 (MQ=255)
                                                                           gATACGCCAGATACCCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAc              <  1:1687856/93‑1 (MQ=255)
                                                                           gATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAc              <  2:2499848/93‑1 (MQ=255)
                                                                           gATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAc              >  2:4459422/1‑93 (MQ=255)
                                                                           gATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAc              >  2:3600945/1‑93 (MQ=255)
                                                                           gATACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAc              <  2:2384032/93‑1 (MQ=255)
                                                                           aaTACGCCAGATAGCCTCCCAGGTGGTTCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAc              <  2:4191041/92‑1 (MQ=255)
                                                                            aTACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAcc             >  2:2028944/1‑93 (MQ=255)
                                                                            aTACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAcc             >  1:1930563/1‑93 (MQ=255)
                                                                             tACGCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAccg            >  1:1017679/1‑93 (MQ=255)
                                                                              aCGACAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAccgc           <  1:1208567/93‑1 (MQ=255)
                                                                                gCCAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAccgccg         <  2:504738/93‑1 (MQ=255)
                                                                                 ccAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACcgccgc        >  1:3652197/1‑93 (MQ=255)
                                                                                  cAGATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCa       >  1:3196975/1‑93 (MQ=255)
                                                                                   aGATAACCTCCCAGGTGGGGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCaa      <  2:2575101/93‑1 (MQ=255)
                                                                                    gATAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAAt     <  1:3686390/93‑1 (MQ=255)
                                                                                     aTAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATg    >  2:3448740/1‑93 (MQ=255)
                                                                                      tAACCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGt   >  2:846231/1‑93 (MQ=255)
                                                                                       aaCCTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTa  <  1:4089091/93‑1 (MQ=255)
                                                                                          |                                                                                         
CATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTA  >  NC_000913/3907298‑3907477

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