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.1% 753.4 / ‑1.1 213coding (2/363 nt)araHABC‑type arabinose transport system permease component
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (129/82);  total (129/84)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 230 total aligned reads displayed.

CAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACGGT  >  NC_000913/1980848‑1981026
                                                                                            |                                                                                       
cAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGTCGCGCGGc                                                                                        <  1:6489037/93‑1 (MQ=255)
cAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGc                                                                                        >  1:3120790/1‑93 (MQ=255)
cAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGc                                                                                        <  2:5787423/93‑1 (MQ=255)
cAGAATAATTCCGGCTATCGCTGATACCAGGCCTGACAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTCCACCCGCCAGACGCGCGGc                                                                                        <  2:4220095/93‑1 (MQ=255)
 aGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcc                                                                                       >  1:6387325/1‑92 (MQ=255)
  gAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcct                                                                                      >  1:5785545/1‑91 (MQ=255)
  gAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcct                                                                                      >  1:767018/1‑91 (MQ=255)
    ataatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcctct                                                                                    <  1:20391/93‑5 (MQ=255)
      aatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTc                                                                                   >  2:1906505/1‑93 (MQ=255)
      aatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTc                                                                                   >  1:368443/1‑93 (MQ=255)
         tCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGtt                                                                                <  1:5584318/93‑1 (MQ=255)
           cGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACCCCCGCCAGACGCGCGGCCTCTTCGTTcc                                                                              <  1:1501172/93‑1 (MQ=255)
           cGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTcc                                                                              <  1:2773925/93‑1 (MQ=255)
            ggCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTccc                                                                             >  2:1312912/1‑93 (MQ=255)
               tATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGcc                                                                          >  2:1761534/1‑93 (MQ=255)
                 tCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTCCACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCaa                                                                        <  1:5077801/93‑1 (MQ=255)
                 tCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCaa                                                                        <  1:5941489/93‑1 (MQ=255)
                  cGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAAt                                                                       <  1:2821886/93‑1 (MQ=255)
                   gCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATc                                                                      <  1:2349387/93‑1 (MQ=255)
                   gCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATc                                                                      >  1:3852266/1‑93 (MQ=255)
                   gCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGCTCCCGCCAATc                                                                      >  2:1091117/1‑93 (MQ=255)
                    cTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCg                                                                     >  1:6310040/1‑93 (MQ=255)
                     tGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCGCGCCAATCGc                                                                    >  1:6434427/1‑93 (MQ=255)
                     tGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGc                                                                    >  2:2324863/1‑93 (MQ=255)
                      gATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACCCCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGcc                                                                   <  2:3734721/93‑1 (MQ=255)
                      gATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGcc                                                                   >  2:6308040/1‑93 (MQ=255)
                      gATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGcc                                                                   <  1:3937346/93‑1 (MQ=255)
                      gATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGcc                                                                   >  2:356330/1‑93 (MQ=255)
                        tACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                                                 <  1:6595979/93‑1 (MQ=255)
                        tACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                                                 <  2:3454276/93‑1 (MQ=255)
                        tACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                                                 >  2:3341481/1‑93 (MQ=255)
                          ccAGGCCTGAGAGAACAAAGATAATAATTTGGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAggg                                                               <  1:1894567/93‑1 (MQ=255)
                           cAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgt                                                              >  2:3915871/1‑93 (MQ=255)
                             ggCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtgt                                                            >  2:2700298/1‑93 (MQ=255)
                             ggCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtgt                                                            >  2:4300500/1‑93 (MQ=255)
                             ggCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtgt                                                            <  2:5104957/93‑1 (MQ=255)
                              gCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGtt                                                           <  2:1139443/93‑1 (MQ=255)
                                cTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTAc                                                         >  2:5545106/1‑93 (MQ=255)
                                cTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTAc                                                         >  2:2102361/1‑93 (MQ=255)
                                cTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTAc                                                         >  2:6556407/1‑93 (MQ=255)
                                cTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTAc                                                         <  1:1299832/93‑1 (MQ=255)
                                 tGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACg                                                        >  1:474178/1‑93 (MQ=255)
                                  gagagaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGa                                                       >  2:2429686/1‑93 (MQ=255)
                                  gagagaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGa                                                       >  2:2780341/1‑93 (MQ=255)
                                  gagagaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCACGGCGTTACGa                                                       >  1:3243404/1‑93 (MQ=255)
                                   agagaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGAc                                                      <  2:6468428/93‑1 (MQ=255)
                                     agaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCa                                                    >  2:1363581/1‑93 (MQ=255)
                                      gaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCaa                                                   >  1:1606646/1‑93 (MQ=255)
                                       aaCAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCaag                                                  >  1:1732001/1‑92 (MQ=255)
                                        aCAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAg                                                 >  2:3510555/1‑93 (MQ=255)
                                         cAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAgg                                                <  2:1563283/93‑1 (MQ=255)
                                          aaaGATAATAATTTTGGTGCGAACAACCGGTACCCCCGCCAGACGCGCGGCCTCTTCGTCCCCGCCAATCGCCAGGGTGTTACGACCAAAggt                                               <  1:2318779/93‑1 (MQ=255)
                                          aaaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggt                                               <  1:3590828/93‑1 (MQ=255)
                                           aaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              <  2:2567407/93‑1 (MQ=255)
                                           aaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              >  1:6258137/1‑93 (MQ=255)
                                           aaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              >  2:1962717/1‑93 (MQ=255)
                                            aGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtgg                                             >  1:701925/1‑93 (MQ=255)
                                            aGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGAGCGGCCTCTTCGTTCCCGCCAATCGCCAGCGCGTTACGACCAAaggagg                                             >  2:4398183/1‑93 (MQ=255)
                                             gATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtggt                                            >  1:4373255/1‑93 (MQ=255)
                                             gATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtggt                                            >  1:5026032/1‑93 (MQ=255)
                                             aataataATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtggt                                            <  2:2445909/92‑1 (MQ=255)
                                                aataatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGtttt                                         >  1:2497598/1‑93 (MQ=255)
                                                aataatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGtttt                                         >  1:2693233/1‑93 (MQ=255)
                                                aataatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGtttt                                         <  1:3512914/93‑1 (MQ=255)
                                                  taatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTAt                                       <  1:4238186/93‑1 (MQ=255)
                                                  taatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTAt                                       >  1:5262413/1‑93 (MQ=255)
                                                  taatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTAt                                       >  2:308920/1‑93 (MQ=255)
                                                    atTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAGAGGTGGTTTTATTc                                     <  1:4052124/93‑1 (MQ=255)
                                                    atTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTc                                     <  2:4935382/93‑1 (MQ=255)
                                                     ttttGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTca                                    >  2:4568792/1‑93 (MQ=255)
                                                     ttttGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTca                                    <  1:4564245/93‑1 (MQ=255)
                                                        tGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCACTCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagca                                 >  1:1103494/1‑93 (MQ=255)
                                                        tGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagca                                 >  2:2998897/1‑93 (MQ=255)
                                                         ggTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagcag                                >  1:3173110/1‑93 (MQ=255)
                                                         ggTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagcag                                >  1:983196/1‑93 (MQ=255)
                                                          gTGCGAACACCCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCagcagc                               <  2:6505647/93‑1 (MQ=255)
                                                          gTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCagcagc                               <  1:2217113/93‑1 (MQ=255)
                                                           tGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAgcagca                              >  2:5383747/1‑93 (MQ=255)
                                                            gCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCaa                             >  2:2167207/1‑93 (MQ=255)
                                                             cGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCaaa                            >  1:4799398/1‑93 (MQ=255)
                                                              gAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAAc                           >  1:961074/1‑93 (MQ=255)
                                                               aacaacCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAAcc                          <  2:6202224/93‑1 (MQ=255)
                                                                     cGGTACACCAGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGCGGTTTTATTCAGCAGCAAACCAAAGAt                    >  2:4825642/1‑93 (MQ=255)
                                                                      ggTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATa                   >  2:2670581/1‑93 (MQ=255)
                                                                       gTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATaa                  <  2:2254721/93‑1 (MQ=255)
                                                                        tACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAAt                 <  1:6041424/93‑1 (MQ=255)
                                                                          cccccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCa               <  2:3413237/91‑1 (MQ=255)
                                                                            cccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGa             <  2:5475188/93‑1 (MQ=255)
                                                                            cccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGa             >  2:6140028/1‑93 (MQ=255)
                                                                             ccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGac            >  2:1231625/1‑93 (MQ=255)
                                                                               gCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGacac          <  2:4718993/93‑1 (MQ=255)
                                                                               gCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGacac          <  1:4470699/93‑1 (MQ=255)
                                                                               cccAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTGTTCAGCAGCAAACCAAAGATAATCAGacac          <  1:3999509/92‑1 (MQ=255)
                                                                                 cAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACAcgc        <  2:4118370/93‑1 (MQ=255)
                                                                                  aGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACAcgcg       >  1:4639730/1‑93 (MQ=255)
                                                                                  aGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACAcgcg       >  2:2970434/1‑93 (MQ=255)
                                                                                   gACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGa      >  2:3627302/1‑93 (MQ=255)
                                                                                     cgcgcgGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACg    >  2:2836371/1‑93 (MQ=255)
                                                                                      gcgcgGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACgg   >  2:5251837/1‑93 (MQ=255)
                                                                                       cgcgGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACGGt  <  1:4299833/93‑1 (MQ=255)
                                                                                            |                                                                                       
CAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACGGT  >  NC_000913/1980848‑1981026

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