Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,556 C→T 12.5% intergenic (‑379/+18) asnB ← / ← umpH asparagine synthetase B/UMP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913699,5560CT12.5% 170.7 / 16.8 64intergenic (‑379/+18)asnB/umpHasparagine synthetase B/UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base C (36/20);  new base T (4/4);  total (40/24)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-03

GTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCG  >  NC_000913/699299‑699826
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                              
gTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCTTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGGCATCAGAGCCAAACACGTCCGCATCAGTGGTTGcctt                                                                                                                                                                                                                                                               >  4:367376/1‑272 (MQ=255)
                                                                                             ggCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACg                                                                                                                                                                                                                                                                         >  4:73672/1‑172 (MQ=255)
                                                                                              gCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTgggg                                                                                                                                                                                                                 >  4:250640/1‑227 (MQ=255)
                                                                                                 aaaTTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGt                                                                                                                                                                                                                                                                     <  2:73673/172‑1 (MQ=255)
                                                                                                  aaTTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGa                                                                                                                                                                                                             <  2:250641/227‑1 (MQ=255)
                                                                                                            aaaGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATCTCGTCGAGCGACGAAACACCAGAAAGCACCAGACTCCtctc                                                                                                                                                   >  2:232366/1‑275 (MQ=255)
                                                                                                                                  gATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGc                                                                                                                              >  2:301958/1‑274 (MQ=255)
                                                                                                                                      cggcgTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGaa                                                                                                                         <  2:367377/275‑1 (MQ=255)
                                                                                                                                                     aTCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTc                                                                                                                                                                                                                                                             >  6:152191/1‑128 (MQ=255)
                                                                                                                                                     aTCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTc                                                                                                                                                                                                                                                             <  5:152191/128‑1 (MQ=255)
                                                                                                                                                          gCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAg                                                                                                                                                            >  4:441345/1‑220 (MQ=255)
                                                                                                                                                              aCAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATc                                                                                                                                                        <  2:441346/220‑1 (MQ=255)
                                                                                                                                                                     cATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCAc                                                                                          >  2:366872/1‑275 (MQ=255)
                                                                                                                                                                                 atatatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACg                                                                                                                                                                                                                                                                         <  7:260576/88‑1 (MQ=255)
                                                                                                                                                                                 atatatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACg                                                                                                                                                                                                                                                                         >  8:260576/1‑88 (MQ=255)
                                                                                                                                                                                       tGCCGGCGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCACTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGCAAGGCCTACTGTCGATATCGTCGAGCGACGAAACACCGGAACGCACCAGAATCGTCTCAAGACATGCCTGGCAGCCGGCAAGAATAGCGGTACGCCGGTTATCGCCAACAATCCCCGTTCCTTCCGAATgcgc                                                                        >  4:314199/1‑275 (MQ=255)
                                                                                                                                                                                       tGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATgcgc                                                                        >  2:73682/1‑275 (MQ=255)
                                                                                                                                                                                       tGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGGCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGACTGGAAGCCGGCCAGACTACCGGTACGCCGGTTATCGCAGACAATCACCGTTTCCTCCAAATgcac                                                                        >  4:31314/1‑273 (MQ=255)
                                                                                                                                                                                              gCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                                                                                                                                                                                                       >  8:747424/1‑141 (MQ=255)
                                                                                                                                                                                              gCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                                                                                                                                                                                                       <  7:747424/141‑1 (MQ=255)
                                                                                                                                                                                                tGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCGGAATACCGGTACGCAGGTTATCGCCGACATACACCGTTTCTTCCGAATGCGCCCGCCtttt                                                               >  4:375011/1‑275 (MQ=255)
                                                                                                                                                                                                          tGATAAGCGTAGCCCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATa                                                                                                                                                                                                                                                        <  6:993164/80‑1 (MQ=255)
                                                                                                                                                                                                          tGATAAGCGTAGCCCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATa                                                                                                                                                                                                                                                        >  5:993164/1‑80 (MQ=255)
                                                                                                                                                                                                           gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTgg                                                                                                                                                                                                                                                                      <  5:201512/65‑1 (MQ=21)
                                                                                                                                                                                                           gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTgg                                                                                                                                                                                                                                                                      >  6:201512/1‑65 (MQ=21)
                                                                                                                                                                                                           gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGt                                                                                                                                                                                                                                                                     <  7:219070/66‑1 (MQ=255)
                                                                                                                                                                                                           gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGt                                                                                                                                                                                                                                                                     >  8:219070/1‑66 (MQ=21)
                                                                                                                                                                                                           gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGa                                                                                                                                                                                                                                        <  6:230286/95‑1 (MQ=255)
                                                                                                                                                                                                           gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGa                                                                                                                                                                                                                                        >  5:230286/1‑95 (MQ=255)
                                                                                                                                                                                                            aTAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTCTTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAACGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAAtgctgc                                                   >  2:96423/1‑275 (MQ=255)
                                                                                                                                                                                                            aTAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCg                                                                                                                                                                                                                                                  <  6:622875/84‑1 (MQ=255)
                                                                                                                                                                                                            aTAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCg                                                                                                                                                                                                                                                  >  5:622875/1‑84 (MQ=255)
                                                                                                                                                                                                            aTAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAgg                                                                                                                                                                                                         >  1:138358/1‑125 (MQ=255)
                                                                                                                                                                                                                gCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAgg                                                                                                                                                                                                         <  7:777557/121‑1 (MQ=255)
                                                                                                                                                                                                                gCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAgg                                                                                                                                                                                                         >  8:777557/1‑121 (MQ=255)
                                                                                                                                                                                                                gCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATa                                                                                                                                                                                                     <  3:138358/125‑1 (MQ=255)
                                                                                                                                                                                                                   tAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGgatg                                            >  1:18425/1‑275 (MQ=255)
                                                                                                                                                                                                                      cgcATCAGGCAGTTTGGCTTTTGCCATCAGAGCTAACCACGTCCGCAGACGTGGTTGCTATTCAGACAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCCGAAAGGAATACTTTCGATATCGTCGAGCGACGAAACACCACAAAGCACCAGAATCGTCTCCAGACCCGCCTGGCAGCCGGCCAGAATCCGCATACGCCGGTCATCGCCGAGAATCACCGTTTCTCCCGAATGCGCCTGCCATTTGTTTCATGTTGACCGgctgaac                                         >  4:140395/1‑273 (MQ=255)
                                                                                                                                                                                                                      cgcATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATGTTTTTTAATGCAGCGCCgttgaca                                         >  4:176794/1‑273 (MQ=255)
                                                                                                                                                                                                                           cAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCCCGGATTATCCACgg                                    >  2:357530/1‑275 (MQ=255)
                                                                                                                                                                                                                             ggCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCa                                                                                                                                                                                                                                                            <  8:724921/57‑1 (MQ=255)
                                                                                                                                                                                                                             ggCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCa                                                                                                                                                                                                                                                            >  7:724921/1‑57 (MQ=255)
                                                                                                                                                                                                                               cAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                                                                                                                                                                                                       >  8:905953/1‑108 (MQ=255)
                                                                                                                                                                                                                               cAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                                                                                                                                                                                                       <  7:905953/108‑1 (MQ=255)
                                                                                                                                                                                                                                            gTCATCAGAGTCAACCACGTCCTCAGACGTGGTTTCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGCGCGACGAACCACCAGAAAGCACCACAATCGTCTCCAGACCTGCCTTGAAGCCGGCAAGAATATCGGTACTCAGGTTATCGCCGACAATCCCCGTTTCTTCCGAATGCCCCTCCCTTTTGTTTAATGCTGCCCGGAGGACCCCCGGGCCGGCCTTGCCaaaa                    >  4:215847/1‑272 (MQ=255)
                                                                                                                                                                                                                                            gTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATACCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCTCCTGCATTTTGTTTAATGCTGCGCGGATGATCCCCGGGCTGGGCTTACCATCa                    >  2:469020/1‑274 (MQ=255)
                                                                                                                                                                                                                                                       cAACCACGTCCGTAGGCGTGGTTGCTATTCAGAGAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGGGCGACGAAACACCAGGAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGGACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTg        <  2:215848/275‑1 (MQ=255)
                                                                                                                                                                                                                                                         aCCACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGa                                                                                                                                     >  7:709500/1‑148 (MQ=255)
                                                                                                                                                                                                                                                          ccACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGatat                                                                                                                                                                                           >  7:1080641/1‑93 (MQ=255)
                                                                                                                                                                                                                                                          ccACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGatat                                                                                                                                                                                           <  8:1080641/93‑1 (MQ=255)
                                                                                                                                                                                                                                                          ccACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtcg                                                                                                                                                                                      >  5:1105047/1‑98 (MQ=255)
                                                                                                                                                                                                                                                          ccACGTCCGCAGACGTGGTTGCTAATCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtcg                                                                                                                                                                                      <  6:1105047/98‑1 (MQ=255)
                                                                                                                                                                                                                                                           cACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAg                                                                                                                                                                                                                                           >  6:951294/1‑44 (MQ=255)
                                                                                                                                                                                                                                                           cACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAg                                                                                                                                                                                                                                           <  5:951294/44‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGc                                                                                                                                  <  8:709500/148‑1 (MQ=255)
                                                                                                                                                                                                                                                              gTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACCTAGAACGGTTTGCGCCCg  >  2:57594/1‑274 (MQ=255)
                                                                                                                                                                                                                                                              gTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCg  >  2:394313/1‑274 (MQ=255)
                                                                                                                                                                                                                                                              gTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCCCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCTCGGATGATCCACGGGCTGGACTTACCAACATAGAACGGTTTGCGCCCg  >  2:408942/1‑274 (MQ=255)
                                                                                                                                                                                                                                                               tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTc                                                                                                                                                                                                                                             >  5:156263/1‑38 (MQ=255)
                                                                                                                                                                                                                                                               tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTc                                                                                                                                                                                                                                             <  6:156263/38‑1 (MQ=255)
                                                                                                                                                                                                                                                               tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac                                                                                                                                                                           <  7:776809/104‑1 (MQ=255)
                                                                                                                                                                                                                                                               tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac                                                                                                                                                                           >  8:776809/1‑104 (MQ=255)
                                                                                                                                                                                                                                                               tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGt                                                                                                                                                      >  8:447288/1‑125 (MQ=255)
                                                                                                                                                                                                                                                               tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGt                                                                                                                                                      <  7:447288/125‑1 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                              
GTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCG  >  NC_000913/699299‑699826

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