Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,980,940 +C coding (1/363 nt)
coding (627/627 nt)
araH ←
araH ←
ABC‑type arabinose transport system permease component
ABC‑type arabinose transport system permease component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,980,9391.C99.4% 564.4 / ‑3.9 161coding (2/363 nt)araHABC‑type arabinose transport system permease component
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (96/64);  total (96/65)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 185 total aligned reads displayed.

CCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACGGT  >  NC_000913/1980847‑1981026
                                                                                             |                                                                                       
ccAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCgg                                                                                          >  1:3210061/1‑93 (MQ=255)
 cAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGc                                                                                        <  1:1017739/93‑1 (MQ=255)
  aGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcc                                                                                       >  2:3608230/1‑92 (MQ=255)
  aGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcc                                                                                       >  1:2365761/1‑92 (MQ=255)
  aGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcc                                                                                       <  2:556648/93‑2 (MQ=255)
  aGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcc                                                                                       >  2:752919/1‑92 (MQ=255)
   gAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcct                                                                                      >  2:334358/1‑91 (MQ=255)
   gAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAAACGGTACACCCGCCAGACGCGCGGcct                                                                                      >  2:1278693/1‑91 (MQ=255)
    aataatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcctt                                                                                     <  2:3610101/93‑4 (MQ=255)
    aataatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcctc                                                                                     <  1:4601873/93‑4 (MQ=255)
    aataatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcctc                                                                                     >  2:3731539/1‑90 (MQ=255)
     ataatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcctct                                                                                    >  2:3706733/1‑89 (MQ=255)
     ataatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcctct                                                                                    >  2:3682582/1‑89 (MQ=255)
     ataatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcctct                                                                                    >  2:5961170/1‑89 (MQ=255)
         ttCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGt                                                                                 >  2:5600883/1‑93 (MQ=255)
          tCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGtt                                                                                <  1:6070598/93‑1 (MQ=255)
          tCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGtt                                                                                <  2:5786803/93‑1 (MQ=255)
            cGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTcc                                                                              <  2:5452719/93‑1 (MQ=255)
            cGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTcc                                                                              <  2:203908/93‑1 (MQ=255)
            cGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTcc                                                                              >  1:3936745/1‑93 (MQ=255)
               cTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGc                                                                           >  2:5855412/1‑93 (MQ=255)
               cTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGc                                                                           >  2:2885867/1‑93 (MQ=255)
               cTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGc                                                                           >  2:2013611/1‑93 (MQ=255)
                tATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGcc                                                                          >  2:1430820/1‑93 (MQ=255)
                 aTCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCa                                                                         >  2:1595784/1‑93 (MQ=255)
                   cGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAAt                                                                       <  2:2625912/93‑1 (MQ=255)
                    gCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATc                                                                      >  1:2762360/1‑93 (MQ=255)
                     cTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCg                                                                     >  1:1628845/1‑93 (MQ=255)
                        aTACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCa                                                                  <  2:3147485/93‑1 (MQ=255)
                           ccAGTCTTGAGAAAAAAAAGTTAATAATGTGGGTGCGAACACCCGGGACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAggg                                                               <  2:1251963/93‑1 (MQ=255)
                           ccAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACCCCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAggg                                                               <  2:2252861/93‑1 (MQ=255)
                           ccAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAggg                                                               <  1:1023901/93‑1 (MQ=255)
                             aGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtg                                                             <  1:3719745/93‑1 (MQ=255)
                             aGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtg                                                             <  2:2371925/93‑1 (MQ=255)
                             aGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtg                                                             >  1:4208261/1‑93 (MQ=255)
                              ggCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtgt                                                            >  2:659317/1‑93 (MQ=255)
                              ggCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtgt                                                            >  2:5566517/1‑93 (MQ=255)
                               gCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGtt                                                           >  1:5522206/1‑93 (MQ=255)
                                  tGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACg                                                        >  2:5171640/1‑93 (MQ=255)
                                  tGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACg                                                        >  2:4043072/1‑93 (MQ=255)
                                  tGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACg                                                        >  1:3119550/1‑93 (MQ=255)
                                  tGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACg                                                        >  1:950454/1‑93 (MQ=255)
                                   gagagaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGa                                                       >  1:2534229/1‑93 (MQ=255)
                                     gagaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGAcc                                                     >  1:4351192/1‑93 (MQ=255)
                                         aCAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAg                                                 >  2:2631918/1‑93 (MQ=255)
                                         aCAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAg                                                 >  1:3897584/1‑93 (MQ=255)
                                         aCAAAGAAAATAATTTTGGTGCGAACAACCGGTACCCCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAg                                                 <  1:4692805/93‑1 (MQ=255)
                                          cAAAGATAATAATTTTGGTGCGAACAACCGTTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAgg                                                >  2:2940874/1‑93 (MQ=255)
                                          aaaaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAgg                                                <  1:1424684/92‑1 (MQ=255)
                                          aaaaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAgg                                                <  2:1965641/92‑1 (MQ=255)
                                           aaaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACTCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggg                                               >  2:4784813/1‑92 (MQ=255)
                                           aaaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggt                                               >  1:5642221/1‑93 (MQ=255)
                                           aaaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggt                                               >  1:5791770/1‑93 (MQ=255)
                                            aaGATAATAATTTTGGTGCGAACACCCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              <  1:4086722/93‑1 (MQ=255)
                                            aaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              >  1:5179736/1‑93 (MQ=255)
                                            aaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              <  2:4857339/93‑1 (MQ=255)
                                            aaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              >  2:5248979/1‑93 (MQ=255)
                                             aGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggttg                                             >  2:4319520/1‑91 (MQ=255)
                                             aGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtgg                                             >  1:4878159/1‑93 (MQ=255)
                                              gATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCCATCGCCAGGGTGTTACGACCAAAggtggt                                            >  1:416301/1‑93 (MQ=255)
                                              gATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtggt                                            >  1:4733509/1‑93 (MQ=255)
                                                taataaTTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGttt                                          >  1:1821231/1‑93 (MQ=255)
                                                   taatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTAt                                       >  1:5922571/1‑93 (MQ=255)
                                                    aatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTAtt                                      >  2:1259403/1‑93 (MQ=255)
                                                     atTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTc                                     <  1:3515912/93‑1 (MQ=255)
                                                       tttGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcag                                   >  1:4676028/1‑93 (MQ=255)
                                                       tttGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcag                                   <  1:4854310/93‑1 (MQ=255)
                                                       tttGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTaag                                   >  1:4500847/1‑93 (MQ=255)
                                                        ttGGTGCGAACACCCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagc                                  <  1:5736348/93‑1 (MQ=255)
                                                        ttGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagc                                  >  1:5599514/1‑93 (MQ=255)
                                                         tGGTCCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagca                                 <  2:3040905/93‑1 (MQ=255)
                                                          ggTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagcag                                >  1:2591291/1‑93 (MQ=255)
                                                           gTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCagcagc                               >  2:2771062/1‑93 (MQ=255)
                                                            tGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAgcagca                              >  2:4966421/1‑93 (MQ=255)
                                                               gAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAAc                           >  2:20501/1‑93 (MQ=255)
                                                                aacaacCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAAcc                          >  1:797838/1‑93 (MQ=255)
                                                                 acaacCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCa                         >  1:607562/1‑93 (MQ=255)
                                                                  caacCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCaa                        >  2:3513552/1‑93 (MQ=255)
                                                                   aacCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCaaa                       <  1:1377573/93‑1 (MQ=255)
                                                                     ccGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGa                     >  2:1752208/1‑93 (MQ=255)
                                                                     ccGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGa                     <  1:2578629/93‑1 (MQ=255)
                                                                      cGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGAAAACCAAAGAt                    >  1:2881262/1‑93 (MQ=255)
                                                                       ggTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATa                   <  1:1160348/93‑1 (MQ=255)
                                                                          acacCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATc                <  1:3656369/93‑1 (MQ=255)
                                                                          acacCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATc                <  1:2126842/93‑1 (MQ=255)
                                                                          acacCCGCCAGACGCGCGGCCTCTTCGTACCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATc                <  2:3004631/93‑1 (MQ=255)
                                                                           cacCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCa               <  1:3745040/93‑1 (MQ=255)
                                                                            ccccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAg              <  2:1407035/92‑1 (MQ=255)
                                                                             cccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGa             <  1:2302979/93‑1 (MQ=255)
                                                                              ccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGac            <  1:663167/93‑1 (MQ=255)
                                                                              ccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGac            >  2:2989982/1‑93 (MQ=255)
                                                                               cGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGaca           <  2:4903087/93‑1 (MQ=255)
                                                                                 ccAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACAcg         >  1:5316370/1‑93 (MQ=255)
                                                                                  cAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACAcgc        <  2:1186280/93‑1 (MQ=255)
                                                                                   aGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGCGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACAcgcg       >  1:3320530/1‑93 (MQ=255)
                                                                                    gACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGa      <  1:5114953/93‑1 (MQ=255)
                                                                                     aCGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGAc     <  2:2269233/93‑1 (MQ=255)
                                                                                     aCGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGAc     <  2:4118096/93‑1 (MQ=255)
                                                                                     aCGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGAc     <  1:4241505/93‑1 (MQ=255)
                                                                                        cgcgGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACGGt  <  1:2953193/93‑1 (MQ=255)
                                                                                             |                                                                                       
CCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACGGT  >  NC_000913/1980847‑1981026

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