Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 232,218 A→G D98G (GAC→GGC)  yafE → putative biotin synthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913232,2180AG100.0% 795.5 / NA 248D98G (GAC→GGC) yafEputative biotin synthesis protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (126/122);  total (126/122)

Only 100 of 253 total aligned reads displayed.

GTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTACGCC  >  NC_000913/232128‑232306
                                                                                          |                                                                                        
gTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCg                                                                                        <  1:5126322/93‑1 (MQ=255)
gTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCg                                                                                        >  2:6035355/1‑93 (MQ=255)
 ttATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGt                                                                                       <  2:3883848/93‑1 (MQ=255)
    ttCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAAt                                                                                    >  2:6491960/1‑93 (MQ=255)
    ttCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAAt                                                                                    <  1:619401/93‑1 (MQ=255)
     tCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATg                                                                                   >  2:2511031/1‑93 (MQ=255)
     tCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATg                                                                                   <  1:3161871/93‑1 (MQ=255)
     tCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATg                                                                                   <  2:3976912/93‑1 (MQ=255)
      cTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGt                                                                                  >  1:513568/1‑93 (MQ=255)
       tGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGtc                                                                                 >  1:2489329/1‑93 (MQ=255)
       tGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGtc                                                                                 >  1:5677071/1‑93 (MQ=255)
        gCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGtct                                                                                >  1:5770711/1‑93 (MQ=255)
        gCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGtct                                                                                >  1:2359754/1‑93 (MQ=255)
        gCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGtct                                                                                >  1:1636683/1‑93 (MQ=255)
         cccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGtctc                                                                               <  2:1977634/93‑1 (MQ=255)
          ccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTcc                                                                              >  1:328702/1‑93 (MQ=255)
           cATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCg                                                                             >  1:5198577/1‑93 (MQ=255)
            aTCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCgg                                                                            <  2:4449813/93‑1 (MQ=255)
            aTCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCgg                                                                            <  2:4279523/93‑1 (MQ=255)
             tCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCggg                                                                           <  1:4139699/93‑1 (MQ=255)
             tCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCggg                                                                           >  2:1346825/1‑93 (MQ=255)
              cACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGt                                                                          >  2:4667677/1‑93 (MQ=255)
               aCTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGTGCGTAATGTCTCCGGGTc                                                                         >  1:5950795/1‑93 (MQ=255)
                cTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCa                                                                        >  2:2937923/1‑93 (MQ=255)
                  ggCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCAcc                                                                      >  1:4217272/1‑93 (MQ=255)
                    cATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCa                                                                    >  2:3126996/1‑93 (MQ=255)
                     atgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAg                                                                   <  2:4303370/93‑1 (MQ=255)
                     atgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAg                                                                   >  2:4089364/1‑93 (MQ=255)
                        atgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTgc                                                                <  2:5986556/93‑1 (MQ=255)
                         tgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTgcg                                                               <  1:3020073/93‑1 (MQ=255)
                          gTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGGgcgc                                                              >  2:2207263/1‑93 (MQ=255)
                                gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                                        >  2:3063687/1‑93 (MQ=255)
                                gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                                        <  2:2567330/93‑1 (MQ=255)
                                gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                                        >  2:4209754/1‑93 (MQ=255)
                                gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACACAGTGCGCGACATc                                                        <  2:5024358/93‑1 (MQ=255)
                                   gcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTgg                                                     >  2:5761266/1‑93 (MQ=255)
                                   gcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTgg                                                     >  2:1723696/1‑93 (MQ=255)
                                    caCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGt                                                    >  2:3645384/1‑93 (MQ=255)
                                    caCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGt                                                    <  1:1053349/93‑1 (MQ=255)
                                     aCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGtt                                                   <  1:688007/93‑1 (MQ=255)
                                      cTGCGAGAAGTGAATAGGATATTGAACCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTa                                                  <  2:3278421/93‑1 (MQ=255)
                                       tGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTAc                                                 <  1:2607731/93‑1 (MQ=255)
                                       tGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTAc                                                 >  2:3543940/1‑93 (MQ=255)
                                        gCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACa                                                <  2:6414120/93‑1 (MQ=255)
                                            gaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                                            >  1:6309951/1‑93 (MQ=255)
                                            gaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                                            >  1:4051563/1‑93 (MQ=255)
                                             aaGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACgg                                           >  1:4693047/1‑93 (MQ=255)
                                             aaGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACgg                                           <  1:5031705/93‑1 (MQ=255)
                                              aGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGt                                          >  2:5431228/1‑93 (MQ=255)
                                              aGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGt                                          <  1:3505406/93‑1 (MQ=255)
                                              aGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGt                                          >  2:5632804/1‑93 (MQ=255)
                                              aGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGt                                          <  2:4766024/93‑1 (MQ=255)
                                               gTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                                         >  1:4884609/1‑93 (MQ=255)
                                                 gAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGa                                       >  2:5859665/1‑93 (MQ=255)
                                                 gAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGa                                       >  2:5299538/1‑93 (MQ=255)
                                                  aaTAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGaa                                      >  1:30504/1‑93 (MQ=255)
                                                  aaTAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGaa                                      >  2:1695489/1‑93 (MQ=255)
                                                  aaTAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGaa                                      >  1:2193521/1‑93 (MQ=255)
                                                  aaTAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGaa                                      <  2:6202144/93‑1 (MQ=255)
                                                   aTAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAg                                     >  2:1993276/1‑93 (MQ=255)
                                                     aGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGGCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCa                                   >  1:5851813/1‑93 (MQ=255)
                                                      ggATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAt                                  <  1:1927433/93‑1 (MQ=255)
                                                       gATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAtt                                 <  1:1525326/93‑1 (MQ=255)
                                                         tatTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAc                               <  2:3680949/93‑1 (MQ=255)
                                                         tatTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAc                               <  2:4856148/93‑1 (MQ=255)
                                                         tatTGAAACCTGGCGCTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAc                               <  2:2158315/93‑1 (MQ=255)
                                                         aatTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAc                               <  2:2811389/92‑1 (MQ=255)
                                                          atTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcg                              >  1:1430071/1‑93 (MQ=255)
                                                           ttGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcgc                             <  1:3835255/93‑1 (MQ=255)
                                                           ttGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcgc                             <  2:1235290/93‑1 (MQ=255)
                                                             gAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa                           <  2:2856008/93‑1 (MQ=255)
                                                             gAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa                           >  1:1160163/1‑93 (MQ=255)
                                                              aaaCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGAt                          <  1:4381231/93‑1 (MQ=255)
                                                               aaCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATa                         >  1:4530432/1‑93 (MQ=255)
                                                               aaCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATa                         >  1:2857020/1‑93 (MQ=255)
                                                               aaCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATa                         >  2:3765061/1‑93 (MQ=255)
                                                               aaCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATa                         <  2:1293107/93‑1 (MQ=255)
                                                                aCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAc                        >  1:5563982/1‑93 (MQ=255)
                                                                aCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAc                        >  1:6445900/1‑93 (MQ=255)
                                                                aCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAc                        >  2:4238501/1‑93 (MQ=255)
                                                                 ccTGTCGGGAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACGTCTGGTTACAGACGGTAGAAGCATTACGCGATAcc                       <  1:2812317/93‑1 (MQ=255)
                                                                 ccTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAcc                       >  2:1522024/1‑93 (MQ=255)
                                                                 ccTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAcc                       <  2:116162/93‑1 (MQ=255)
                                                                  cTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACct                      >  2:76786/1‑93 (MQ=255)
                                                                   tGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctc                     <  1:2064352/93‑1 (MQ=255)
                                                                    ggCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctct                    <  2:5555735/93‑1 (MQ=255)
                                                                      cGGTAGGCTGATTGTGATGGGCGTAATGGCTCCGGGTCACCCAGAGCGCGCCATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCa                  >  1:3557650/1‑93 (MQ=255)
                                                                        gTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACg                >  2:2689575/1‑93 (MQ=255)
                                                                           ggCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTAc             >  1:6541637/1‑93 (MQ=255)
                                                                             cTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGa           >  2:3488612/1‑93 (MQ=255)
                                                                               gATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGaaa         >  1:1234375/1‑93 (MQ=255)
                                                                                gttGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCCCGTACGAAAc        <  1:4807723/92‑1 (MQ=255)
                                                                                aTTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAAc        >  2:2607193/1‑93 (MQ=255)
                                                                                 ttGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACt       >  2:2316773/1‑93 (MQ=255)
                                                                                  tgtgATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTa      <  2:4535613/93‑1 (MQ=255)
                                                                                   gtgATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTAc     >  1:931853/1‑93 (MQ=255)
                                                                                   gtgATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTAc     <  2:2367027/93‑1 (MQ=255)
                                                                                   gtgATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTAc     <  2:3424419/93‑1 (MQ=255)
                                                                                     gATGGGCGTAATGTCTCCGGCTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACTAAACTACGc   >  1:521611/1‑93 (MQ=255)
                                                                                     gATGGGCGTAATGTCTCCGGCTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTACGc   <  1:3905433/93‑1 (MQ=255)
                                                                                          |                                                                                        
GTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTACGCC  >  NC_000913/232128‑232306

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