Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,454,344 C→T G254G (GGC→GGT yjfF → putative sugar ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,454,3440CT88.0% 232.3 / ‑5.3 76G254G (GGC→GGTyjfFputative sugar ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (29/37);  minor base G (9/0);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 2.29e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGCGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTT  >  NC_000913/4454196‑4454488
                                                                                                                                                       |                                                                                                                                                
ccACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTTTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGGAGGGGGGGGACTGGGTGCTATCGCCTCAATGGGAATTGgggg                                                                                                                                                >  2:455170/1‑151 (MQ=255)
 cACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTggtggt                                                                                                                                               <  1:202177/151‑1 (MQ=255)
 cACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTggtggt                                                                                                                                               <  2:430037/151‑1 (MQ=255)
  aCTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTa                                                                                                                                              <  1:502679/151‑1 (MQ=255)
     cGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACg                                                                                                                                            <  1:513448/150‑1 (MQ=255)
     cGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACg                                                                                                                                            <  1:415830/150‑1 (MQ=255)
        agcaCCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATAACCAGGCCGGGTATGCGCTGGCGGGCGGTGGGGTGGGAATGGGTGCTTTCGCCCCCGGGGGAATTTGGGGGACGCctt                                                                                                                                        >  2:472173/1‑151 (MQ=255)
            ccacTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGGAGGGGTGGGAATGGATGCTATCGCCTCCGGGGGAATTGGGGGTACGCtttt                                                                                                                                       >  2:341039/1‑148 (MQ=255)
                tATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTTTTCTCGATTTATACCCAGGCCGGATATGCGCCGGCGGGGGGTGGAGGGGGAAAGGGAGCTATCGCCTCCGTGGCAAATGGGGGTACGCcttt                                                                                                                                       >  2:923077/1‑144 (MQ=255)
                     gCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGt                                                                                                                           <  2:825314/151‑1 (MQ=255)
                      cATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGtt                                                                                                                          <  2:1142529/151‑1 (MQ=255)
                      cATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGtt                                                                                                                          <  2:384400/151‑1 (MQ=255)
                        tttATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGGGGAGGGGGGGGAATGGGTGCTAACGCCCCCGGGGGAAATTGGGGGAAGCTTTTGGGGGGGGGGGtttg                                                                                                                        >  2:512873/1‑149 (MQ=255)
                         ttATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGGGGGACTTGATGCTATCGCCCCCGTGGTAATTGGGGGGACGCTTTTTAGCGGGGGGGTTgg                                                                                                                        >  2:183039/1‑150 (MQ=255)
                           atatGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGaa                                                                                                                      <  1:573630/150‑1 (MQ=255)
                           atatGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAAc                                                                                                                     <  2:80936/151‑1 (MQ=255)
                              tGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGGGGGCGGAGGGGTGGGAACGGGTGCTTTCGCCCCACTGGGAATTGGGGGGAAGCCTTTTGGCGGGGGGGTTGGGAccg                                                                                                                   >  1:188355/1‑148 (MQ=255)
                               gCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTa                                                                                                                 <  2:349233/151‑1 (MQ=255)
                                 tctcCACCGGACTGGCAACGCTGGCGGGGGTTGTCGTCTCGCTTTATACCCAGGCCGGATATGCGCTGGCGGGGGTTGGTGGGGGAATGGGTGCTTTTGCCTCAGTGGGGGATTGTGGGACGCCTTTGAGCGGGGGCGTTGGGACGGGAtt                                                                                                               >  2:524539/1‑151 (MQ=255)
                                 tctcCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGGTATGCGCTGGCGGGGGGAGGGGTGGGAATGGGATGCAATGCCCCCGTGGGAAATTGGGGGTCGCTTTTTTGGGGGGGGGTTTGAACgggttt                                                                                                               >  2:937194/1‑147 (MQ=255)
                                 tctcCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGt                                                                                                                  <  1:207962/148‑1 (MQ=255)
                                    ccACCGGACTGGCAACGCTGGCGGGGATTGTTTTCTCGATTTATACCCCGGCCGGGGATGCGCCGGCGGTGGTAGGGGGTGAAATGGATGCCATTGCCGCAGTGGGACATTGGTGGACGCATTTGAGCCGTGGGGGTGGAACCGTCttttg                                                                                                            >  1:1058643/1‑148 (MQ=255)
                                     cACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAggg                                                                                                           <  2:1136876/151‑1 (MQ=255)
                                      aCCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGa                                                                                                          <  2:239030/151‑1 (MQ=255)
                                      aCCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGGAGGGGGGGGAATGGGTTCTATCGCCCCAGTGGGAAATGGTGGGACGCTTTTTAGCGGGGGCGTTGGAACGGGATTAgggg                                                                                                          >  1:1069677/1‑150 (MQ=255)
                                           aCTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCCGGCCGGATATGCGCTGGCGGGCGGAGGGGGGGGAATGGGTGCTATCGCCTCAGTGGGAATTGGTGGGACGCTTTTTGGCGGGGGGGTTGGGACGGGATTTGGGGCGCtt                                                                                                     >  2:1066931/1‑151 (MQ=255)
                                             tGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGGGGAGGGGGGGAAATGGATGCTATCGCCTCCGGGGGAATTTGGGGTACGCTTTTGAGCGGGGGGGGTGGGACGGGATTTGGGGCGCtttt                                                                                                   >  2:200952/1‑151 (MQ=255)
                                              ggCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCtttt                                                                                                   <  1:1075258/150‑1 (MQ=255)
                                               gCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTtg                                                                                                 <  2:188357/151‑1 (MQ=255)
                                                 aaCGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTtggc                                                                                               <  1:970194/151‑1 (MQ=255)
                                                    gCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTtggcgtg                                                                                            <  2:192624/151‑1 (MQ=255)
                                                      tGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAAATGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGGATTAGGGACGCTTTTtggcgtgg                                                                                           >  2:892153/1‑150 (MQ=255)
                                                       ggCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTtggcgtgg                                                                                           <  2:561954/149‑1 (MQ=255)
                                                        gCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGa                                                                                        <  2:1036192/151‑1 (MQ=255)
                                                          ggggATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGAtt                                                                                      <  2:60249/151‑1 (MQ=255)
                                                          ggggATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGAt                                                                                       <  2:1154185/150‑1 (MQ=255)
                                                          ggggATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGGGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGAtt                                                                                      >  2:723915/1‑151 (MQ=255)
                                                             gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTa                                                                                                                 <  1:852916/121‑1 (MQ=255)
                                                             gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTa                                                                                                                 >  2:852930/1‑121 (MQ=255)
                                                                  tcttctCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGGAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACt                                                                              >  2:1061033/1‑151 (MQ=255)
                                                                        cGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTc                                                                         <  2:1035645/150‑1 (MQ=255)
                                                                        cGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCa                                                                        <  2:787981/151‑1 (MQ=255)
                                                                           tttATACCCAGGCCGGATATGCGCTGGCGGGCGGAGGCGTGGAACTGGATGCTATCGCCTCAGGGGGAATTGGTGGGACGCTTTTGAGCGGTGGCGTTGGAACGGCATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGAccg                                                                             >  1:559617/1‑141 (MQ=255)
                                                                            ttATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGAc                                                                     <  2:898366/150‑1 (MQ=255)
                                                                              ataCCCAGGCCGGATATGCGCTGGCGGGGGTAGGTGTGGAACTGGATGCTATCGCCTCAGGGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGGATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACtt                                                                   >  1:542021/1‑150 (MQ=255)
                                                                              ataCCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTa                                                                  <  1:109701/151‑1 (MQ=255)
                                                                                 cccAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACAt                                                               >  1:851675/1‑151 (MQ=255)
                                                                                   cAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGGGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATa                                                              >  1:53324/1‑150 (MQ=255)
                                                                                   cAGGCCGGATATGCGCTGGCGGGCGTAGGGGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGGATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATa                                                              >  1:733819/1‑150 (MQ=255)
                                                                                      gCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACt                                                          <  2:812959/151‑1 (MQ=255)
                                                                                         ggATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGGAATTGGTGGGACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTg                                                       >  2:536478/1‑151 (MQ=255)
                                                                                             atGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATgg                                                   <  1:665713/151‑1 (MQ=255)
                                                                                             atGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGGACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATgg                                                   >  1:1006377/1‑151 (MQ=255)
                                                                                              tGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATgg                                                   >  2:860179/1‑150 (MQ=255)
                                                                                               gcgcTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCa                                                 <  2:518566/151‑1 (MQ=255)
                                                                                                cgcTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATgg                                                   >  1:222902/1‑148 (MQ=255)
                                                                                                     gcgggcgTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGa                                           <  1:64613/151‑1 (MQ=255)
                                                                                                       gggcgTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGt                                         <  2:1110252/151‑1 (MQ=255)
                                                                                                         gcgTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTc                                       <  2:1108245/151‑1 (MQ=255)
                                                                                                             aGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCtgg                                   >  1:405992/1‑151 (MQ=255)
                                                                                                                    ggAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACg                             <  1:505531/151‑1 (MQ=255)
                                                                                                                    ggAACTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACg                             <  2:1075381/151‑1 (MQ=255)
                                                                                                                      aaCTGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGa                            >  2:1054482/1‑150 (MQ=255)
                                                                                                                       aCTGGATGCTATCGCCTCAGTGGTAATTGGCGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGaaa                          <  2:1160918/151‑1 (MQ=255)
                                                                                                                         tGGATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAAt                        >  2:260131/1‑151 (MQ=255)
                                                                                                                           gATGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCg                      <  1:961541/151‑1 (MQ=255)
                                                                                                                             tGCTATCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGcc                    >  2:921402/1‑151 (MQ=255)
                                                                                                                                  tCGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCgg               >  1:192950/1‑151 (MQ=255)
                                                                                                                                   cGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGt              <  1:1030597/151‑1 (MQ=255)
                                                                                                                                   cGCCTCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGt              >  1:603677/1‑151 (MQ=255)
                                                                                                                                       tCAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTAttt          >  2:796030/1‑151 (MQ=255)
                                                                                                                                        cAGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTAttt          <  2:836945/150‑1 (MQ=255)
                                                                                                                                         aGTGGTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTttg        >  2:374108/1‑151 (MQ=255)
                                                                                                                                             ggTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTttgttg     <  2:3270/151‑1 (MQ=255)
                                                                                                                                              gTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTtgttgt    >  1:951390/1‑151 (MQ=255)
                                                                                                                                              gTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTtgttgt    >  2:79178/1‑151 (MQ=255)
                                                                                                                                              gTAATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTtgttgt    >  2:428801/1‑151 (MQ=255)
                                                                                                                                                  aTTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGttt  >  2:147661/1‑150 (MQ=255)
                                                                                                                                                  aTTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGttt  >  2:811562/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                
CCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGCGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTT  >  NC_000913/4454196‑4454488

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