Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,376,959 C→T 100% N256N (AAC→AAT elaD → protease, capable of cleaving an AMC‑ubiquitin model substrate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,376,9590CT98.6% 668.5 / 0.2 213N256N (AAC→AATelaDprotease, capable of cleaving an AMC‑ubiquitin model substrate
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base T (104/106);  total (107/106)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 215 total aligned reads displayed.

ATGAAAATATTAAGAGCGATCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAACTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATATTTTTACATGCCGATC  >  CP009273/2376814‑2377104
                                                                                                                                                 |                                                                                                                                                 
  gAAAATATTAAGAGCGATCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCAt                                                                                                                                               >  1:2579678/1‑148 (MQ=255)
     aaTATTAAGAGCGATCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGc                                                                                                                                            <  2:2631109/148‑1 (MQ=255)
            agagCGATCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACt                                                                                                                                     <  1:1069906/148‑1 (MQ=255)
            agagCGATCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACt                                                                                                                                     <  1:463227/148‑1 (MQ=255)
            agagCGATCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACt                                                                                                                                     <  2:2714704/148‑1 (MQ=255)
            agagCGATCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCAGTGCGTGCACt                                                                                                                                     <  2:1459682/148‑1 (MQ=255)
             gagCGATCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTa                                                                                                                                    <  2:295489/148‑1 (MQ=255)
             gagCGATCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTa                                                                                                                                    <  2:2106254/148‑1 (MQ=255)
              agCGATCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAg                                                                                                                                   <  2:1245191/148‑1 (MQ=255)
              agCGATCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAg                                                                                                                                   >  2:2767013/1‑148 (MQ=255)
                  aTCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTa                                                                                                                               <  2:2900107/148‑1 (MQ=255)
                  aTCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTa                                                                                                                               >  1:591588/1‑148 (MQ=255)
                   tCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTAt                                                                                                                              >  2:512286/1‑148 (MQ=255)
                    ccATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATg                                                                                                                             >  2:1006449/1‑148 (MQ=255)
                    ccATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATg                                                                                                                             <  1:1837771/148‑1 (MQ=255)
                       tcatGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATa                                                                                                                          >  2:2421900/1‑148 (MQ=255)
                             ccGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTCTGATAAAGaac                                                                                                                    >  2:1848142/1‑147 (MQ=255)
                              cGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGaaaa                                                                                                                   <  1:1904565/148‑1 (MQ=255)
                                 aTGAAGTGCAGTTATTCCACATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTc                                                                                                                <  2:2473711/148‑1 (MQ=255)
                                  tGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCa                                                                                                               >  1:3140587/1‑148 (MQ=255)
                                  tGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCa                                                                                                               >  1:2739084/1‑148 (MQ=255)
                                  tGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCa                                                                                                               >  1:1459866/1‑148 (MQ=255)
                                   gAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCAc                                                                                                              >  1:594026/1‑148 (MQ=255)
                                     aGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACtt                                                                                                            >  2:2642529/1‑148 (MQ=255)
                                          aGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCg                                                                                                       >  1:203386/1‑148 (MQ=255)
                                            ttATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCgtg                                                                                                     <  1:2861560/148‑1 (MQ=255)
                                            ttATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCgtg                                                                                                     <  1:2419542/148‑1 (MQ=255)
                                               ttCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCa                                                                                                  >  1:2817589/1‑148 (MQ=255)
                                                tCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCAt                                                                                                 <  1:2598282/148‑1 (MQ=255)
                                                  cccATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTa                                                                                               <  2:731089/148‑1 (MQ=255)
                                                   ccATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTaa                                                                                              <  2:3111541/148‑1 (MQ=255)
                                                   ccATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTaa                                                                                              >  2:1109018/1‑148 (MQ=255)
                                                    cATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAAt                                                                                             <  2:788745/148‑1 (MQ=255)
                                                     aTAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATa                                                                                            <  1:2787070/148‑1 (MQ=255)
                                                        aaTACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGtt                                                                                         <  1:1779677/148‑1 (MQ=255)
                                                           aCTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTa                                                                                      <  1:1739035/148‑1 (MQ=255)
                                                           aCTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTa                                                                                      <  1:485156/148‑1 (MQ=255)
                                                             tGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAAt                                                                                    <  1:2636127/148‑1 (MQ=255)
                                                                cGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCt                                                                                 <  1:1046265/148‑1 (MQ=255)
                                                                   acacTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGaa                                                                              <  1:2185557/148‑1 (MQ=255)
                                                                     acTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACt                                                                            <  2:3219394/148‑1 (MQ=255)
                                                                      cTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTc                                                                           <  1:2352837/148‑1 (MQ=255)
                                                                      cTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTc                                                                           >  1:971132/1‑148 (MQ=255)
                                                                       tGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCa                                                                          <  1:363445/148‑1 (MQ=255)
                                                                             ttAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGgag                                                                    >  2:1518477/1‑148 (MQ=255)
                                                                               aGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGaa                                                                  <  2:380417/148‑1 (MQ=255)
                                                                                gTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCAGTAATAGTTTTAATTCTGAACTCATGGGAGaaa                                                                 >  1:2260122/1‑148 (MQ=255)
                                                                                  ttCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATg                                                               >  2:1744836/1‑148 (MQ=255)
                                                                                   tCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGt                                                              <  2:1813781/148‑1 (MQ=255)
                                                                                      cTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCAGTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGa                                                           <  1:228772/148‑1 (MQ=255)
                                                                                       tACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGaa                                                          <  1:3299236/148‑1 (MQ=255)
                                                                                       tACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGaa                                                          <  2:683524/148‑1 (MQ=255)
                                                                                            aaaaTAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGaataata                                                     <  2:895764/148‑1 (MQ=255)
                                                                                            aaaaTAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGaataata                                                     >  1:2805955/1‑148 (MQ=255)
                                                                                              aaTAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGaataataat                                                   <  1:11789/148‑1 (MQ=255)
                                                                                                tAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAatat                                                 >  1:3095778/1‑148 (MQ=255)
                                                                                                tAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAatat                                                 >  1:1434633/1‑148 (MQ=255)
                                                                                                  gTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATaa                                               <  2:2318209/148‑1 (MQ=255)
                                                                                                    aaataaaAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATaaaa                                             <  2:3488388/148‑1 (MQ=255)
                                                                                                    aaataaaAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATaaaa                                             <  2:2579678/148‑1 (MQ=255)
                                                                                                    aaataaaAAAAATAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATaaaa                                             <  1:379142/148‑1 (MQ=255)
                                                                                                         aaaaaaaTAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCa                                        >  1:1143537/1‑148 (MQ=255)
                                                                                                          aaaaaaTAATAAACTACAAATAAAATGCGTCATATTCAATAAATTGCGTGCACTAGGCTATGATAAATAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCAt                                       >  1:2539865/1‑148 (MQ=255)
                                                                                                           aaaaaTAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCAtt                                      <  1:1685192/148‑1 (MQ=255)
                                                                                                           aaaaaTAATAAACTACAAATAAAAAGCGTCATATACAATACATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAATCGTGTCATTAATAGTAATAATTCTTAAATAATGGGAGAAATGTCGAATAATAATATAAAAGATCAtt                                      >  1:3170717/1‑148 (MQ=255)
                                                                                                            aaaaTAATAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCAttt                                     >  2:2444620/1‑148 (MQ=255)
                                                                                                              aataatAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTaa                                   <  1:2652780/148‑1 (MQ=255)
                                                                                                              aataatAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTaa                                   >  1:1875771/1‑148 (MQ=255)
                                                                                                              aataatAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTaa                                   >  2:888816/1‑148 (MQ=255)
                                                                                                              aataatAAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTaa                                   >  2:249400/1‑148 (MQ=255)
                                                                                                               ataataAACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTaaa                                  >  2:3386424/1‑148 (MQ=255)
                                                                                                                taataaACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAAt                                 >  1:436796/1‑148 (MQ=255)
                                                                                                                taataaACTACAAATAAAATGCGTCATATTCAACTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAAt                                 >  1:1640536/1‑148 (MQ=255)
                                                                                                                  ataaACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGa                               >  2:1977379/1‑148 (MQ=255)
                                                                                                                   taaACTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGaa                              >  1:3259817/1‑148 (MQ=255)
                                                                                                                    acaCTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAAc                             <  1:3342244/146‑1 (MQ=255)
                                                                                                                    aaaCTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAAc                             >  1:3367196/1‑148 (MQ=255)
                                                                                                                     aaCTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAAcc                            >  1:2722616/1‑148 (MQ=255)
                                                                                                                     aaCTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAAcc                            <  2:1070580/148‑1 (MQ=255)
                                                                                                                      aCTACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCa                           >  2:3327066/1‑148 (MQ=255)
                                                                                                                        tACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCaga                         >  1:2464630/1‑148 (MQ=255)
                                                                                                                        tACAAATAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCaga                         <  2:2796273/148‑1 (MQ=255)
                                                                                                                            aaTAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGata                     <  2:2280315/148‑1 (MQ=255)
                                                                                                                            aaTAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGata                     <  2:836922/148‑1 (MQ=255)
                                                                                                                            aaTAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGata                     <  2:3294177/148‑1 (MQ=255)
                                                                                                                              tAAAATGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGataat                   >  1:742475/1‑148 (MQ=255)
                                                                                                                               aaaaTGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGataata                  >  2:1305535/1‑148 (MQ=255)
                                                                                                                                aaaTGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAatat                 >  2:2946730/1‑148 (MQ=255)
                                                                                                                                 aaTGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATAtt                >  2:692486/1‑148 (MQ=255)
                                                                                                                                  aagCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCTTGTCATTAATATTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATAttt               <  2:3027964/146‑1 (MQ=255)
                                                                                                                                   tGCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATAtttt              >  2:949542/1‑148 (MQ=255)
                                                                                                                                    gCGTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATAttttt             >  1:128341/1‑148 (MQ=255)
                                                                                                                                      gTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATATTTTTAc           <  2:2002656/148‑1 (MQ=255)
                                                                                                                                      gTCATATTCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATATTTTTAc           <  1:668669/148‑1 (MQ=255)
                                                                                                                                      ctCATATTAAATTCATTGGGTTCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATATTTTTAc           <  2:2622100/147‑1 (MQ=255)
                                                                                                                                             tCAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATATTTTTACATGCCGa    >  2:2159903/1‑148 (MQ=255)
                                                                                                                                              cAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATATTTTTACATGCCGAt   <  2:1542181/148‑1 (MQ=255)
                                                                                                                                              cAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATATTTTTACATGCCGAt   <  1:306780/148‑1 (MQ=255)
                                                                                                                                              cAATTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATATTTTTACATGCCGAt   >  1:289977/1‑148 (MQ=255)
                                                                                                                                               aaCTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATATTTTTACATGCCGATc  >  1:2311139/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ATGAAAATATTAAGAGCGATCCATCATGGCCGAATGAAGTCCAGTTATTCCCCATAAATACTGGCGGACACTGGATATTAGTTTCGCTACAGAAAATAGTAAATAAAAAAAATAATAAACTACAAATAAAATGCGTCATATTCAACTCATTGCGTGCACTAGGCTATGATAAAGAAAATTCACTTAAGCGTGTCATTAATAGTTTTAATTCTGAACTCATGGGAGAAATGTCGAATAATAATATAAAAGTTCATTTAAATGAACCAGAGATAATATTTTTACATGCCGATC  >  CP009273/2376814‑2377104

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