Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,220,272 T→C 8.0% D559D (GAT→GAC aceK → isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,220,2720TC 8.0% 246.1 / 7.6 75D559D (GAT→GACaceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base T (31/38);  new base C (3/3);  total (34/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 = 9.86e-01

GATGTTTTCCCGGAAGAGTTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAATTTTTCTCCGCAATGTATCGAGGGTTATCCGTAAAGCCAAAGCTTTC  >  NC_000913/4220129‑4220396
                                                                                                                                               |                                                                                                                            
gATGTTTTCCCGGAAGAGTTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGttt                                                                                                                          <  1:27560/148‑1 (MQ=255)
 aTGTTTTCCCGGAAGAGTTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTa                                                                                                                         <  2:14816/148‑1 (MQ=255)
  tGTTTTCCCGGAAGAGTTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTAt                                                                                                                        <  4:480072/148‑1 (MQ=255)
   gTTTTCCCGGAAGAGTTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATg                                                                                                                       <  2:575513/148‑1 (MQ=255)
          cGGAAGAGTTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCg                                                                                                                <  1:159328/148‑1 (MQ=255)
                      gCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAg                                                                                                    >  1:206069/1‑148 (MQ=255)
                                 tGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGttt                                                                                                                          >  2:482479/1‑115 (MQ=255)
                                 tGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGttt                                                                                                                          <  1:482479/115‑1 (MQ=255)
                                           cgcgTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGtt                                                                                                                           <  1:105270/104‑1 (MQ=255)
                                           cgcgTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGtt                                                                                                                           >  2:105270/1‑104 (MQ=255)
                                            gcgTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGgag                                                                              <  3:16773/148‑1 (MQ=255)
                                             cgTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGAt                                                                                                                              >  2:465276/1‑99 (MQ=255)
                                             cgTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGAt                                                                                                                              <  1:465276/99‑1 (MQ=255)
                                             cgTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAg                                                                                                    <  1:548786/125‑1 (MQ=255)
                                             cgTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAg                                                                                                    >  2:548786/1‑125 (MQ=255)
                                             cgTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCTTACGGTATGGGgaga                                                                             <  4:122526/148‑1 (MQ=255)
                                                aTTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGt                                                                                                                            >  2:810619/1‑98 (MQ=255)
                                                aTTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGt                                                                                                                            <  1:810619/98‑1 (MQ=255)
                                                   ggTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCttt                                                                       <  2:206069/148‑1 (MQ=255)
                                                     tCCGCTGTTTGAAGCGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCttttt                                                                     >  4:1093454/1‑148 (MQ=255)
                                                      ccGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCaa                                                                                                      >  4:517013/1‑114 (MQ=255)
                                                      ccGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCaa                                                                                                      <  3:517013/114‑1 (MQ=255)
                                                        gCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTg                                                                   <  3:609776/147‑1 (MQ=255)
                                                        gCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTg                                                                   >  4:609776/1‑147 (MQ=255)
                                                            tttGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATc                                                                                                                 >  1:234577/1‑97 (MQ=255)
                                                            tttGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATc                                                                                                                 <  2:234577/97‑1 (MQ=255)
                                                                aaGAGATGCACGCCGACCTGTTCCGCGCTGATTGCTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCt                                                          >  1:901430/1‑148 (MQ=255)
                                                                aaGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGTAGATGt                                                                                                                            <  4:560217/82‑1 (MQ=255)
                                                                aaGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGt                                                                                                                            >  3:560217/1‑82 (MQ=255)
                                                                aaGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTa                                                                                                                         >  1:66259/1‑85 (MQ=255)
                                                                aaGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTa                                                                                                                         <  2:66259/85‑1 (MQ=255)
                                                                aaGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGAttt                                                                                                >  1:383807/1‑110 (MQ=255)
                                                                aaGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGAttt                                                                                                <  2:383807/110‑1 (MQ=255)
                                                                   agaTGCACGCCGACCTGTTCCGCGCTGATTGCTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTcc                                                       <  2:901430/148‑1 (MQ=255)
                                                                   agaTGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCttttt                                                                     >  3:79984/1‑134 (MQ=255)
                                                                   agaTGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCttttt                                                                     <  4:79984/134‑1 (MQ=255)
                                                                   agaTGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTcc                                                       <  3:1093454/148‑1 (MQ=255)
                                                                    gaTGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATgg                                                                                  >  1:449828/1‑120 (MQ=255)
                                                                    gaTGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATgg                                                                                  <  2:449828/120‑1 (MQ=255)
                                                                          cGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTc                                                        >  4:1032853/1‑140 (MQ=255)
                                                                          cGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTc                                                        <  3:1032853/140‑1 (MQ=255)
                                                                               aCCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTGTCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCAtata                                                  <  1:288928/141‑1 (MQ=255)
                                                                               aCCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGAGTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCAtata                                                  >  2:288928/1‑141 (MQ=255)
                                                                                   gTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTa                                                                                                                   >  1:998659/1‑72 (MQ=255)
                                                                                   gTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTa                                                                                                                   <  2:998659/72‑1 (MQ=255)
                                                                                    ttCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTg                                                                   >  4:443506/1‑119 (MQ=255)
                                                                                    ttCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTg                                                                   <  3:443506/119‑1 (MQ=255)
                                                                                     tCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATg                                                                                   <  1:156211/102‑1 (MQ=255)
                                                                                     tCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGAGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATg                                                                                   >  2:156211/1‑102 (MQ=255)
                                                                                         cgcTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATgggg                                                                                <  1:191076/101‑1 (MQ=255)
                                                                                         cgcTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATgggg                                                                                >  2:191076/1‑101 (MQ=255)
                                                                                             gATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTg                                                                   >  4:648340/1‑110 (MQ=255)
                                                                                             gATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTg                                                                   <  3:648340/110‑1 (MQ=255)
                                                                                               ttACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATg                                                                           <  1:349099/100‑1 (MQ=255)
                                                                                               ttACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATg                                                                           >  2:349099/1‑100 (MQ=255)
                                                                                                                aaaaCCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCt                                                          <  3:323979/100‑1 (MQ=255)
                                                                                                                aaaaCCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCt                                                          >  4:323979/1‑100 (MQ=255)
                                                                                                                   aCCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTaa                                                              <  3:753080/93‑1 (MQ=255)
                                                                                                                   aCCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTaa                                                              >  4:753080/1‑93 (MQ=255)
                                                                                                                             tGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAAtt                                               <  1:721574/98‑1 (MQ=255)
                                                                                                                             tGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAAtt                                               >  2:721574/1‑98 (MQ=255)
                                                                                                                              gAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAATTTTTCTCCGCAATGTATCGAGGGTTATCCGTAAAGCCAAAGCTTTc  <  4:597717/142‑1 (MQ=255)
                                                                                                                              gAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAATTTTTCTCCGCAATGTATCGAGGGTTATCCGTAAAGCCAAAGCTTTc  >  3:597717/1‑142 (MQ=255)
                                                                                                                              gAAGGGCATGTGGAAGACGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAAtt                                               >  3:112937/1‑97 (MQ=255)
                                                                                                                              gAAGGGCATGTGGAAGACGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAAtt                                               >  3:779560/1‑97 (MQ=255)
                                                                                                                              gAAGGGCATGTGGAAGACGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAAtt                                               <  4:779560/97‑1 (MQ=255)
                                                                                                                              gAAGGGCATGTGGAAGACGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAAtt                                               >  1:487978/1‑97 (MQ=255)
                                                                                                                              gAAGGGCATGTGGAAGACGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAAtt                                               <  4:112937/97‑1 (MQ=255)
                                                                                                                              gAAGGGCATGTGGAAGACGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAAtt                                               <  2:487978/97‑1 (MQ=255)
                                                                                                                               aaGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAATTTTTCTCCGCAATGTATCGAg                           >  1:613931/1‑116 (MQ=255)
                                                                                                                               aaGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAATTTTTCTCCGCAATGTATCGAg                           <  2:613931/116‑1 (MQ=255)
                                                                                                                                       gtgGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCttt                                                                       >  2:130022/1‑64 (MQ=255)
                                                                                                                                       gtgGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCttt                                                                       <  1:130022/64‑1 (MQ=255)
                                                                                                                                             gATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCttt                                                                       <  3:666958/58‑1 (MQ=255)
                                                                                                                                             gATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCttt                                                                       >  4:666958/1‑58 (MQ=255)
                                                                                                                                             gATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAATTTTTCTCCGCAATGTATCGAGGGTTATCCGTAAAGCCAAAGCtt    <  1:857029/125‑1 (MQ=255)
                                                                                                                                             gATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAATTTTTCTCCGCAATGTATCGAGGGTTATCCGTAAAGCCAAAGCtt    >  2:857029/1‑125 (MQ=255)
                                                                                                                                               tGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAATTTTTCTCCGCAATGTAt                               >  1:582575/1‑96 (MQ=255)
                                                                                                                                               tGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAATTTTTCTCCGCAATGTAt                               <  2:582575/96‑1 (MQ=255)
                                                                                                                                               |                                                                                                                            
GATGTTTTCCCGGAAGAGTTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAACCGCATACGTGAAGGGCATGTGGAAGATGTTTATGCGTATCGGCGCAGGCAAAGATTTAGCGTACGGTATGGGGAGATGCTTTTTTGAGTAAAGCTTCCATATAATTTTTCTCCGCAATGTATCGAGGGTTATCCGTAAAGCCAAAGCTTTC  >  NC_000913/4220129‑4220396

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