Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 73,807 C→T 6.4% D239N (GAT→AAT)  thiP ← thiamine/thiamine pyrophosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091373,8070CT 6.4% 192.8 / 2.1 62D239N (GAT→AAT) thiPthiamine/thiamine pyrophosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base T (2/2);  total (32/30)
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.99e-01

CTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGG  >  NC_000913/73662‑73933
                                                                                                                                                 |                                                                                                                              
cTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTa                                                                                                                             >  5:19041/1‑149 (MQ=255)
  aTGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGt                                                                                                                           >  4:25999/1‑149 (MQ=255)
           aTCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCTATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACt                                                                                                                    >  6:79664/1‑147 (MQ=255)
           aTCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACt                                                                                                                    >  4:333273/1‑147 (MQ=255)
                    ggcgcgCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCACTCGCTGCCTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCCGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGAt                                                                                                         <  4:115323/147‑1 (MQ=255)
                     tCGCGACAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGCAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATa                                                                                                        <  3:241484/149‑1 (MQ=255)
                       gcgcCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAAACGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGa                                                                                                          >  1:55065/1‑145 (MQ=255)
                         gcCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGAtt                                                                                                    >  7:33317/1‑149 (MQ=255)
                         gcCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAAAGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGAtt                                                                                                    <  5:295941/149‑1 (MQ=255)
                           cAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGa                                                                                                      >  6:351228/1‑145 (MQ=255)
                           cAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGa                                                                                                      <  5:351228/145‑1 (MQ=255)
                                ttGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGATCAGCGCCAGCATTGCCCCGAGGGCAGTATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGca                                                                                             >  6:199524/1‑148 (MQ=255)
                                  gcagcaGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATc                                                                                                                                <  4:245044/112‑1 (MQ=255)
                                  gcagcaGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATc                                                                                                                                >  3:245044/1‑112 (MQ=255)
                                   cagcagCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGa                                                                                          <  3:333273/149‑1 (MQ=255)
                                    agcagcGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATTGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGa                                                                                          >  1:4418/1‑148 (MQ=255)
                                         cGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTgg                                                                                    <  2:25837/149‑1 (MQ=255)
                                              tgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGCGGTccg                                                                                >  4:332991/1‑146 (MQ=255)
                                                ccGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGc                                                                                                  >  3:173609/1‑128 (MQ=255)
                                                ccGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCt                                                                              <  8:33317/148‑1 (MQ=255)
                                                     cgcAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGAc                                                                        <  1:120390/149‑1 (MQ=255)
                                                      gcAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGAcc                                                                       >  6:187112/1‑149 (MQ=255)
                                                        aaTGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATc                                                                                                                                <  3:132203/90‑1 (MQ=255)
                                                        aaTGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATc                                                                                                                                >  4:132203/1‑90 (MQ=255)
                                                              cTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCa                                                               <  6:233651/149‑1 (MQ=255)
                                                              cTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCa                                                               >  6:253062/1‑149 (MQ=255)
                                                                     aaTCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTTGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAg                                                        <  5:79664/149‑1 (MQ=255)
                                                                     aaTCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATTGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAg                                                        <  7:153304/149‑1 (MQ=255)
                                                                     aaTCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATTGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAg                                                        >  8:153304/1‑149 (MQ=255)
                                                                     aaTCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATa                                                         <  5:187112/148‑1 (MQ=255)
                                                                     aaTCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAg                                                        >  5:276046/1‑149 (MQ=255)
                                                                       tCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCa                                                      <  8:79329/149‑1 (MQ=255)
                                                                         gCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACg                                                    <  2:55065/149‑1 (MQ=255)
                                                                          cTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACgg                                                   >  6:198663/1‑149 (MQ=255)
                                                                             aCTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTg                                                 <  5:198663/148‑1 (MQ=255)
                                                                             aCTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTg                                                 <  5:253062/148‑1 (MQ=255)
                                                                              cTCAACAGCACCAGCCCGAGGCATCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGc                                               <  5:199524/149‑1 (MQ=255)
                                                                               tCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATa                                                         <  4:267618/138‑1 (MQ=255)
                                                                               tCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATa                                                         >  3:267618/1‑138 (MQ=255)
                                                                                 aaCAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGc                                                                                              >  2:101085/1‑99 (MQ=255)
                                                                                 aaCAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCg                                              >  4:349714/1‑147 (MQ=255)
                                                                                 aaCAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGaa                                            <  3:349714/149‑1 (MQ=255)
                                                                                      cACCAGCCCGCGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATTGCACGGTGGCGAAGCTgg                                       >  3:76859/1‑149 (MQ=255)
                                                                                      cACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTgg                                       >  8:233429/1‑149 (MQ=255)
                                                                                         cAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCg                                     >  1:209844/1‑148 (MQ=255)
                                                                                               gAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACaga                              <  7:233429/149‑1 (MQ=255)
                                                                                               gAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCGGGCGAAACaga                              <  6:19041/149‑1 (MQ=255)
                                                                                                    agcaCACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATTGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATa                         <  2:4418/149‑1 (MQ=255)
                                                                                                    agcaCACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCAAGGTGGCGAAGCTGGCGAAACAGAGCAt                          >  4:14177/1‑148 (MQ=255)
                                                                                                         acCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGAt                    <  3:75372/149‑1 (MQ=255)
                                                                                                          ccATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTg                                        >  6:134784/1‑128 (MQ=255)
                                                                                                          ccATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTg                                        <  5:134784/128‑1 (MQ=255)
                                                                                                           cATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGAt                    >  2:149098/1‑147 (MQ=255)
                                                                                                           cATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGAt                    >  5:332873/1‑147 (MQ=255)
                                                                                                           cATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATaa                  <  2:209844/149‑1 (MQ=255)
                                                                                                               tGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGa                     <  4:7319/142‑1 (MQ=255)
                                                                                                               tGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGa                     >  3:7319/1‑142 (MQ=255)
                                                                                                                 gAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGc                                      >  6:267249/1‑123 (MQ=255)
                                                                                                                 gAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGc                                      <  5:267249/123‑1 (MQ=255)
                                                                                                                    cagcGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCaa         <  6:112673/149‑1 (MQ=255)
                                                                                                                     agcGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCACCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAAc        <  7:301849/149‑1 (MQ=255)
                                                                                                                        gcCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACcgg     >  8:141204/1‑149 (MQ=255)
                                                                                                                         ccAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACcggc    >  7:121102/1‑149 (MQ=255)
                                                                                                                         ccAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACcggc    >  3:248402/1‑149 (MQ=255)
                                                                                                                         ccAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACcggc    >  6:179384/1‑149 (MQ=255)
                                                                                                                           aGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACcggcgg  <  1:149098/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                              
CTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGG  >  NC_000913/73662‑73933

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: ‑