Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 888,315 C→T G243D (GGC→GAC)  b0847 ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913888,3150CT99.6% 949.7 / NA 259G243D (GGC→GAC) b0847putative transport protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (132/126);  minor base G (0/1);  total (132/127)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 267 total aligned reads displayed.

ACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATCACC  >  NC_000913/888225‑888404
                                                                                          |                                                                                         
 ccAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCgg                                                                                        <  1:4667074/93‑1 (MQ=255)
    aCGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGttt                                                                                     >  2:181736/1‑93 (MQ=255)
         aTTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGa                                                                                >  1:2179086/1‑93 (MQ=255)
          tttCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGATTACGTTCAATGTAGCAGTCGGTTTGTCGAt                                                                               <  2:4687773/93‑1 (MQ=255)
          tttCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGAt                                                                               >  2:3819022/1‑93 (MQ=255)
          tttCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGAt                                                                               <  1:5486518/93‑1 (MQ=255)
          tttCATCGCACATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGGCGGTTTGTCGAt                                                                               <  2:1046393/93‑1 (MQ=255)
           ttCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATa                                                                              <  1:1102772/93‑1 (MQ=255)
           ttCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATa                                                                              >  2:5370011/1‑93 (MQ=255)
            tCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATaa                                                                             >  1:2686272/1‑93 (MQ=255)
            tCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATaa                                                                             >  1:5857867/1‑93 (MQ=255)
              aTCGCCCGTTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTAGGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAAt                                                                           <  1:303853/93‑1 (MQ=255)
              aTCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAAt                                                                           >  1:5133939/1‑93 (MQ=255)
              aTCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAAt                                                                           <  1:4982770/93‑1 (MQ=255)
               tCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATa                                                                          <  2:3478990/93‑1 (MQ=255)
               tCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATa                                                                          >  2:2642630/1‑93 (MQ=255)
                cGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATAc                                                                         <  1:3998662/93‑1 (MQ=255)
                 gCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATAcc                                                                        >  1:5425428/1‑93 (MQ=255)
                 gCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATAcc                                                                        >  1:1774106/1‑93 (MQ=255)
                 gCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGGTTGTCGATAAATAcc                                                                        >  2:5933853/1‑93 (MQ=255)
                  cccATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATAccc                                                                       >  2:3024895/1‑93 (MQ=255)
                  cccATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATAccc                                                                       >  1:5544914/1‑93 (MQ=255)
                   ccATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCa                                                                      >  1:4427711/1‑93 (MQ=255)
                    cATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAg                                                                     <  2:2365659/93‑1 (MQ=255)
                    cATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAg                                                                     <  2:4067354/93‑1 (MQ=255)
                    cATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAg                                                                     >  2:1273418/1‑93 (MQ=255)
                    cATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAg                                                                     <  1:277830/93‑1 (MQ=255)
                    cATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAg                                                                     <  1:2494102/93‑1 (MQ=255)
                     aTTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGt                                                                    <  1:3101977/93‑1 (MQ=255)
                     aTTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGt                                                                    >  2:3704/1‑93 (MQ=255)
                     aTTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAACACCCAGt                                                                    <  2:6196337/93‑1 (MQ=255)
                       ttGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTc                                                                  <  2:232438/93‑1 (MQ=255)
                        tGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCa                                                                 <  1:5274532/93‑1 (MQ=255)
                         gTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCAc                                                                <  2:5490997/93‑1 (MQ=255)
                         gTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCAc                                                                <  2:337487/93‑1 (MQ=255)
                           aGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGc                                                              >  2:5762069/1‑93 (MQ=255)
                           aGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGc                                                              <  1:3973939/93‑1 (MQ=255)
                            gCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCa                                                             >  1:3319984/1‑93 (MQ=255)
                            gCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCa                                                             <  2:2578681/93‑1 (MQ=255)
                                ggCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTTGATAAATACCCAGTTCACGCAGAtt                                                         <  2:2924301/93‑1 (MQ=255)
                                 gCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAttt                                                        >  2:4266424/1‑93 (MQ=255)
                                   aTCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAttttt                                                      >  1:613344/1‑93 (MQ=255)
                                    tCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTg                                                     >  1:1520877/1‑93 (MQ=255)
                                    tCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTg                                                     >  1:2065035/1‑93 (MQ=255)
                                     cACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGc                                                    <  2:3062002/93‑1 (MQ=255)
                                      aCCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAATTCACGCAGATTTTTGcc                                                   >  1:127035/1‑93 (MQ=255)
                                       ccGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCg                                                  >  1:790431/1‑93 (MQ=255)
                                        cGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGt                                                 >  1:5666558/1‑93 (MQ=255)
                                         gTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTc                                                >  2:4322987/1‑93 (MQ=255)
                                           cTGGATTTGCCAGAATCCCATTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCgg                                              <  2:3868096/93‑1 (MQ=255)
                                           cTGGATGTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCgg                                              <  2:2239003/93‑1 (MQ=255)
                                            tGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGt                                             >  2:5885205/1‑93 (MQ=255)
                                             ggATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTc                                            <  2:4865211/93‑1 (MQ=255)
                                               aTTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCa                                          >  1:5731338/1‑93 (MQ=255)
                                                tttGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAg                                         <  2:5064034/93‑1 (MQ=255)
                                                tttGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAg                                         >  2:5810614/1‑93 (MQ=255)
                                                tttGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAg                                         <  2:1140469/93‑1 (MQ=255)
                                                tttGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAg                                         >  1:5093725/1‑93 (MQ=255)
                                                   gCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGcc                                      >  2:2295458/1‑93 (MQ=255)
                                                   gCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGcc                                      <  1:3791656/93‑1 (MQ=255)
                                                    ccAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGcca                                     >  1:1030055/1‑93 (MQ=255)
                                                     cAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGccac                                    >  1:985466/1‑93 (MQ=255)
                                                       gAATCTCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGccacca                                  <  1:380523/93‑1 (MQ=255)
                                                       gAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGccacca                                  <  2:3182260/93‑1 (MQ=255)
                                                       gAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGccacca                                  <  1:4598604/93‑1 (MQ=255)
                                                       gAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGccacca                                  <  1:5638776/93‑1 (MQ=255)
                                                       gAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGccacca                                  <  1:5795865/93‑1 (MQ=255)
                                                       gAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGccacca                                  <  1:944679/93‑1 (MQ=255)
                                                       gAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTCGCCGTCGGTCCAGGccacca                                  <  2:2679041/93‑1 (MQ=255)
                                                         aTCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGt                                <  2:1553531/93‑1 (MQ=255)
                                                         aTCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGt                                >  2:75729/1‑93 (MQ=255)
                                                           cccGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTc                              <  1:2949859/93‑1 (MQ=255)
                                                             cGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCg                            <  2:3286954/93‑1 (MQ=255)
                                                              gTTACGTCGAATACGTTCAATGTGGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCgg                           >  1:5995136/1‑93 (MQ=255)
                                                              gTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCgg                           >  1:3389161/1‑93 (MQ=255)
                                                              gTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCgg                           <  2:5441963/93‑1 (MQ=255)
                                                               ttACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCggg                          >  1:4292107/1‑93 (MQ=255)
                                                                tACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGc                         <  2:916930/93‑1 (MQ=255)
                                                                tACGTAGGATACCTTCAATGGAGCAGTAGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGTTCCAGGCCACCAGTTCCGGGc                         <  1:1966860/93‑1 (MQ=255)
                                                                 aCGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGcc                        >  2:3070339/1‑93 (MQ=255)
                                                                   gTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGa                      >  2:5212554/1‑93 (MQ=255)
                                                                   gTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGa                      >  2:709457/1‑93 (MQ=255)
                                                                    tCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGAc                     >  1:1542838/1‑93 (MQ=255)
                                                                    tCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGAc                     <  1:3637212/93‑1 (MQ=255)
                                                                     cGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGAcg                    >  2:2991074/1‑93 (MQ=255)
                                                                        aTACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGa                 >  2:5144032/1‑93 (MQ=255)
                                                                         tACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGAt                >  2:4468047/1‑93 (MQ=255)
                                                                         tACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGAt                <  2:2815798/93‑1 (MQ=255)
                                                                         tACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGAt                <  1:5231690/93‑1 (MQ=255)
                                                                          aCGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATa               <  1:2352610/93‑1 (MQ=255)
                                                                          aCGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATa               <  1:6257709/93‑1 (MQ=255)
                                                                          aCGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATa               <  1:1682632/93‑1 (MQ=255)
                                                                          aCGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATa               <  2:562783/93‑1 (MQ=255)
                                                                            gTTCAATGTAGCAGTCGGTTCGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAgg             <  1:4205070/93‑1 (MQ=255)
                                                                               cAATGTAGCAGTCGTTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGgcgc          <  2:1956083/93‑1 (MQ=255)
                                                                                aaTGTAGCAGTCGGTTTGTTGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGcgcg         >  2:6121033/1‑93 (MQ=255)
                                                                                  tGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGa       >  1:4782385/1‑93 (MQ=255)
                                                                                  tGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGa       <  2:6261953/93‑1 (MQ=255)
                                                                                   gTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGAt      <  1:3229573/93‑1 (MQ=255)
                                                                                    tAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATc     >  2:3540913/1‑93 (MQ=255)
                                                                                          |                                                                                         
ACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATCACC  >  NC_000913/888225‑888404

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