Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 894,727 C→A 5.6% A315D (GCC→GAC)  potF → putrescine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913894,7270CA 5.6% 293.9 / 3.9 85A315D (GCC→GAC) potFputrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (40/40);  new base A (2/3);  total (42/43)
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.30e-01

GATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGAC  >  NC_000913/894584‑894874
                                                                                                                                               |                                                                                                                                                   
gATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAAc                                                                                                                                                 <  8:123128/148‑1 (MQ=255)
gATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAAc                                                                                                                                                 <  7:78340/148‑1 (MQ=255)
     cAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCGATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTTGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAAt                                                                                                                                           <  4:341169/149‑1 (MQ=255)
      aaaaGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGACTTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTAGCCCAGCGCCAAt                                                                                                                                           >  8:171408/1‑148 (MQ=255)
      aaaaGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTACGCCAACGCCAAt                                                                                                                                           >  2:373890/1‑148 (MQ=255)
       aaaGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATa                                                                                                                                          >  2:34111/1‑148 (MQ=255)
           aagGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAagca                                                                                                                                     <  7:56203/149‑1 (MQ=255)
             ggggCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGGAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAagcagc                                                                                                                                   <  2:87450/149‑1 (MQ=255)
               ggCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCa                                                                                                                                 <  3:238851/149‑1 (MQ=255)
                        gTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTACGCCAACGCCAATAAAGCAGCCACGCCGCTgg                                                                                                                        <  1:373890/149‑1 (MQ=255)
                         ttCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGt                                                                                                                       <  2:265635/149‑1 (MQ=255)
                              tGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGa                                                                                                                     <  1:34111/146‑1 (MQ=255)
                                   tATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTg                                                                                                                         <  4:3472/137‑1 (MQ=255)
                                   tATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTg                                                                                                                         >  3:3472/1‑137 (MQ=255)
                                   tATTCGCCATGCCTGCGGATGCAAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGATGCGCCCGGATGTAGTACCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCgg                                                                                                               >  2:101887/1‑147 (MQ=255)
                                     ttCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGc                                                                                                                                                     >  7:234057/1‑107 (MQ=255)
                                     ttCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGc                                                                                                                                                     <  8:234057/107‑1 (MQ=255)
                                     ttCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGt                                                                                                           >  4:252375/1‑149 (MQ=255)
                                           aTGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCAGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTga                                                                                                     <  6:406338/149‑1 (MQ=255)
                                               cTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAAc                                                                                                 >  4:66039/1‑149 (MQ=255)
                                               cTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAAc                                                                                                 <  2:81144/149‑1 (MQ=255)
                                                tGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCTTATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAAc                                                                                                 >  1:64484/1‑148 (MQ=255)
                                                 gCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAAcc                                                                                                >  8:268578/1‑148 (MQ=255)
                                                  cGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGAGGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCa                                                                                              <  8:243507/149‑1 (MQ=255)
                                                   ggATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGa                                                                                                                     >  2:175017/1‑125 (MQ=255)
                                                   ggATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGa                                                                                                                     <  1:175017/125‑1 (MQ=255)
                                                    gATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAgg                                                                                            >  7:112238/1‑149 (MQ=255)
                                                        ccAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCAAACGCCAATAAAGCAGCTACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTAtt                                                                                        >  5:399084/1‑149 (MQ=255)
                                                          aaaaaCAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTa                                                                                      >  1:152482/1‑149 (MQ=255)
                                                                aaaGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGCCCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGcc                                                                                >  3:67945/1‑149 (MQ=255)
                                                                aaaGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAgg                                                                                            >  5:427429/1‑137 (MQ=255)
                                                                aaaGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAgg                                                                                            <  6:427429/137‑1 (MQ=255)
                                                                          ccTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGtt                                                                      >  5:101557/1‑149 (MQ=255)
                                                                           cTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTc                                                                     >  2:267490/1‑149 (MQ=255)
                                                                            tATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCg                                                                    >  1:353506/1‑149 (MQ=255)
                                                                             aTCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGt                                                                   >  5:271779/1‑149 (MQ=255)
                                                                             aTCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGAGTTCTATGCCAACGCCAATAAAGCAGCCACGCAGCTGGTGAGTGCGGAACTCCGTGAGAACCCAGTTAATTATCCGCCTGCGGATGTTctt                                                                   >  2:185025/1‑147 (MQ=255)
                                                                              tCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGt                                                                   <  1:267490/148‑1 (MQ=255)
                                                                              tCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTg                                                                  <  6:214803/149‑1 (MQ=255)
                                                                                aGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGc                                                                 <  4:333657/148‑1 (MQ=255)
                                                                                  ttCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAAccc                                                                                               >  7:410664/1‑116 (MQ=255)
                                                                                  ttCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAAccc                                                                                               <  8:410664/116‑1 (MQ=255)
                                                                                    ccTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAg                                                             <  2:289944/148‑1 (MQ=255)
                                                                                     cTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCt                                                           >  2:376994/1‑149 (MQ=255)
                                                                                        aaTTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGt                                                         <  2:152482/148‑1 (MQ=255)
                                                                                           tACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCAc                                                     >  2:84545/1‑149 (MQ=255)
                                                                                            aCCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCAct                                                    >  8:419969/1‑149 (MQ=255)
                                                                                             ccTGCTGCGCCCGGATGTAGTAGCGCATAATTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCAGCCTGCGGATGTTCGTGCGAAGCTGTTCActc                                                   >  6:95054/1‑149 (MQ=255)
                                                                                                    cgcCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGaaa                                              <  1:253403/147‑1 (MQ=255)
                                                                                                    cgcCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGaaa                                              <  6:365215/147‑1 (MQ=255)
                                                                                                      cccGGATGTAGTAGCGCTTATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGc                                          <  2:64484/149‑1 (MQ=255)
                                                                                                      cccGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGc                                          <  1:376994/149‑1 (MQ=255)
                                                                                                       ccGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGACAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTg                                           >  8:320979/1‑147 (MQ=255)
                                                                                                       ccGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGACAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTg                                           >  8:315011/1‑147 (MQ=255)
                                                                                                        cGGATGTCGTAGCGCATATTTCCGACCATGTGTTCTTTGCCACCGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAg                                        <  2:394835/149‑1 (MQ=255)
                                                                                                        cGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAg                                        <  8:112238/149‑1 (MQ=255)
                                                                                                        cGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAg                                        <  2:50549/149‑1 (MQ=255)
                                                                                                                gtagCGCATATTTCCGACCATGTGTTCTATGACAACGCCAATCAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGa                                 <  3:276599/148‑1 (MQ=255)
                                                                                                                   gcgcATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAAt                             <  4:67945/149‑1 (MQ=255)
                                                                                                                    cgcATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATc                            <  6:101557/149‑1 (MQ=255)
                                                                                                                      cATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGa                          >  3:20752/1‑149 (MQ=255)
                                                                                                                      cATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGa                          >  8:229000/1‑149 (MQ=255)
                                                                                                                      cATATTTCCGACCATGTGTTCTATGACAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGa                          <  7:320979/149‑1 (MQ=255)
                                                                                                                      cATATTTCCGACCATGTGTTCTATGACAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGa                          <  7:315011/149‑1 (MQ=255)
                                                                                                                       atatTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGAc                         >  2:301424/1‑149 (MQ=255)
                                                                                                                        tatTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGa                                      <  8:31716/135‑1 (MQ=255)
                                                                                                                        tatTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGa                                      >  7:31716/1‑135 (MQ=255)
                                                                                                                        tatTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGAcc                        <  3:66039/149‑1 (MQ=255)
                                                                                                                        tatTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGAcc                        >  4:52702/1‑149 (MQ=255)
                                                                                                                        tatTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGAcc                        <  7:268578/149‑1 (MQ=255)
                                                                                                                         atTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAATCGACCgt                      >  6:65849/1‑149 (MQ=255)
                                                                                                                           ttCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTc                                                       >  2:14826/1‑115 (MQ=255)
                                                                                                                           ttCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTc                                                       <  1:14826/115‑1 (MQ=255)
                                                                                                                            tCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCgtg                     <  1:96295/148‑1 (MQ=255)
                                                                                                                              cGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTgc                  <  5:95054/149‑1 (MQ=255)
                                                                                                                              cGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAATCGACCGTGTgc                  <  5:65849/148‑1 (MQ=255)
                                                                                                                              cGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCGGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGtgtg                   <  7:183029/148‑1 (MQ=255)
                                                                                                                                  cATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCAc              >  8:17296/1‑149 (MQ=255)
                                                                                                                                   aTGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCAcc             >  4:118236/1‑149 (MQ=255)
                                                                                                                                     gtgtTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTg                                           >  4:378581/1‑117 (MQ=255)
                                                                                                                                     gtgtTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTg                                           <  3:378581/117‑1 (MQ=255)
                                                                                                                                     gtgtTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCcg           >  7:140969/1‑149 (MQ=255)
                                                                                                                                          cTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGt      >  7:197593/1‑149 (MQ=255)
                                                                                                                                           tATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTg     >  4:256114/1‑149 (MQ=255)
                                                                                                                                             tGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTgg    <  3:256114/148‑1 (MQ=255)
                                                                                                                                              gCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGAc  >  3:370809/1‑149 (MQ=255)
                                                                                                                                              gCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGAc  >  7:272396/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGAC  >  NC_000913/894584‑894874

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: