Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 681,448 G→A 5.6% A413T (GCT→ACT)  djlC → J domain‑containing HscC co‑chaperone; Hsc56

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913681,4480GA 5.6% 229.8 / 2.8 71A413T (GCT→ACT) djlCJ domain‑containing HscC co‑chaperone; Hsc56
Reads supporting (aligned to +/- strand):  ref base G (30/37);  new base A (2/2);  total (32/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 = 8.98e-01

GCGTGGTGGAAATCGAGAGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGC  >  NC_000913/681301‑681582
                                                                                                                                                   |                                                                                                                                      
gCGTGGTGGAAATCGAGAGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGc                                                                                                                                       <  1:156380/149‑1 (MQ=255)
 cGTGGTGGAAATCGAGAGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCt                                                                                                                                      <  4:341489/149‑1 (MQ=255)
      tggAAATCGAGAGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATca                                                                                                                                 >  4:504764/1‑148 (MQ=255)
       ggAAATCGAGAGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCa                                                                                                                                >  2:425008/1‑149 (MQ=255)
       ggAAATCGAGAGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATACTc                                                                                                                                     >  4:239837/1‑144 (MQ=255)
       ggAAATCGAGAGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATACTc                                                                                                                                     <  3:239837/144‑1 (MQ=255)
       ggAAATCGAGAGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATACTc                                                                                                                                     >  2:484288/1‑144 (MQ=255)
       ggAAATCGAGAGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATACTc                                                                                                                                     <  1:484288/144‑1 (MQ=255)
           aTCGAGAGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTag                                                                                                                            >  2:442121/1‑149 (MQ=255)
             cGAGAGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTagag                                                                                                                          <  1:366660/149‑1 (MQ=255)
              gagaGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAgg                                                                                                                         >  2:555769/1‑149 (MQ=255)
              gagaGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAgg                                                                                                                         <  2:442125/149‑1 (MQ=255)
                         gCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCAt                                                                                                                               <  1:16519/132‑1 (MQ=255)
                         gCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCAt                                                                                                                               >  2:16519/1‑132 (MQ=255)
                             gaCAGAGCAAGATTCACCGTTCCTGGATTTGTATTGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGa                                                                                                          >  2:107892/1‑149 (MQ=255)
                                    cAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCtttt                                                                                                   >  4:330332/1‑149 (MQ=255)
                                    cAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCtttt                                                                                                   <  4:261423/149‑1 (MQ=255)
                                    cAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCtttt                                                                                                   >  3:170121/1‑149 (MQ=255)
                                     aaGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTg                                                                                                  >  4:429425/1‑149 (MQ=255)
                                     aaGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTg                                                                                                  >  1:596702/1‑149 (MQ=255)
                                      aGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCCTTGTCACCCGGATGCTTTTGa                                                                                                 <  2:61195/149‑1 (MQ=255)
                                           cACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGaaa                                                                                            <  1:173797/149‑1 (MQ=255)
                                           cACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGaaa                                                                                            <  3:107653/149‑1 (MQ=255)
                                             ccGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGt                                                                                          <  3:330332/149‑1 (MQ=255)
                                                  ccTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAAc                                                                                     <  1:425008/149‑1 (MQ=255)
                                                    tGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTc                                                                                   <  1:442121/149‑1 (MQ=255)
                                                    tGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTc                                                                                   >  4:313209/1‑149 (MQ=255)
                                                      gATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGa                                                                                 <  4:491031/149‑1 (MQ=255)
                                                        tttGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGc                                                                                                                                       >  1:308671/1‑93 (MQ=255)
                                                        tttGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGc                                                                                                                                       <  2:308671/93‑1 (MQ=255)
                                                            tATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCgga                                                                           >  1:484515/1‑149 (MQ=255)
                                                            tATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCgga                                                                           <  3:429425/149‑1 (MQ=255)
                                                                gggTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGAt                                                                       <  1:149269/149‑1 (MQ=255)
                                                                    aGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAg                                                                   <  1:555769/149‑1 (MQ=255)
                                                                      cGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGAt                                                                 <  4:376477/149‑1 (MQ=255)
                                                                        gAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGt                                                               >  3:362496/1‑149 (MQ=255)
                                                                        gAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGt                                                               >  1:331795/1‑149 (MQ=255)
                                                                              gtagtaCTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGa                                                         <  2:596702/149‑1 (MQ=255)
                                                                              gtagtaCTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGa                                                         <  4:345082/149‑1 (MQ=255)
                                                                              gtagtaCTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGa                                                         >  4:424424/1‑149 (MQ=255)
                                                                               tagtaCTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAg                                                        >  1:465830/1‑149 (MQ=255)
                                                                                agtaCTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGc                                                       <  2:484515/149‑1 (MQ=255)
                                                                                agtaCTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGc                                                       <  4:239310/149‑1 (MQ=255)
                                                                                 gtaCTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGcc                                                      >  1:36598/1‑149 (MQ=255)
                                                                                   aCTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCaa                                                    <  2:311678/149‑1 (MQ=255)
                                                                                    cTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAac                                                   <  4:170121/149‑1 (MQ=255)
                                                                                     ttGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGTCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCggag                                                                          <  1:85863/125‑1 (MQ=255)
                                                                                     ttGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGTCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCggag                                                                          >  2:85863/1‑125 (MQ=255)
                                                                                         aGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGa                                                                                                 >  1:38743/1‑98 (MQ=255)
                                                                                         aGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGa                                                                                                 <  2:38743/98‑1 (MQ=255)
                                                                                         aGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTg                                              >  2:306703/1‑149 (MQ=255)
                                                                                              ttCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATGGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAAt                                         <  4:86583/149‑1 (MQ=255)
                                                                                                  gAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCttt                                                                                                    >  4:103294/1‑86 (MQ=255)
                                                                                                  gAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCttt                                                                                                    <  3:103294/86‑1 (MQ=255)
                                                                                                  gAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATg                                     >  3:378059/1‑149 (MQ=255)
                                                                                                  gAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATg                                     >  4:412845/1‑149 (MQ=255)
                                                                                                   aaaCAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATga                                    <  3:313209/149‑1 (MQ=255)
                                                                                                    aaCAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCtttt                                                                                                   <  1:534963/85‑1 (MQ=255)
                                                                                                    aaCAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCtttt                                                                                                   >  2:534963/1‑85 (MQ=255)
                                                                                                     aCAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATgaag                                  >  1:318972/1‑149 (MQ=255)
                                                                                                          agaTGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCAt                                      <  4:563310/140‑1 (MQ=255)
                                                                                                          agaTGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCAt                                      >  3:563310/1‑140 (MQ=255)
                                                                                                           aaTGCTGGCCCGTATTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGATTATGACGAAAGTGAAAGTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGaaaaaa                            >  4:507454/2‑149 (MQ=255)
                                                                                                                    ccGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAgg                   <  2:318972/149‑1 (MQ=255)
                                                                                                                      gTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGt                 <  1:306703/149‑1 (MQ=255)
                                                                                                                             tACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAAc          >  2:268766/1‑149 (MQ=255)
                                                                                                                             tACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAAc          <  4:111710/149‑1 (MQ=255)
                                                                                                                                aaGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTATCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACcac       >  2:59540/1‑149 (MQ=255)
                                                                                                                                 aGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACcaca      <  1:416616/149‑1 (MQ=255)
                                                                                                                                    cTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACagag   <  4:568738/149‑1 (MQ=255)
                                                                                                                                     aggAAGACGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGTCGAAAGTAAAACGCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGc  <  4:154386/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                      
GCGTGGTGGAAATCGAGAGTCGGTAGCGAGACAGAGCAAGATTCACCGTTCCTGGATTTGTATCGGGTAGCGGAACAGGTAGTACTTGAAGCGTTTCCGAAACAAGAGATGCTGGCCCGTCTTAATACAAGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGC  >  NC_000913/681301‑681582

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑