Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,037,787 A→G 6.4% V107A (GTA→GCA)  recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,037,7870AG 6.4% 286.9 / 7.5 94V107A (GTA→GCA) recJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base A (49/38);  new base G (3/3);  total (53/41)
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.87e-01

ACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTT  >  NC_000913/3037640‑3037934
                                                                                                                                                   |                                                                                                                                                   
aCCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTAc                                                                                                                                                    <  6:130361/149‑1 (MQ=255)
 ccGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTAc                                                                                                                                                    >  3:289925/1‑148 (MQ=255)
  cGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCa                                                                                                                                                  <  6:345189/149‑1 (MQ=255)
     gATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGt                                                                                                                                               <  2:293432/149‑1 (MQ=255)
          ccAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAgg                                                                                                                                                <  2:82068/143‑1 (MQ=255)
          ccAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAgg                                                                                                                                                >  1:82068/1‑143 (MQ=255)
          ccAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCg                                                                                                                                          <  4:148989/149‑1 (MQ=255)
          ccAACGAGCGAGCGTGCTCAACCCCCGCAGGGGAGGAAATACCGTTATCCACCGGGACACTTAACTGCGCGCCACGGGCATGGGCCTGCTCGCCCACTTCCGGGCTTAAGCCGTAACCGTCTTCGACACGGTGTGGTCCCAGGTAGTc                                                                                                                                           >  1:131813/1‑148 (MQ=255)
              cgagcgagCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGatat                                                                                                                                      <  2:240356/149‑1 (MQ=255)
                   gagcgTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATtgctg                                                                                                                                 >  3:344876/1‑149 (MQ=255)
                    agcgTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATtgctg                                                                                                                                 >  6:45929/1‑148 (MQ=255)
                         gCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAAcc                                                                                                                            >  5:12745/1‑148 (MQ=255)
                          cTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCaa                                                                                                                          <  5:82736/149‑1 (MQ=255)
                           tCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGTGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAag                                                                                                                         >  2:288238/1‑149 (MQ=255)
                           tCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCACCCAag                                                                                                                         >  8:367415/1‑149 (MQ=255)
                           tCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAag                                                                                                                         >  3:281127/1‑149 (MQ=255)
                           tCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAag                                                                                                                         >  1:130803/1‑149 (MQ=255)
                               cccccGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAagcgag                                                                                                                     >  2:267571/1‑149 (MQ=255)
                                ccccGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAagcgagc                                                                                                                    >  3:30561/1‑149 (MQ=255)
                                ccccGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAagcgagc                                                                                                                    >  2:78140/1‑149 (MQ=255)
                                 cccGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAagcgagcg                                                                                                                   <  6:382756/149‑1 (MQ=255)
                                  ccGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAAcc                                                                                                                            <  7:51451/139‑1 (MQ=255)
                                  ccGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAAcc                                                                                                                            >  8:51451/1‑139 (MQ=255)
                                     cATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAagcgtgctggc                                                                                                                >  6:315315/1‑144 (MQ=255)
                                         ggaggaAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATtgctgc                                                                                                                                >  6:175505/1‑128 (MQ=255)
                                         ggaggaAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATtgctgc                                                                                                                                <  5:175505/128‑1 (MQ=255)
                                         ggaggaAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCa                                                                                                           <  4:61297/149‑1 (MQ=255)
                                          gaggaAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAg                                                                                                          >  1:119914/1‑149 (MQ=255)
                                           aggaAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAg                                                                                                          >  1:91016/1‑148 (MQ=255)
                                           aggaAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGc                                                                                                         >  3:45427/1‑149 (MQ=255)
                                                  aCCGTTATCCACCGGGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGACACCAAGCGAGCGCATCGCCAGCCCGCt                                                                                                    >  6:210052/1‑147 (MQ=255)
                                                   ccGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCt                                                                                                    >  8:392181/1‑146 (MQ=255)
                                                   ccGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTa                                                                                                  <  2:119914/148‑1 (MQ=255)
                                                    cGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTa                                                                                                  >  6:234286/1‑147 (MQ=255)
                                                         tCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCCGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGt                                                                                           <  8:36313/149‑1 (MQ=255)
                                                         tCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTGCCGGGCTTAAGCCGTACCCGTGTTCGAAACGGTTTCGGACCAGGTAGTCGATATTGCTGCAAGCAAGCGAGCCCATCGCCAGCACGCGTAGCGACGt                                                                                           >  6:276673/1‑149 (MQ=255)
                                                          ccACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTg                                                                                          >  1:193683/1‑149 (MQ=255)
                                                          ccACCCTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTg                                                                                          >  5:341968/1‑149 (MQ=255)
                                                           cACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGc                                                                                         <  6:12745/149‑1 (MQ=255)
                                                           cACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGAGATTGCTGCAACCAAGCGCGCGCATCGCCAGCACGCTTAGAGACGTGc                                                                                         >  6:319524/1‑149 (MQ=255)
                                                            aCCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCt                                                                                        <  2:193683/149‑1 (MQ=255)
                                                               gTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGt                                                                                     >  7:33026/1‑149 (MQ=255)
                                                               gTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGt                                                                                     >  8:232651/1‑149 (MQ=255)
                                                                    aaTTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCgcgc                                                                                >  5:285295/1‑149 (MQ=255)
                                                                         aCTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCg                                                                           >  7:271640/1‑149 (MQ=255)
                                                                           tGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCGGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGc                                                                         >  1:171937/1‑149 (MQ=255)
                                                                           tGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTgg                                                                                      <  8:217881/136‑1 (MQ=255)
                                                                           tGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTgg                                                                                      >  7:217881/1‑136 (MQ=255)
                                                                             ctcgcCACGGGCATTGGCCTGATCGACCACTTCCGTGCTTAAGCCGTAACCGTATTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGc                                                                         <  6:142243/145‑1 (MQ=255)
                                                                               cgcCACGGGCATGGGCCTGATCGACCACTTCCGGGCTAAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCg                                                                     <  8:145421/149‑1 (MQ=255)
                                                                                gcCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCg                                                                     >  1:248489/1‑148 (MQ=255)
                                                                                gcCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGa                                                                    <  4:289925/149‑1 (MQ=255)
                                                                                 ccACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGaa                                                                   <  1:37986/149‑1 (MQ=255)
                                                                                         aTGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGa                                                           <  8:271640/149‑1 (MQ=255)
                                                                                             gCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACca                                                        >  8:84939/1‑148 (MQ=255)
                                                                                               cTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACaa                                                     <  1:267571/149‑1 (MQ=255)
                                                                                                 gATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACaat                                                    >  7:40250/1‑148 (MQ=255)
                                                                                                  aTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCCCaataa                                                  >  7:319677/1‑149 (MQ=255)
                                                                                                  aTCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACaataa                                                  <  5:283656/149‑1 (MQ=255)
                                                                                                   tCGACCACTTCCGGGCTTTAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCCGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGGGCTGGTCGCGCCGTCGGCGTCGAAATGACCGACCACaataat                                                 >  4:225536/1‑149 (MQ=255)
                                                                                                   tCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACaataat                                                 >  6:83494/1‑149 (MQ=255)
                                                                                                   tCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACaataat                                                 <  1:151543/149‑1 (MQ=255)
                                                                                                    cGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATc                                                >  1:381982/1‑149 (MQ=255)
                                                                                                      aCCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCcg                                              >  7:204368/1‑149 (MQ=255)
                                                                                                      aCCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGATGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAAAAATCcg                                              >  2:229101/1‑149 (MQ=255)
                                                                                                         aCTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGt                                           <  8:328555/149‑1 (MQ=255)
                                                                                                               gggCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTc                                     >  6:369707/1‑149 (MQ=255)
                                                                                                               gggCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTc                                     <  3:15396/149‑1 (MQ=255)
                                                                                                                     aaGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGaa                                >  3:202923/1‑148 (MQ=255)
                                                                                                                     aaGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGaa                                >  1:160962/1‑148 (MQ=255)
                                                                                                                     aaGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGaa                                >  8:162037/1‑148 (MQ=255)
                                                                                                                        ccGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGc                            <  2:70327/149‑1 (MQ=255)
                                                                                                                           tAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGtt                         <  5:64315/149‑1 (MQ=255)
                                                                                                                           tAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGtt                         >  2:307036/1‑149 (MQ=255)
                                                                                                                                gTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTaaa                    <  4:30561/149‑1 (MQ=255)
                                                                                                                                 tCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCTCGAAAAGCGTTGTAAAg                   <  7:377217/149‑1 (MQ=255)
                                                                                                                                   ttCGAAACGGTTTGGTGCCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTg                        <  4:243204/142‑1 (MQ=255)
                                                                                                                                   ttCGAAACGGTTTGGTGCCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTg                        >  3:238249/1‑142 (MQ=255)
                                                                                                                                   ttCGAAACGGTTTGGTGCCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTg                        >  3:243204/1‑142 (MQ=255)
                                                                                                                                   ttCGAAACGGTTTGGTGCCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTg                        <  4:238249/142‑1 (MQ=255)
                                                                                                                                   ttCGAAACGGTTTGGTGCCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTg                        >  5:194150/1‑142 (MQ=255)
                                                                                                                                   ttCGAAACGGTTTGGTGCCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTg                        <  6:194150/142‑1 (MQ=255)
                                                                                                                                     cGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAgg                  <  6:285295/146‑1 (MQ=255)
                                                                                                                                     cGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGa                 <  4:45427/147‑1 (MQ=255)
                                                                                                                                     cGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGAtc               >  7:36943/1‑149 (MQ=255)
                                                                                                                                      gAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGAtc               <  7:84939/148‑1 (MQ=255)
                                                                                                                                       aaaCGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGAtct              >  7:360214/1‑148 (MQ=255)
                                                                                                                                           ggTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACg         <  5:319524/149‑1 (MQ=255)
                                                                                                                                            gTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCAAAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACgg        <  2:358989/149‑1 (MQ=255)
                                                                                                                                             tttGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGc       >  3:328890/1‑149 (MQ=255)
                                                                                                                                               tGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGcc      <  5:234286/148‑1 (MQ=255)
                                                                                                                                                  tACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCtttt  <  4:262865/149‑1 (MQ=255)
                                                                                                                                                  tACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCtttt  <  2:248489/149‑1 (MQ=255)
                                                                                                                                                  tACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCtttt  >  1:45325/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTT  >  NC_000913/3037640‑3037934

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