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.C97.5% 660.9 / 7.6 200coding (2/363 nt)araHABC‑type arabinose transport system permease component
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base C (116/79);  total (117/83)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 220 total aligned reads displayed.

CCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACGGT  >  NC_000913/1980847‑1981026
                                                                                              |                                                                                       
ccAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCgg                                                                                          <  1:2751125/93‑1 (MQ=255)
 cAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGc                                                                                        <  1:3908920/93‑1 (MQ=255)
   gAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcct                                                                                      >  1:3845756/1‑91 (MQ=255)
     ataatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGcctct                                                                                    <  2:1062356/93‑5 (MQ=255)
       aatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTc                                                                                   >  2:55761/1‑93 (MQ=255)
         ttCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGt                                                                                 <  2:2666040/93‑1 (MQ=255)
         ttCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGt                                                                                 >  2:457881/1‑93 (MQ=255)
            cGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTcc                                                                              >  1:1960990/1‑93 (MQ=255)
              gCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCg                                                                            >  1:5201597/1‑93 (MQ=255)
                tATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGcc                                                                          <  2:479475/93‑1 (MQ=255)
                 aTCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCa                                                                         >  1:4884561/1‑93 (MQ=255)
                  tCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCaa                                                                        >  2:4418333/1‑93 (MQ=255)
                  tCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCaa                                                                        <  1:2430039/93‑1 (MQ=255)
                  tCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCaa                                                                        >  2:5079176/1‑93 (MQ=255)
                   cGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAAt                                                                       <  2:5745686/93‑1 (MQ=255)
                    gCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATc                                                                      <  2:2566636/93‑1 (MQ=255)
                     cTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCg                                                                     >  1:1393441/1‑93 (MQ=255)
                       gATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGcc                                                                   >  1:108449/1‑93 (MQ=255)
                         tACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                                                 >  1:4385656/1‑93 (MQ=255)
                         tACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                                                 <  1:2597653/93‑1 (MQ=255)
                         tACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                                                 >  1:6278190/1‑93 (MQ=255)
                         tACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                                                 <  2:1468603/93‑1 (MQ=255)
                          aCCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAgg                                                                <  1:1535876/93‑1 (MQ=255)
                          aCCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAgg                                                                >  1:5258603/1‑93 (MQ=255)
                           ccAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCATCGTTCCCGCCAATCGCCAggg                                                               >  1:4617591/1‑93 (MQ=255)
                            cAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgt                                                              >  2:2693590/1‑93 (MQ=255)
                             aGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtg                                                             >  1:5678333/1‑93 (MQ=255)
                             aGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtg                                                             >  1:4044084/1‑93 (MQ=255)
                             aGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtg                                                             >  1:2652599/1‑93 (MQ=255)
                             aGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtg                                                             <  1:5172479/93‑1 (MQ=255)
                              ggCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtgt                                                            <  2:4980025/93‑1 (MQ=255)
                              ggCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtgt                                                            <  1:6376803/93‑1 (MQ=255)
                                  tGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACg                                                        >  1:3578101/1‑93 (MQ=255)
                                  tGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACg                                                        >  2:163535/1‑93 (MQ=255)
                                  tGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACg                                                        >  2:4142417/1‑93 (MQ=255)
                                    agagaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGAc                                                      <  1:2327392/93‑1 (MQ=255)
                                     gagaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGAcc                                                     >  1:803427/1‑93 (MQ=255)
                                      agaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCTCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCa                                                    >  2:3436398/1‑93 (MQ=255)
                                      agaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCa                                                    >  2:1678351/1‑93 (MQ=255)
                                       gaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCaa                                                   >  2:5105484/1‑93 (MQ=255)
                                       gaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCaa                                                   >  2:62808/1‑93 (MQ=255)
                                       gaACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCaa                                                   >  1:3331365/1‑93 (MQ=255)
                                        aaCAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCaaa                                                  <  1:5530814/93‑1 (MQ=255)
                                        aaCAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCaaa                                                  <  1:6009585/93‑1 (MQ=255)
                                         aCAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAg                                                 >  1:5790577/1‑93 (MQ=255)
                                          cAAAGATAATATTTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAgg                                                <  2:1362988/93‑1 (MQ=255)
                                           aaaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggt                                               >  1:5895155/1‑93 (MQ=255)
                                           aaaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggg                                               >  1:5604390/1‑92 (MQ=255)
                                           aaaGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggg                                               >  1:290463/1‑92 (MQ=255)
                                            aaGATAATATTTTTGGTGCGAACACCCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtg                                              <  1:2291492/93‑1 (MQ=255)
                                              gATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtggt                                            <  1:3435023/93‑1 (MQ=255)
                                              gATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtggt                                            >  1:2065524/1‑93 (MQ=255)
                                              gATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtggt                                            >  2:3164774/1‑93 (MQ=255)
                                               ataataATTTTGGTGCTAACACCCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGtt                                           <  1:184015/93‑1 (MQ=255)
                                               ataataATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGtt                                           <  1:538422/93‑1 (MQ=255)
                                                  ataatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTa                                        >  2:4625342/1‑93 (MQ=255)
                                                  ataatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTa                                        <  1:365796/93‑1 (MQ=255)
                                                    aatTTTGGTGCTACAACCCGGTCCACCCGCCAGCCGCGCGGCCTCTTCGTCCCCGCCAATCGCCAGCGTGTTACGGCCAAAGGTGGTTTTAtt                                      <  2:3546745/93‑1 (MQ=255)
                                                    aatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTAtt                                      >  2:3778803/1‑93 (MQ=255)
                                                    aatTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTAtt                                      <  2:6464608/93‑1 (MQ=255)
                                                     atTTTGGTGCGAACAACCGTTACACCCGCCAGACGCGCGTCCTCTTCGTCCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTc                                     <  2:6021742/93‑1 (MQ=255)
                                                      ttttGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTca                                    >  2:717090/1‑93 (MQ=255)
                                                       tttGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcag                                   >  1:2378816/1‑93 (MQ=255)
                                                          ggTGCGAACACCCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagcag                                <  1:2958495/93‑1 (MQ=255)
                                                          ggTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagcag                                >  2:4046146/1‑93 (MQ=255)
                                                          ggTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagcag                                >  2:5545702/1‑93 (MQ=255)
                                                           gTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCagcagc                               >  1:2024807/1‑93 (MQ=255)
                                                            tGCGAACAACCGGTACCCCCGCCAGACGCGCGCCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAgcagca                              <  1:1675713/93‑1 (MQ=255)
                                                              cGAACAACCGGTACCCCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCaaa                            <  1:3026602/93‑1 (MQ=255)
                                                               gAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAAc                           >  1:5016181/1‑93 (MQ=255)
                                                                aacaacCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAAcc                          >  2:175230/1‑93 (MQ=255)
                                                                 acaacCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCa                         >  1:2152387/1‑93 (MQ=255)
                                                                   agccGGTACCCCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCaaa                       <  1:737233/91‑1 (MQ=255)
                                                                   aacCGGTACCCCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCaaa                       <  2:5164957/93‑1 (MQ=255)
                                                                   aacCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCaaa                       >  1:1281227/1‑93 (MQ=255)
                                                                    acCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAg                      >  2:5442183/1‑93 (MQ=255)
                                                                    acCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAg                      >  1:4746156/1‑93 (MQ=255)
                                                                    acCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAg                      >  1:6247062/1‑93 (MQ=255)
                                                                     cgggTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGa                     <  2:1727599/91‑1 (MQ=255)
                                                                      cGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGAt                    <  2:3705647/93‑1 (MQ=255)
                                                                       ggTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATa                   >  1:631205/1‑93 (MQ=255)
                                                                       ggTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGCGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATa                   >  2:201348/1‑93 (MQ=255)
                                                                        gTACCCCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATaa                  <  2:1197636/93‑1 (MQ=255)
                                                                        gTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATaa                  >  1:3543748/1‑93 (MQ=255)
                                                                         tACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAAt                 <  2:2029028/93‑1 (MQ=255)
                                                                         tACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAAt                 <  2:3019879/93‑1 (MQ=255)
                                                                          acccccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATc                <  1:3732822/90‑1 (MQ=255)
                                                                           cacCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCa               <  1:4318267/93‑1 (MQ=255)
                                                                             cccGCCAGACGCGCGGCCTCTTCGTTCCCGCCATTCGCCGGGGTGTTACGACCAAAGGTGGTTTTATTCAGCGGCAAACCAAAGATAATCAGa             >  1:4332732/1‑93 (MQ=255)
                                                                             cccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAATCAAAGATAATCAGa             >  1:4523122/1‑93 (MQ=255)
                                                                             cccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGa             <  2:4509613/93‑1 (MQ=255)
                                                                             cccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGa             >  1:2998162/1‑93 (MQ=255)
                                                                              ccGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGac            <  1:2568721/93‑1 (MQ=255)
                                                                                 ccAGACGCGCGCCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACAcg         <  2:5510759/93‑1 (MQ=255)
                                                                                   aGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACAcgcg       >  2:332215/1‑93 (MQ=255)
                                                                                   aGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACAcgcg       >  1:6129645/1‑93 (MQ=255)
                                                                                    gACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGa      >  1:2806587/1‑93 (MQ=255)
                                                                                     aCGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGAc     <  2:2872579/93‑1 (MQ=255)
                                                                                        cgcgGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACGGt  <  2:5362314/93‑1 (MQ=255)
                                                                                        cgcgGCCTCTTCGTACCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACGGt  <  2:5859911/93‑1 (MQ=255)
                                                                                              |                                                                                       
CCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAACCAAAGATAATCAGACACGCGACGGT  >  NC_000913/1980847‑1981026

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