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.6% 946.7 / ‑3.9 281coding (2/363 nt)araHABC‑type arabinose transport system permease component
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (167/113);  total (167/114)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

Only 100 of 317 total aligned reads displayed.

CAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACGGT  >  NC_000913/1980848‑1981026
                                                                                            |                                                                                       
cAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGc                                                                                        >  1:8791529/1‑93 (MQ=255)
 aGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcc                                                                                       >  2:7111836/1‑92 (MQ=255)
 aGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcc                                                                                       <  1:1542981/93‑2 (MQ=255)
 aGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcc                                                                                       >  1:7357300/1‑92 (MQ=255)
  gAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcct                                                                                      <  1:4414763/93‑3 (MQ=255)
  gAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcct                                                                                      >  1:843395/1‑91 (MQ=255)
  gAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcct                                                                                      >  1:2579586/1‑91 (MQ=255)
   aataatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcctc                                                                                     <  1:368790/93‑4 (MQ=255)
    ataatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcctct                                                                                    >  2:271917/1‑89 (MQ=255)
       atTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGTTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCg                                                                                  <  1:7903247/93‑1 (MQ=255)
       atTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCg                                                                                  >  2:3370504/1‑93 (MQ=255)
        ttCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGt                                                                                 >  1:1885782/1‑93 (MQ=255)
         tCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGtt                                                                                <  1:8380973/93‑1 (MQ=255)
           cGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTcc                                                                              <  2:8634803/93‑1 (MQ=255)
              cTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGc                                                                           >  1:4651554/1‑93 (MQ=255)
              cTATCGCGGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGc                                                                           <  2:5558878/93‑1 (MQ=255)
               tATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGcc                                                                          >  2:6097998/1‑93 (MQ=255)
               tATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGcc                                                                          <  1:7106101/93‑1 (MQ=255)
               tATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGcc                                                                          >  1:558201/1‑93 (MQ=255)
                aTCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGATGCGCGGCCTCTTCGTTCCCGCCa                                                                         >  2:5914752/1‑93 (MQ=255)
                aTCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCa                                                                         >  2:2072903/1‑93 (MQ=255)
                 tCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCaa                                                                        >  1:111641/1‑93 (MQ=255)
                     tGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGc                                                                    >  1:1763626/1‑93 (MQ=255)
                      gATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGcc                                                                   >  1:8618659/1‑93 (MQ=255)
                      gATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGcc                                                                   <  2:4345/93‑1 (MQ=255)
                       aTACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCTAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCa                                                                  <  2:4223908/93‑1 (MQ=255)
                       aTACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCa                                                                  >  1:3117298/1‑93 (MQ=255)
                       aTACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCa                                                                  <  2:5255849/93‑1 (MQ=255)
                        tACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                                                 >  2:5668190/1‑93 (MQ=255)
                        tACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                                                 >  1:5695695/1‑93 (MQ=255)
                        aaCCAGGCCTGAGAGAACCAAGATAATAATTTTGGTGCGACCACCCGGTCCACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                                                 <  1:3433836/92‑1 (MQ=255)
                           cAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgt                                                              <  2:1232419/93‑1 (MQ=255)
                                cTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTAc                                                         <  2:3074506/93‑1 (MQ=255)
                                cTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTAc                                                         >  1:2311128/1‑93 (MQ=255)
                                 tGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACg                                                        >  2:2312440/1‑93 (MQ=255)
                                 tGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACg                                                        >  1:1196514/1‑93 (MQ=255)
                                 tGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACg                                                        >  1:6309342/1‑93 (MQ=255)
                                  gagagaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGa                                                       >  2:6494383/1‑93 (MQ=255)
                                   agagaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGAc                                                      <  1:4285282/93‑1 (MQ=255)
                                   agagaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGAc                                                      >  1:4780708/1‑93 (MQ=255)
                                   agagaACAAAAATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGAc                                                      <  1:3011626/93‑1 (MQ=255)
                                       aaaaaaGATAATAATTTTGGTGCGAACAACCGGTACCCCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCaaa                                                  <  1:6910409/90‑1 (MQ=255)
                                       aaCAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCaaa                                                  <  1:1010010/93‑1 (MQ=255)
                                        aCAAAGATAATAATTTTGGTGCGAGATGCCGGTACACCCGCCAGACGCGCGGACGCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAg                                                 <  2:2917855/93‑1 (MQ=255)
                                        aCAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAg                                                 >  2:177515/1‑93 (MQ=255)
                                        aCAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAg                                                 >  2:4145490/1‑93 (MQ=255)
                                         aaaaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAgg                                                <  2:6705486/92‑1 (MQ=255)
                                           aaGATAATAGTTTGGGGGCGAACAACCGCTGCCCCCGCCGGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              <  1:4550146/93‑1 (MQ=255)
                                           aaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              <  2:7687189/93‑1 (MQ=255)
                                           aaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              <  1:7246364/93‑1 (MQ=255)
                                           aaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              <  1:183471/93‑1 (MQ=255)
                                           aaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              >  1:5758392/1‑93 (MQ=255)
                                           aaGAGAATAATTGTGGTGCGAACAACCGGTACCCCCGCCAGCCGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              <  1:1486410/93‑1 (MQ=255)
                                            aGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtgg                                             >  2:2672416/1‑93 (MQ=255)
                                            aGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtgg                                             >  1:4885446/1‑93 (MQ=255)
                                            aGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtgg                                             <  2:5792818/93‑1 (MQ=255)
                                            aGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggggg                                             >  1:1972736/1‑93 (MQ=255)
                                             gATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCCATCGCCAGGGTGTTACGACCAAAggtggt                                            >  1:658721/1‑93 (MQ=255)
                                             gATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtggt                                            >  1:5878930/1‑93 (MQ=255)
                                              ataataATTTTGGTGCGAACAACCGGTACCCCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGtt                                           <  1:6329541/93‑1 (MQ=255)
                                              ataataATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtgggt                                           >  1:7707058/1‑91 (MQ=255)
                                              ataataATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGtt                                           <  1:3247814/93‑1 (MQ=255)
                                              ataataATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGtt                                           <  2:8063121/93‑1 (MQ=255)
                                               taataaTTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGttt                                          >  2:466407/1‑93 (MQ=255)
                                               taataaTTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGttt                                          <  2:8308781/93‑1 (MQ=255)
                                               aaataatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGttt                                          <  1:7476478/92‑1 (MQ=255)
                                                 ataatTTTGGTGCGAACACCCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTa                                        <  2:4557708/93‑1 (MQ=255)
                                                  taatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTAt                                       <  1:2409714/93‑1 (MQ=255)
                                                   aatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTAtt                                      >  2:7773999/1‑93 (MQ=255)
                                                   aatTTGGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTAtt                                      <  1:2097565/93‑1 (MQ=255)
                                                    atTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTc                                     >  1:202046/1‑93 (MQ=255)
                                                       ttGGTGAGAACAGCCGGTACACCCGCCAGACGGGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagc                                  <  2:686998/93‑1 (MQ=255)
                                                        tGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagca                                 <  2:2872261/93‑1 (MQ=255)
                                                        tGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagca                                 >  1:6438778/1‑93 (MQ=255)
                                                             ccaacaacCGCTCCACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCaaa                            <  2:2142834/91‑1 (MQ=255)
                                                               aacaacCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGTGTGTTACGACCAAAGGTGGTTTATTCAGCAGCAAACCa                         >  1:5214175/1‑93 (MQ=255)
                                                               aacaacCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAAcc                          >  2:7463967/1‑93 (MQ=255)
                                                                   acCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAg                      >  1:4053231/1‑93 (MQ=255)
                                                                    ccGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGa                     <  1:6763211/93‑1 (MQ=255)
                                                                    ccGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGa                     <  1:7629526/93‑1 (MQ=255)
                                                                     cGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGAt                    >  2:7352854/1‑93 (MQ=255)
                                                                      ggTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATa                   >  2:823888/1‑93 (MQ=255)
                                                                      ggTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATa                   <  1:5110752/93‑1 (MQ=255)
                                                                      ggTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATa                   >  2:3025318/1‑93 (MQ=255)
                                                                       gTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATaa                  >  1:2724688/1‑93 (MQ=255)
                                                                        tACACCCGCCAGACGCGCGGCCTCTTCGTTGCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAAt                 >  1:8045444/1‑93 (MQ=255)
                                                                         acacCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATc                <  1:5519094/93‑1 (MQ=255)
                                                                           acCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAg              >  2:6237006/1‑93 (MQ=255)
                                                                            ccccccAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGa             <  2:3848946/93‑1 (MQ=255)
                                                                            cccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGa             >  2:7871767/1‑93 (MQ=255)
                                                                             ccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGac            <  1:7834261/93‑1 (MQ=255)
                                                                              cGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGaca           >  2:1624597/1‑93 (MQ=255)
                                                                               gCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGacac          >  1:6041345/1‑93 (MQ=255)
                                                                               gCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGacac          <  1:2900024/93‑1 (MQ=255)
                                                                               gCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGacac          >  1:161096/1‑93 (MQ=255)
                                                                                ccAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACAcg         >  1:4168362/1‑93 (MQ=255)
                                                                                 cAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACAcgc        <  2:8715764/93‑1 (MQ=255)
                                                                                   gACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGa      >  2:5038458/1‑93 (MQ=255)
                                                                                       cgcgGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACggg  >  2:1935110/1‑92 (MQ=255)
                                                                                       cgcgGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACGGt  <  1:8271236/93‑1 (MQ=255)
                                                                                            |                                                                                       
CAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACGGT  >  NC_000913/1980848‑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: