Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 147,857 T→A 5.6% I297N (ATT→AAT)  yadD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913147,8570TA 5.6% 378.6 / 10.4 107I297N (ATT→AAT) yadDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base T (46/55);  new base A (3/3);  total (49/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

Only 100 of 110 total aligned reads displayed.

TTTGAAGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAGGCGGCTACC  >  NC_000913/147709‑148004
                                                                                                                                                    |                                                                                                                                                   
tttGAAGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAAt                                                                                                                                                     <  3:776404/149‑1 (MQ=255)
    aaGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAAc                                                                                                                                                 <  3:28431/149‑1 (MQ=255)
    aaGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAAc                                                                                                                                                 >  1:47760/1‑149 (MQ=255)
         aaaGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTAt                                                                                                                                            <  2:41077/149‑1 (MQ=255)
         aaaGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTAt                                                                                                                                            >  1:823921/1‑149 (MQ=255)
         aaaGGGATCGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTAt                                                                                                                                            >  1:928426/1‑149 (MQ=255)
          aaGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTTAATTACCTTAtc                                                                                                                                           <  4:893318/149‑2 (MQ=255)
          aaGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTAtt                                                                                                                                           <  4:498154/149‑1 (MQ=255)
           aGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTAttt                                                                                                                                          >  2:268295/1‑149 (MQ=255)
                ttGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGtt                                                                                                                                     >  1:149809/1‑149 (MQ=255)
                  gagaAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTAc                                                                                                                                   <  2:879309/149‑1 (MQ=255)
                  gagaAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTAc                                                                                                                                   <  3:131563/149‑1 (MQ=255)
                   agaAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTAcc                                                                                                                                  >  4:687469/1‑149 (MQ=255)
                   agaAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCCCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTAcc                                                                                                                                  <  4:751977/149‑1 (MQ=255)
                   agaAGGGGATTCAGCAGGGAAGACAGGAAGTAAGGCAGGAATTCGCCCAGCGCATTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTccc                                                                                                                                  <  2:657890/149‑1 (MQ=255)
                       ggggATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTa                                                                                                                                             >  1:498713/1‑134 (MQ=255)
                       ggggATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTa                                                                                                                                             <  2:498713/134‑1 (MQ=255)
                           aTTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTAt                                                                                                                          <  2:149809/149‑1 (MQ=255)
                           aTTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTAt                                                                                                                          <  1:1069581/149‑1 (MQ=255)
                           aTTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTAt                                                                                                                          >  4:1123671/1‑149 (MQ=255)
                             tCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGa                                                                                                                        <  3:545387/149‑1 (MQ=255)
                                gcagGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAg                                                                                                                     >  4:521314/1‑149 (MQ=255)
                                     gAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGAc                                                                                                                <  1:268295/149‑1 (MQ=255)
                                      aaGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACg                                                                                                               <  1:618113/149‑1 (MQ=255)
                                        gACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTg                                                                                                             <  3:458100/149‑1 (MQ=255)
                                         aCAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTGTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTgg                                                                                                            <  4:783217/149‑1 (MQ=255)
                                          cAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTggg                                                                                                           <  4:315670/149‑1 (MQ=255)
                                           aGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAAc                                                                                                                                                 >  2:517336/1‑110 (MQ=255)
                                           aGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAAc                                                                                                                                                 <  1:517336/110‑1 (MQ=255)
                                           aGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGtt                                                                                                                            >  2:1009334/1‑131 (MQ=255)
                                           aGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGtt                                                                                                                            <  1:1009334/131‑1 (MQ=255)
                                                gtaagtCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAAt                                                                                                     >  2:1042006/1‑149 (MQ=255)
                                                    gtCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATtaat                                                                                                 <  1:489365/149‑1 (MQ=255)
                                                     tCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCtgtg                                                                                                                              >  2:124134/1‑119 (MQ=255)
                                                     tCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCtgtg                                                                                                                              <  1:124134/119‑1 (MQ=255)
                                                       aGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTgg                                                                                                            >  1:832431/1‑135 (MQ=255)
                                                       aGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTgg                                                                                                            <  2:832431/135‑1 (MQ=255)
                                                       aGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAAc                                                                                              <  1:855914/149‑1 (MQ=255)
                                                         gAATTCGCCAAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTg                                                                                            <  1:1021067/149‑1 (MQ=255)
                                                            ttCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGgcg                                                                                         >  3:593766/1‑149 (MQ=255)
                                                                     cGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGcc                                                                                >  2:12726/1‑149 (MQ=255)
                                                                          tctGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAgc                                                                           <  3:669785/149‑1 (MQ=255)
                                                                          tctGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAgc                                                                           <  1:391100/149‑1 (MQ=255)
                                                                            tGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTaa                                                                                                                                        <  1:723276/86‑1 (MQ=255)
                                                                            tGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTaa                                                                                                                                        >  2:723276/1‑86 (MQ=255)
                                                                            tGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAgcgc                                                                         >  4:52467/1‑149 (MQ=255)
                                                                                tAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTgtgt                                                                                                                             <  2:691854/93‑1 (MQ=255)
                                                                                tAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTgtgt                                                                                                                             >  1:691854/1‑93 (MQ=255)
                                                                                   atgAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTaa                                                                  <  3:687469/147‑1 (MQ=255)
                                                                                   aGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTaa                                                                  <  4:336401/149‑1 (MQ=255)
                                                                                    ggAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACATATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTa                                                                             <  4:665419/137‑1 (MQ=255)
                                                                                    ggAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTa                                                                             >  3:665419/1‑137 (MQ=255)
                                                                                    ggAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAg                                                                 >  4:617420/1‑149 (MQ=255)
                                                                                    ggAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAg                                                                 <  1:1097801/149‑1 (MQ=255)
                                                                                     gAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGt                                                                <  3:52467/149‑1 (MQ=255)
                                                                                              gggAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGc                                                         >  1:846129/1‑147 (MQ=255)
                                                                                              gggAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGc                                                         <  2:846129/147‑1 (MQ=255)
                                                                                              gggAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCa                                                       >  1:356489/1‑149 (MQ=255)
                                                                                                   gACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGc                                                     >  3:997044/1‑146 (MQ=255)
                                                                                                   gACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCg                                                  >  3:27662/1‑149 (MQ=255)
                                                                                                      gTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATTACAGa                                                                                                                    <  2:152598/80‑1 (MQ=255)
                                                                                                      gTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGa                                                                                                                    >  1:152598/1‑80 (MQ=255)
                                                                                                      gTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGAcc                                               >  3:596011/1‑149 (MQ=255)
                                                                                                         gCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAg                                                                 <  3:469874/128‑1 (MQ=255)
                                                                                                         gCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAg                                                                 >  4:469874/1‑128 (MQ=255)
                                                                                                           agagATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTg                                                                                                             <  3:298605/82‑1 (MQ=255)
                                                                                                           agagATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTg                                                                                                             >  4:298605/1‑82 (MQ=255)
                                                                                                           agagATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCa                                                                       <  2:960510/120‑1 (MQ=255)
                                                                                                           agagATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCa                                                                       >  1:960510/1‑120 (MQ=255)
                                                                                                            gagaTGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTttg                                         >  2:19248/1‑149 (MQ=255)
                                                                                                             agaTGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTttgt                                        <  2:706695/149‑1 (MQ=255)
                                                                                                              gaTGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTttgtt                                       <  1:12726/149‑1 (MQ=255)
                                                                                                              gaTGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTttgtt                                       >  1:83068/1‑149 (MQ=255)
                                                                                                                  gCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTaa                                                                                                                                        <  3:361880/48‑1 (MQ=255)
                                                                                                                  gCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTaa                                                                                                                                        >  4:361880/1‑48 (MQ=255)
                                                                                                                  gCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCg                                   >  3:487246/1‑149 (MQ=255)
                                                                                                                           ccTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGa                          <  3:1123671/149‑1 (MQ=255)
                                                                                                                            ctctTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAAt                                                                                                     <  4:822542/73‑1 (MQ=255)
                                                                                                                            ctctTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAAt                                                                                                     >  3:822542/1‑73 (MQ=255)
                                                                                                                            ctctTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAg                         >  3:303445/1‑149 (MQ=255)
                                                                                                                             tctTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTttgtt                                       >  3:60038/1‑134 (MQ=255)
                                                                                                                             tctTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTttgtt                                       <  4:60038/134‑1 (MQ=255)
                                                                                                                                 gCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTATGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCg                    >  3:752466/1‑149 (MQ=255)
                                                                                                                                     agaTTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCaa                <  1:1042006/149‑1 (MQ=255)
                                                                                                                                      cttaaaTAAGGTTAATTACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGa                                  <  3:250076/125‑1 (MQ=255)
                                                                                                                                         aaaTAAGGTTAATTACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCa            >  4:1082909/3‑149 (MQ=255)
                                                                                                                                          gATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGccgcctggc           <  4:303445/149‑8 (MQ=255)
                                                                                                                                          gATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAg           <  4:487246/149‑1 (MQ=255)
                                                                                                                                           aTAAGGTTAATTACCTTATTTAAGTTACCTTTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTa                                                            <  4:906542/99‑1 (MQ=255)
                                                                                                                                           aTAAGGTTAATTACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTa                                                            >  3:906542/1‑99 (MQ=255)
                                                                                                                                           aTAAGGTTAATTACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAgg          >  2:441374/1‑149 (MQ=255)
                                                                                                                                           aTAAGGTAATTAACCTTATTTAAGTTACCTGTGTTAAGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTc                                    <  2:184862/123‑1 (MQ=255)
                                                                                                                                           aTAAGGTAATTAACCTTATTTAAGTTACCTGTGTTAAGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTc                                    >  1:184862/1‑123 (MQ=255)
                                                                                                                                            tAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAggc         >  3:199384/1‑149 (MQ=255)
                                                                                                                                               ggTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAg                              <  2:208064/125‑1 (MQ=255)
                                                                                                                                               ggTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAg                              >  1:208064/1‑125 (MQ=255)
                                                                                                                                               ggTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAggcggc      >  3:294050/1‑149 (MQ=255)
                                                                                                                                                   aTTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAGGCGGCTAcc  <  2:83068/149‑1 (MQ=255)
                                                                                                                                                    ttAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTAc                                                           >  4:59441/1‑91 (MQ=255)
                                                                                                                                                    ttAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTAc                                                           <  3:59441/91‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTTGAAGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAGGCGGCTACC  >  NC_000913/147709‑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: ‑