Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,460,981 C→A 5.7% A30D (GCT→GAT)  paaJ → 3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,460,9810CA 5.7% 239.7 / 2.5 70A30D (GCT→GAT) paaJ3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base C (29/37);  new base A (2/2);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

CTGTTTCGCGGTAAATCTCACCGCATCGGCGGCACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCT  >  NC_000913/1460840‑1461122
                                                                                                                                             |                                                                                                                                             
cTGTTTCGCGGTAAATCTCACCGCATCGGCGGCACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATc                                                                                                                                        <  1:184137/149‑1 (MQ=255)
 tGTTTCGCGGTAAATCTCACCGCATCGGCGGCACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATcc                                                                                                                                       >  4:291694/1‑149 (MQ=255)
     tCGCGGTAAATCTCACCGCATCGGCGGCACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCttt                                                                                                                                   >  4:337718/1‑149 (MQ=255)
        cgGTAAATCTCACCGCATCGGCGGCACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGA‑GGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCgg                                                                                                                               <  4:240845/149‑1 (MQ=255)
             aaTCTCACCGCATCGGCGGCACCATTACAGGAGAAGCCTGATGCGTTAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGATAATATCCCTTTTCGGGAAc                                                                                                                           >  5:319130/1‑149 (MQ=255)
             aaTCTCACCGCATCGGCGGCACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAAc                                                                                                                           <  3:249695/149‑1 (MQ=255)
             aaTCTCACCGCATCGGCGGCAACATTACTGGAGAATCGTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGaa                                                                                                                            <  3:163374/148‑1 (MQ=255)
                    ccGCATCGGCGGCACCATTAGAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCAAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGt                                                                                                                    <  2:350165/149‑1 (MQ=255)
                    ccGCATCGGCGGCACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGt                                                                                                                    >  5:178979/1‑149 (MQ=255)
                     cGCATCGGCGGCACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATcc                                                                                                                                       >  8:22978/1‑129 (MQ=255)
                     cGCATCGGCGGCACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATcc                                                                                                                                       <  7:22978/129‑1 (MQ=255)
                           ggcggcACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGaaa                                                                                                              <  6:178979/148‑1 (MQ=255)
                           ggcggcACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAAc                                                                                                             >  8:191070/1‑149 (MQ=255)
                            gcggcCCCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGCTGCTCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAAcc                                                                                                            >  3:195285/1‑149 (MQ=255)
                            gcggcACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAAcc                                                                                                            >  6:231710/1‑149 (MQ=255)
                               gcACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCcgc                                                                                                         <  8:277957/149‑1 (MQ=255)
                                 aCCATTACAGGAGAAGCCTGATGCGTGAAGCCTGTATTTGTGACGGAATTCGTACGGCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGt                                                                                                       >  2:163027/1‑149 (MQ=255)
                                   cATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGATGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGtct                                                                                                     <  8:241464/149‑1 (MQ=255)
                                   cATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGATGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGtct                                                                                                     <  8:247124/149‑1 (MQ=255)
                                        cAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATg                                                                                                >  5:54500/1‑149 (MQ=255)
                                        cAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATg                                                                                                >  4:179925/1‑149 (MQ=255)
                                           gagaAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCg                                                                                              >  2:264921/1‑148 (MQ=255)
                                           gagaAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCg                                                                                              >  5:128022/1‑148 (MQ=255)
                                            agaAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGa                                                                                            <  1:310109/149‑1 (MQ=255)
                                            agaAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGa                                                                                            <  3:337718/149‑1 (MQ=255)
                                              aaGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAgt                                                                                          <  6:128022/149‑1 (MQ=255)
                                                gCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGTCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCATCGAGACGCgggtttc                                                                                        <  7:299423/149‑6 (MQ=255)
                                                  cTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTAt                                                                                      <  7:241018/149‑1 (MQ=255)
                                                     aTGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGGGCGAAACCCGCGTCTCGATGCGGAGTGTATCga                                                                                   <  3:291694/149‑1 (MQ=255)
                                                        cGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCgatga                                                                                <  5:80994/149‑1 (MQ=255)
                                                          tGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCATTTGCGTGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGatgatg                                                                              >  4:52283/1‑149 (MQ=255)
                                                           gAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGAtgt                                                                             >  1:181319/1‑149 (MQ=255)
                                                              gCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGAt                                                                          <  5:149094/149‑1 (MQ=255)
                                                                 tttATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCt                                                                       >  5:332777/1‑149 (MQ=255)
                                                                 tttATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCt                                                                       >  8:145059/1‑149 (MQ=255)
                                                                   tATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTAACGATGATGTGATCCt                                                                       <  7:119341/147‑1 (MQ=255)
                                                                       tgtgACGGAATTCGTACGCCAATTGGTCGCAACGGCGGGGCATTATCAAGTGGTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGGGTATCGATGATGTGATCCTCGGCt                                                                  <  1:323576/148‑1 (MQ=255)
                                                                         tgACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCtgt                                                                <  3:179925/148‑1 (MQ=255)
                                                                           aCGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGc                                                              >  8:317427/1‑148 (MQ=255)
                                                                           aCGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGc                                                              <  1:264921/148‑1 (MQ=255)
                                                                             ggAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCGGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTTATCATGGGCTGTGCTaa                                                           <  2:244069/149‑1 (MQ=255)
                                                                              gAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTa                                                                                       <  8:15682/120‑1 (MQ=255)
                                                                                    gTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGc                                                     <  2:181319/148‑1 (MQ=255)
                                                                                       cGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGga                                                 <  1:61733/149‑1 (MQ=255)
                                                                                       cGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGga                                                 <  1:122951/149‑1 (MQ=255)
                                                                                          cAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGgag                                                >  8:270054/1‑147 (MQ=255)
                                                                                           aaTTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCGGATGATCTGGCGGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCGGTGCTAATCAGGCGGGagaag                                             <  4:195285/149‑1 (MQ=255)
                                                                                                 tCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAAcc                                       <  6:201130/149‑1 (MQ=255)
                                                                                                 tCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAAcc                                       <  7:145059/149‑1 (MQ=255)
                                                                                                    cTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGt                                     >  6:172736/1‑148 (MQ=255)
                                                                                                           ggggCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGc                             <  1:361718/149‑1 (MQ=255)
                                                                                                            gggCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGcc                            <  5:172736/149‑1 (MQ=255)
                                                                                                            gggCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGcc                            <  7:191070/149‑1 (MQ=255)
                                                                                                              gCAGTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCTTCTGGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTCACGTAGCCCg                          <  4:38771/149‑1 (MQ=255)
                                                                                                                aTTATCAAGTGTTGGGGCTGATGAGCTGGCGTATATCCCTTTGCGGGAACTGCTGGTGCGAATCCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGATTTGCTAATCAGTCGGGAGAAGATAAACTTAACGTATCCCGGa                        <  6:319130/149‑1 (MQ=255)
                                                                                                                    tCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAAc                                  <  4:145048/135‑1 (MQ=255)
                                                                                                                    tCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAAc                                  >  3:145048/1‑135 (MQ=255)
                                                                                                                    tCAAGTGTTCGGGCTGATGATCTGGATGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCAGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTATGCTGTGCTAATAAGGAGTGAGAAGATAACCGTAACGTAGCCCGGATGGc                    >  5:344600/1‑149 (MQ=255)
                                                                                                                     cAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTTCTAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGAGGGAGAAGATAACCGTAACGTAGCCCTGATg                      >  6:156246/1‑146 (MQ=255)
                                                                                                                     cAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGc                    >  5:192851/1‑148 (MQ=255)
                                                                                                                       aGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCgatgat                                                                               <  2:263855/87‑1 (MQ=255)
                                                                                                                       aGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCgatgat                                                                               >  1:263855/1‑87 (MQ=255)
                                                                                                                           ttCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGGCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTATTCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTa             >  2:263455/1‑149 (MQ=255)
                                                                                                                           atCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCTGATGGCGACTTTa             >  2:109903/2‑149 (MQ=255)
                                                                                                                            tCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTAc            >  1:905/1‑149 (MQ=255)
                                                                                                                            tCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTAc            >  1:94176/1‑149 (MQ=255)
                                                                                                                             cGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTAc            <  5:231710/148‑1 (MQ=255)
                                                                                                                                   gatgatCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGagaaga                                            <  2:235325/110‑1 (MQ=255)
                                                                                                                                   gatgatCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGagaaga                                            >  1:235325/1‑110 (MQ=255)
                                                                                                                                      gatCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCt  >  8:214023/1‑149 (MQ=255)
                                                                                                                                       agctgcctgcTATCACTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCt  <  6:54500/146‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                             
CTGTTTCGCGGTAAATCTCACCGCATCGGCGGCACCATTACAGGAGAAGCCTGATGCGTGAAGCCTTTATTTGTGACGGAATTCGTACGCCAATTGGTCGCTACGGCGGGGCATTATCAAGTGTTCGGGCTGATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCT  >  NC_000913/1460840‑1461122

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑