Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 147,859 A→T 5.3% N298Y (AAC→TAC)  yadD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913147,8590AT 5.3% 390.0 / 7.9 111N298Y (AAC→TAC) yadDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base A (47/57);  new base T (3/3);  total (51/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

AAGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAGGCGGCTACC  >  NC_000913/147713‑148004
                                                                                                                                                   |                                                                                                                                                 
aaGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAAc                                                                                                                                                 >  1:47760/1‑149 (MQ=255)
aaGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAAc                                                                                                                                                 <  3:28431/149‑1 (MQ=255)
    gaAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTa                                                                                                                                             <  2:47760/149‑1 (MQ=255)
     aaaGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTAt                                                                                                                                            <  2:41077/149‑1 (MQ=255)
     aaaGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTAt                                                                                                                                            >  1:823921/1‑149 (MQ=255)
     aaaGGGATCGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTAt                                                                                                                                            >  1:928426/1‑149 (MQ=255)
      aaGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TTAATTACCTTAtc                                                                                                                                           <  4:893318/149‑2 (MQ=255)
      aaGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTAtt                                                                                                                                           <  4:498154/149‑1 (MQ=255)
         ggATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTaa                                                                                                                                        <  2:928426/149‑1 (MQ=255)
            ttGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGtt                                                                                                                                     >  1:149809/1‑149 (MQ=255)
              gagaAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTAc                                                                                                                                   <  2:879309/149‑1 (MQ=255)
              gagaAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTAc                                                                                                                                   <  3:131563/149‑1 (MQ=255)
               agaAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTAcc                                                                                                                                  >  4:687469/1‑149 (MQ=255)
               agaAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCCCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTAcc                                                                                                                                  <  4:751977/149‑1 (MQ=255)
               agaAGGGGATTCAGCAGGGAAGACAGGAAGTAAGGCAGGAATTCGCCCAGCGCATTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTccc                                                                                                                                  <  2:657890/149‑1 (MQ=255)
                   ggggATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTa                                                                                                                                             <  2:498713/134‑1 (MQ=255)
                   ggggATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTa                                                                                                                                             >  1:498713/1‑134 (MQ=255)
                       aTTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTAt                                                                                                                          >  4:1123671/1‑149 (MQ=255)
                       aTTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTAt                                                                                                                          <  2:149809/149‑1 (MQ=255)
                       aTTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTAt                                                                                                                          <  1:1069581/149‑1 (MQ=255)
                        ttCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATg                                                                                                                         >  2:1097801/1‑149 (MQ=255)
                         tCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGa                                                                                                                        <  3:545387/149‑1 (MQ=255)
                            gcagGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAg                                                                                                                     >  4:521314/1‑149 (MQ=255)
                                 gAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGAc                                                                                                                <  1:268295/149‑1 (MQ=255)
                                  aaGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACg                                                                                                               <  1:618113/149‑1 (MQ=255)
                                    gACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTg                                                                                                             <  3:458100/149‑1 (MQ=255)
                                     aCAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTGTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTgg                                                                                                            <  4:783217/149‑1 (MQ=255)
                                       aGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAAc                                                                                                                                                 <  1:517336/110‑1 (MQ=255)
                                       aGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAAc                                                                                                                                                 >  2:517336/1‑110 (MQ=255)
                                       aGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGtt                                                                                                                            <  1:1009334/131‑1 (MQ=255)
                                            gtaagtCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAAt                                                                                                     >  2:1042006/1‑149 (MQ=255)
                                                gtCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATtaat                                                                                                 <  1:489365/149‑1 (MQ=255)
                                                 tCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCtgtg                                                                                                                              >  2:124134/1‑119 (MQ=255)
                                                 tCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCtgtg                                                                                                                              <  1:124134/119‑1 (MQ=255)
                                                   aGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTgg                                                                                                            >  1:832431/1‑135 (MQ=255)
                                                   aGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTgg                                                                                                            <  2:832431/135‑1 (MQ=255)
                                                   aGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAAc                                                                                              <  1:855914/149‑1 (MQ=255)
                                                     gAATTCGCCAAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTg                                                                                            <  1:1021067/149‑1 (MQ=255)
                                                        ttCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGgcg                                                                                         >  3:593766/1‑149 (MQ=255)
                                                                 cGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGcc                                                                                >  2:12726/1‑149 (MQ=255)
                                                                      tctGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAgc                                                                           <  3:669785/149‑1 (MQ=255)
                                                                      tctGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAgc                                                                           <  1:391100/149‑1 (MQ=255)
                                                                        tGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTaa                                                                                                                                        >  2:723276/1‑86 (MQ=255)
                                                                        tGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTaa                                                                                                                                        <  1:723276/86‑1 (MQ=255)
                                                                        tGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAgcgc                                                                         >  4:52467/1‑149 (MQ=255)
                                                                            tAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTgtgt                                                                                                                             <  2:691854/93‑1 (MQ=255)
                                                                              aaGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTa                                                                   >  3:336401/1‑149 (MQ=255)
                                                                               aGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTaa                                                                  <  4:336401/149‑1 (MQ=255)
                                                                                ggAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACATATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTa                                                                             <  4:665419/137‑1 (MQ=255)
                                                                                ggAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTa                                                                             >  3:665419/1‑137 (MQ=255)
                                                                                ggAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAg                                                                 >  4:617420/1‑149 (MQ=255)
                                                                                ggAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAg                                                                 <  1:1097801/149‑1 (MQ=255)
                                                                                     gTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTa                                                            <  4:593766/149‑1 (MQ=255)
                                                                                        tcGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGc                                                         >  1:706695/1‑149 (MQ=255)
                                                                                          gggAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGc                                                         <  2:846129/147‑1 (MQ=255)
                                                                                          gggAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGc                                                         >  1:846129/1‑147 (MQ=255)
                                                                                          gggAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCa                                                       >  1:356489/1‑149 (MQ=255)
                                                                                               gACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGc                                                     >  3:997044/1‑146 (MQ=255)
                                                                                                  gTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATTACAGa                                                                                                                    <  2:152598/80‑1 (MQ=255)
                                                                                                  gTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGa                                                                                                                    >  1:152598/1‑80 (MQ=255)
                                                                                                  gTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGAcc                                               >  3:596011/1‑149 (MQ=255)
                                                                                                     gCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAg                                                                 >  4:469874/1‑128 (MQ=255)
                                                                                                     gCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAg                                                                 <  3:469874/128‑1 (MQ=255)
                                                                                                       agagATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTg                                                                                                             >  4:298605/1‑82 (MQ=255)
                                                                                                       agagATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTg                                                                                                             <  3:298605/82‑1 (MQ=255)
                                                                                                       agagATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCa                                                                       <  2:960510/120‑1 (MQ=255)
                                                                                                       agagATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCa                                                                       >  1:960510/1‑120 (MQ=255)
                                                                                                        gagaTGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTttg                                         >  2:19248/1‑149 (MQ=255)
                                                                                                         agaTGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTttgt                                        <  2:706695/149‑1 (MQ=255)
                                                                                                          gaTGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTttgtt                                       <  3:617420/149‑1 (MQ=255)
                                                                                                          gaTGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTttgtt                                       >  1:83068/1‑149 (MQ=255)
                                                                                                          gaTGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTttgtt                                       <  1:12726/149‑1 (MQ=255)
                                                                                                              gCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTaa                                                                                                                                        >  4:361880/1‑48 (MQ=255)
                                                                                                              gCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTaa                                                                                                                                        <  3:361880/48‑1 (MQ=255)
                                                                                                              gCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCg                                   >  3:487246/1‑149 (MQ=255)
                                                                                                                       ccTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGa                          <  3:1123671/149‑1 (MQ=255)
                                                                                                                        ctctTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAAt                                                                                                     <  4:822542/73‑1 (MQ=255)
                                                                                                                        ctctTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAAt                                                                                                     >  3:822542/1‑73 (MQ=255)
                                                                                                                        ctctTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAg                         >  3:303445/1‑149 (MQ=255)
                                                                                                                         tctTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTttgtt                                       <  4:60038/134‑1 (MQ=255)
                                                                                                                         tctTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTttgtt                                       >  3:60038/1‑134 (MQ=255)
                                                                                                                             gCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTATGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCg                    >  3:752466/1‑149 (MQ=255)
                                                                                                                                 agaTTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCaa                <  1:1042006/149‑1 (MQ=255)
                                                                                                                                  cttaaaTAAGGTTAATT‑ACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGa                                  >  4:250076/6‑130 (MQ=255)
                                                                                                                                  cttaaaTAAGG‑TTAATTACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGa                                  <  3:250076/125‑1 (MQ=255)
                                                                                                                                     aaaTAAGG‑TTAATTACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCa            >  4:1082909/3‑149 (MQ=255)
                                                                                                                                      gATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGccgcctggc           <  4:303445/149‑8 (MQ=255)
                                                                                                                                      gATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAg           <  4:487246/149‑1 (MQ=255)
                                                                                                                                       aTAAGG‑TTAATTACCTTATTTAAGTTACCTTTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTa                                                            <  4:906542/99‑1 (MQ=255)
                                                                                                                                       aTAAGG‑TTAATTACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTa                                                            >  3:906542/1‑99 (MQ=255)
                                                                                                                                       aTAAGG‑TTAATTACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAgg          >  2:441374/1‑149 (MQ=255)
                                                                                                                                       aTAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTAAGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTc                                    <  2:184862/123‑1 (MQ=255)
                                                                                                                                       aTAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTAAGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTc                                    >  1:184862/1‑123 (MQ=255)
                                                                                                                                        tAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAggc         >  3:199384/1‑149 (MQ=255)
                                                                                                                                           gg‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAg                              <  2:208064/125‑1 (MQ=255)
                                                                                                                                           gg‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAg                              >  1:208064/1‑125 (MQ=255)
                                                                                                                                           gg‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAggcggc      >  3:294050/1‑149 (MQ=255)
                                                                                                                                                 ttAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTAc                                                           <  3:59441/91‑1 (MQ=255)
                                                                                                                                                  tAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGc                                                         >  4:530141/1‑92 (MQ=255)
                                                                                                                                                  tAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGc                                                         <  3:530141/92‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AAGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGG‑TAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAGGCGGCTACC  >  NC_000913/147713‑148004

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