Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,559 A→G 11.8% intergenic (‑382/+15) 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,5590AG11.8% 195.4 / 15.6 68intergenic (‑382/+15)asnB/umpHasparagine synthetase B/UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base A (40/20);  new base G (4/4);  total (44/24)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-03

GTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGAT  >  NC_000913/699299‑699831
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                
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                                                                                                                                                                                                                                                                              >  8:260576/1‑88 (MQ=255)
                                                                                                                                                                                 atatatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACg                                                                                                                                                                                                                                                                              <  7:260576/88‑1 (MQ=255)
                                                                                                                                                                                       tGCCGGCGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCACTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGCAAGGCCTACTGTCGATATCGTCGAGCGACGAAACACCGGAACGCACCAGAATCGTCTCAAGACATGCCTGGCAGCCGGCAAGAATAGCGGTACGCCGGTTATCGCCAACAATCCCCGTTCCTTCCGAATgcgc                                                                             >  4:314199/1‑275 (MQ=255)
                                                                                                                                                                                       tGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATgcgc                                                                             >  2:73682/1‑275 (MQ=255)
                                                                                                                                                                                       tGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGGCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGACTGGAAGCCGGCCAGACTACCGGTACGCCGGTTATCGCAGACAATCACCGTTTCCTCCAAATgcac                                                                             >  4:31314/1‑273 (MQ=255)
                                                                                                                                                                                              gCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                                                                                                                                                                                                            <  7:747424/141‑1 (MQ=255)
                                                                                                                                                                                              gCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                                                                                                                                                                                                            >  8:747424/1‑141 (MQ=255)
                                                                                                                                                                                                tGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCGGAATACCGGTACGCAGGTTATCGCCGACATACACCGTTTCTTCCGAATGCGCCCGCCtttt                                                                    >  4:375011/1‑275 (MQ=255)
                                                                                                                                                                                                          tGATAAGCGTAGCCCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATa                                                                                                                                                                                                                                                             >  5:993164/1‑80 (MQ=255)
                                                                                                                                                                                                          tGATAAGCGTAGCCCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATa                                                                                                                                                                                                                                                             <  6:993164/80‑1 (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                                                                                                                                                                                                                                                  <  6:156263/38‑1 (MQ=255)
                                                                                                                                                                                                                                                               tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTc                                                                                                                                                                                                                                                  >  5:156263/1‑38 (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)
                                                                                                                                                                                                                                                                  gCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCa                                                                                                                                 >  6:421100/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                  gCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTttcttc                                                                                      >  2:454570/1‑191 (MQ=255)
                                                                                                                                                                                                                                                                   cAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGCCCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGACTGCGCCTGCATTTTGTTTAATGCTTCCCGGATTCTCCCCGGGCTGGGCTTTCCAACATAGAACGGttttgttcttttaat  >  4:428617/1‑262 (MQ=255)
                                                                                                                                                                                                                                                                   cAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGAt  >  4:98040/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                
GTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGAT  >  NC_000913/699299‑699831

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