Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 147,855 A→T 5.6% V296V (GTA→GTT)  yadD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913147,8550AT 5.6% 379.7 / 8.7 107V296V (GTA→GTT) yadDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base A (46/55);  new base T (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 = 9.90e-01

Only 100 of 110 total aligned reads displayed.

GTTTGAAGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAGGCGGC  >  NC_000913/147708‑148000
                                                                                                                                                   |                                                                                                                                                 
gTTTGAAGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTaa                                                                                                                                                  <  1:403409/149‑1 (MQ=255)
 tttGAAGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAAt                                                                                                                                                 <  3:776404/149‑1 (MQ=255)
     aaGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAAc                                                                                                                                             <  3:28431/149‑1 (MQ=255)
     aaGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAAc                                                                                                                                             >  1:47760/1‑149 (MQ=255)
         gaAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTa                                                                                                                                         <  2:47760/149‑1 (MQ=255)
          aaaGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTAt                                                                                                                                        >  1:823921/1‑149 (MQ=255)
          aaaGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTAt                                                                                                                                        <  2:41077/149‑1 (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)
                      aaGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTa                                                                                                                                                   <  3:364239/126‑1 (MQ=255)
                      aaGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTa                                                                                                                                                   >  4:364239/1‑126 (MQ=255)
                        ggggATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTa                                                                                                                                         <  2:498713/134‑1 (MQ=255)
                        ggggATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTa                                                                                                                                         >  1:498713/1‑134 (MQ=255)
                            aTTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTAt                                                                                                                      >  4:1123671/1‑149 (MQ=255)
                            aTTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTAt                                                                                                                      <  1:1069581/149‑1 (MQ=255)
                            aTTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTAt                                                                                                                      <  2:149809/149‑1 (MQ=255)
                             ttCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATg                                                                                                                     >  2:1097801/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                                                                                                                                             <  1:517336/110‑1 (MQ=255)
                                            aGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAAc                                                                                                                                             >  2:517336/1‑110 (MQ=255)
                                            aGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGtt                                                                                                                        >  2:1009334/1‑131 (MQ=255)
                                            aGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGtt                                                                                                                        <  1:1009334/131‑1 (MQ=255)
                                                     gtCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATtaat                                                                                             <  1:489365/149‑1 (MQ=255)
                                                      tCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCtgtg                                                                                                                          >  2:124134/1‑119 (MQ=255)
                                                      tCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCtgtg                                                                                                                          <  1:124134/119‑1 (MQ=255)
                                                        aGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTgg                                                                                                        <  2:832431/135‑1 (MQ=255)
                                                        aGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTgg                                                                                                        >  1:832431/1‑135 (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                                                                                                                                    >  2:723276/1‑86 (MQ=255)
                                                                             tGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTaa                                                                                                                                    <  1:723276/86‑1 (MQ=255)
                                                                                 tAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTgtgt                                                                                                                         <  2:691854/93‑1 (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)
                                                                                          gTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTa                                                        <  4:593766/149‑1 (MQ=255)
                                                                                             tcGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGc                                                     >  1:706695/1‑149 (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                                                                                                 >  3:822542/1‑73 (MQ=255)
                                                                                                                             ctctTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAAt                                                                                                 <  4:822542/73‑1 (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)
                                                                                                                                            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)
                                                                                                                                                   |                                                                                                                                                 
GTTTGAAGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCTTGCTGAGATTGATAAGGTAATTAACCTTATTTAAGTTACCTGTGTTATGACAGATGACGTGGGGTAAATTAATAACTGGCGCCATCAGCCGTAGCGCCAGTTAAGTATTACGCCAGCTCGACCATTTTGTTGTCGATCAGGCGAGCATCGCCAAGCCAGGCGGC  >  NC_000913/147708‑148000

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: ‑