Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,657 C→T 5.4% V224I (GTT→ATT)  umpH ← UMP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913699,6570CT 5.4% 237.5 / 2.3 74V224I (GTT→ATT) umpHUMP phosphatase
Reads supporting (aligned to +/- strand):  ref base C (31/39);  new base T (2/2);  total (33/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 = 4.42e-01

GCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTAC  >  NC_000913/699512‑699803
                                                                                                                                                 |                                                                                                                                                  
gcgcATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAacac                                                                                                                                                  <  7:361993/148‑1 (MQ=255)
gcgcATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGGCGAGCGACGAAacac                                                                                                                                                  <  1:40334/148‑1 (MQ=255)
 cgcATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGTCGTGGTTGCTTATCAGCTTACGTCGATTTCCGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCa                                                                                                                                                <  4:257378/149‑1 (MQ=255)
 cgcATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCa                                                                                                                                                <  8:43540/149‑1 (MQ=255)
   cATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGa                                                                                                                                              <  3:86995/149‑1 (MQ=255)
       aGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAg                                                                                                                                           >  6:259056/1‑148 (MQ=255)
       aGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGc                                                                                                                                          <  7:386942/149‑1 (MQ=255)
        ggCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCa                                                                                                                                         >  1:239431/1‑149 (MQ=255)
        ggCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCa                                                                                                                                         <  3:303957/149‑1 (MQ=255)
         gCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCAc                                                                                                                                        <  8:147764/149‑1 (MQ=255)
           aGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCa                                                                                                                                      <  5:259056/149‑1 (MQ=255)
             tttGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAg                                                                                                                                     >  7:64880/1‑148 (MQ=255)
               tGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGaa                                                                                                                                   <  6:172120/148‑1 (MQ=255)
                 gCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATc                                                                                                                                 <  6:213623/148‑1 (MQ=255)
                 gCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATc                                                                                                                                 <  6:206781/148‑1 (MQ=255)
                  cGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATAGTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGt                                                                                                                               <  3:288854/149‑1 (MQ=255)
                  cGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGt                                                                                                                               <  8:237970/149‑1 (MQ=255)
                      tGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGtctc                                                                                                                            <  7:185273/148‑1 (MQ=255)
                          atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAATACCAGAAAGCACCAGAATCGTCTCCAGAc                                                                                                                       >  7:46508/1‑149 (MQ=255)
                          atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAATACCAGAAAGCACCAGAATCGTCTCCAGAc                                                                                                                       >  7:52661/1‑149 (MQ=255)
                           tcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGtctc                                                                                                                            >  8:117516/1‑143 (MQ=255)
                           tcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGtctc                                                                                                                            <  7:117516/143‑1 (MQ=255)
                            caGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCCTCTCCAGAcct                                                                                                                     >  7:88751/1‑149 (MQ=255)
                               agCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGGTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCCGAATCGTCTCCAGAcctgcc                                                                                                                  >  5:36010/1‑149 (MQ=255)
                                gCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGAcctgcct                                                                                                                 >  8:188981/1‑149 (MQ=255)
                                      cACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGc                                                                                                                                          <  8:313596/118‑1 (MQ=255)
                                      cACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGc                                                                                                                                          >  7:313596/1‑118 (MQ=255)
                                         gTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCgg                                                                                                        >  3:118921/1‑149 (MQ=255)
                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGaa                                                                                                             >  6:100413/1‑143 (MQ=255)
                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGaa                                                                                                             <  5:100413/143‑1 (MQ=255)
                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGc                                                                                                       <  2:184874/149‑1 (MQ=255)
                                            cGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGc                                                                                                       >  3:362204/1‑147 (MQ=255)
                                            cGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGGAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGGCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCa                                                                                                     <  2:355363/149‑1 (MQ=255)
                                             gCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGc                                                                                                       >  5:125300/1‑146 (MQ=255)
                                               aGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGaa                                                                                                  <  3:307179/149‑1 (MQ=255)
                                               aGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGaa                                                                                                  <  2:352314/149‑1 (MQ=255)
                                                gACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGaa                                                                                                  >  4:160268/1‑148 (MQ=255)
                                                   gTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATc                                                                                              <  7:269819/149‑1 (MQ=255)
                                                      gTTGCTATTCAGATAACGTCGATTTCAGCGACCGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCgg                                                                                            >  3:2642/1‑148 (MQ=255)
                                                                 gATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTAt                                                                                >  4:172142/1‑149 (MQ=255)
                                                                    aaCGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGc                                                                             >  1:219555/1‑149 (MQ=255)
                                                                       gTCGATTTCAGCGACTGCCGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGa                                                                          >  7:323439/1‑149 (MQ=255)
                                                                        tCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGAc                                                                                                                       >  3:254108/1‑103 (MQ=255)
                                                                        tCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGAc                                                                                                                       <  4:254108/103‑1 (MQ=255)
                                                                        tCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCg                                                                              >  8:368222/1‑144 (MQ=255)
                                                                        tCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCg                                                                              <  7:368222/144‑1 (MQ=255)
                                                                            tttCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAAt                                                                      >  3:54264/1‑148 (MQ=255)
                                                                                 gCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGt                                                                <  6:125300/149‑1 (MQ=255)
                                                                                 gCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGt                                                                <  2:45916/149‑1 (MQ=255)
                                                                                               aaTCCAGCTTGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATAGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATgc                                                   <  5:99473/148‑1 (MQ=255)
                                                                                               aaTCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATgcg                                                  >  3:18485/1‑149 (MQ=255)
                                                                                                      cTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAATACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCAt                                           <  8:52661/149‑1 (MQ=255)
                                                                                                      cTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAATACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCAt                                           <  8:46508/149‑1 (MQ=255)
                                                                                                         gggCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCGAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCAttt                                         <  8:273296/148‑1 (MQ=255)
                                                                                                         gggCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCAtttt                                        <  2:239431/149‑1 (MQ=255)
                                                                                                           gCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGt                                      <  4:200085/149‑1 (MQ=255)
                                                                                                              gAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTa                                   <  4:54264/149‑1 (MQ=255)
                                                                                                               aaaGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTaa                                  <  7:188981/149‑1 (MQ=255)
                                                                                                               aaaGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTaa                                  <  4:118921/149‑1 (MQ=255)
                                                                                                               aaaGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTaa                                  <  4:2642/149‑1 (MQ=255)
                                                                                                                   gCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAAtgct                              >  7:283990/1‑149 (MQ=255)
                                                                                                                     aTACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCTCCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAAtgctgc                            >  3:95159/1‑149 (MQ=255)
                                                                                                                         tGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGcgcg                         >  3:121432/1‑148 (MQ=255)
                                                                                                                           tCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCgg                        <  4:121432/147‑1 (MQ=255)
                                                                                                                             gATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGgatga                    <  2:219555/149‑1 (MQ=255)
                                                                                                                             gATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCGGCATTTTGTTTAATGCTGCGCGgatga                    >  7:204204/1‑149 (MQ=255)
                                                                                                                              atatCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTATAATGCTGCGCGgatgat                   <  1:190919/149‑1 (MQ=255)
                                                                                                                               tatCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCTAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATc                  >  8:384628/1‑149 (MQ=255)
                                                                                                                                   gtcgAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACg              >  6:205286/1‑149 (MQ=255)
                                                                                                                                      gAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACgggc           >  2:310989/1‑149 (MQ=255)
                                                                                                                                        gCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACgggctg         <  4:362204/149‑1 (MQ=255)
                                                                                                                                          gacgaAACACCAGAAAGCACCAGGATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACgggctggg       <  1:310989/149‑1 (MQ=255)
                                                                                                                                          gacgaAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACgggctggg       >  5:83453/1‑149 (MQ=255)
                                                                                                                                            cgaAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAAtg                                <  5:358728/122‑1 (MQ=255)
                                                                                                                                            cgaAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAAtg                                >  6:358728/1‑122 (MQ=255)
                                                                                                                                             gaAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCAGTTTGTTTAATGCTGCGCGGATGATCCACgggctgggct     >  1:167596/1‑148 (MQ=255)
                                                                                                                                               atcacCAGATAGCACCAGAATCGTCTCCATACCTGCCTGGAAGCCGGCCATAATATCGGTTCGCAGGTTAGCGCCGACAATCACCGTTTCTTCCTAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTTGGCTTAc  <  6:36010/147‑1 (MQ=255)
                                                                                                                                               accacCAGAAAGCACCAGAATCGGCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTAc  <  8:370927/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTAC  >  NC_000913/699512‑699803

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