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→GGT) yjfFputative 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

CCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGCGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTT  >  NC_000913/4454196‑4454488
                                                                                                                                                       |                                                                                                                                                
ccACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTTTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGGAGGGGG‑GGGACTGGGTGCTATCGCCTCAAT‑GGGA‑ATTGgggg                                                                                                                                                >  2:455170/1‑151 (MQ=255)
 cACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTggtggt                                                                                                                                               <  1:202177/151‑1 (MQ=255)
 cACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTggtggt                                                                                                                                               <  2:430037/151‑1 (MQ=255)
  aCTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTa                                                                                                                                              <  1:502679/151‑1 (MQ=255)
     cGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACg                                                                                                                                            <  1:513448/150‑1 (MQ=255)
     cGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACg                                                                                                                                            <  1:415830/150‑1 (MQ=255)
        agcaCCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATAACCAGGCCGGGTATGCGCTGGCGGGCGGTGGGGT‑GGGAATGGGTGCTTTCGCCCCCGG‑GGGA‑ATTTGGGGGACGCctt                                                                                                                                        >  2:472173/1‑151 (MQ=255)
            ccacTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGGAGGGGT‑GGGAATGGATGCTATCGCCTCCGG‑GGGA‑ATTGGGGGTACGCtttt                                                                                                                                       >  2:341039/1‑148 (MQ=255)
                tATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTTTTCTCGATTTATACCCAGGCCGGATATGCGCCGGCGGGGGGTGGAGG‑GGGAAAGGGAGCTATCGCCTCCGT‑GGCA‑AATGGGGGTACGCcttt                                                                                                                                       >  2:923077/1‑144 (MQ=255)
                     gCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGt                                                                                                                           <  2:825314/151‑1 (MQ=255)
                      cATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGtt                                                                                                                          <  2:1142529/151‑1 (MQ=255)
                      cATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGtt                                                                                                                          <  2:384400/151‑1 (MQ=255)
                        tttATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGGGGAGGGGG‑GGGAATGGGTGCTAACGCCCCCGG‑GGGA‑AATTGGGGGAAGCTTTTGGGGGGGGGGGtttg                                                                                                                        >  2:512873/1‑149 (MQ=255)
                         ttATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGG‑GGGACTTGATGCTATCGCCCCCGT‑GGTA‑ATTGGGGGGACGCTTTTTAGCGGGGGGGTTgg                                                                                                                        >  2:183039/1‑150 (MQ=255)
                           atatGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGaa                                                                                                                      <  1:573630/150‑1 (MQ=255)
                           atatGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAAc                                                                                                                     <  2:80936/151‑1 (MQ=255)
                              tGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGGGGGCGGAGGGGT‑GGGAACGGGTGCTTTCGCCCCACT‑GGGA‑ATTGGGGGGAAGCCTTTTGGCGGGGGGGTTGGGAccg                                                                                                                   >  1:188355/1‑148 (MQ=255)
                               gCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTa                                                                                                                 <  2:349233/151‑1 (MQ=255)
                                 tctcCACCGGACTGGCAACGCTGGCGGGGGTTGTCGTCTCGCTTTATACCCAGGCCGGATATGCGCTGGCGGGGGTTGGTGG‑GGGAATGGGTGCTTTTGCCTCAGT‑GGGG‑GATTGTGGGACGCCTTTGAGCGGGGGCGTTGGGACGGGAtt                                                                                                               >  2:524539/1‑151 (MQ=255)
                                 tctcCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGGTATGCGCTGGCGGGGGGAGGGGTGGGAATGGGATGCAAT‑GCCCCCGTGGGAA‑ATTGGGGGT‑CGCTTTTTTGGGGGGGGGTTTGAACgggttt                                                                                                               >  2:937194/1‑147 (MQ=255)
                                 tctcCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGt                                                                                                                  <  1:207962/148‑1 (MQ=255)
                                    ccACCGGACTGGCAACGCTGGCGGGGATTGTTTTCTCGATTTATACCCCGGCCGGGGATGCGCCGGCGGTGGTAGGGGG‑TGAAATGGATGCCATTGCCGCAGT‑GGGACATTGGTGG‑ACGCATTTGAGCCGTGGGGGTGGAACCGTCttttg                                                                                                            >  1:1058643/1‑148 (MQ=255)
                                     cACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAggg                                                                                                           <  2:1136876/151‑1 (MQ=255)
                                      aCCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGa                                                                                                          <  2:239030/151‑1 (MQ=255)
                                      aCCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGGAGGGGG‑GGGAATGGGTTCTATCGCCCCAGT‑GGGA‑AATGGTGGGACGCTTTTTAGCGGGGGCGTTGGAACGGGATTAgggg                                                                                                          >  1:1069677/1‑150 (MQ=255)
                                           aCTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCCGGCCGGATATGCGCTGGCGGGCGGAGGGGG‑GGGAATGGGTGCTATCGCCTCAGT‑GGGA‑ATTGGTGGGACGCTTTTTGGCGGGGGGGTTGGGACGGGATTTGGGGCGCtt                                                                                                     >  2:1066931/1‑151 (MQ=255)
                                             tGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGGGGAGGGGG‑GGAAATGGATGCTATCGCCTCCGG‑GGGA‑ATTTGGGGTACGCTTTTGAGCGGGGGGGGTGGGACGGGATTTGGGGCGCtttt                                                                                                   >  2:200952/1‑151 (MQ=255)
                                              ggCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCtttt                                                                                                   <  1:1075258/150‑1 (MQ=255)
                                               gCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTtg                                                                                                 <  2:188357/151‑1 (MQ=255)
                                                 aaCGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTtggc                                                                                               <  1:970194/151‑1 (MQ=255)
                                                    gCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTtggcgtg                                                                                            <  2:192624/151‑1 (MQ=255)
                                                      tGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAAATGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGGATTAGGGACGCTTTTtggcgtgg                                                                                           >  2:892153/1‑150 (MQ=255)
                                                       ggCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTtggcgtgg                                                                                           <  2:561954/149‑1 (MQ=255)
                                                        gCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGa                                                                                        <  2:1036192/151‑1 (MQ=255)
                                                          ggggATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGAtt                                                                                      <  2:60249/151‑1 (MQ=255)
                                                          ggggATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGAt                                                                                       <  2:1154185/150‑1 (MQ=255)
                                                          ggggATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGGGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGAtt                                                                                      >  2:723915/1‑151 (MQ=255)
                                                             gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTa                                                                                                                 <  1:852916/121‑1 (MQ=255)
                                                             gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTa                                                                                                                 >  2:852930/1‑121 (MQ=255)
                                                                  tcttctCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGGA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACt                                                                              >  2:1061033/1‑151 (MQ=255)
                                                                        cGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTc                                                                         <  2:1035645/150‑1 (MQ=255)
                                                                        cGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCa                                                                        <  2:787981/151‑1 (MQ=255)
                                                                           tttATACCCAGGCCGGATATGCGCTGGCGGGCGGAGGCGT‑GGAACTGGATGCTATCGCCTCAGG‑GGGA‑ATTGGTGGGACGCTTTTGAGCGGTGGCGTTGGAACGGCATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGAccg                                                                             >  1:559617/1‑141 (MQ=255)
                                                                            ttATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGAc                                                                     <  2:898366/150‑1 (MQ=255)
                                                                              ataCCCAGGCCGGATATGCGCTGGCGGGGGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGG‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGGATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACtt                                                                   >  1:542021/1‑150 (MQ=255)
                                                                              ataCCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTa                                                                  <  1:109701/151‑1 (MQ=255)
                                                                                 cccAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACAt                                                               >  1:851675/1‑151 (MQ=255)
                                                                                   cAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGGGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATa                                                              >  1:53324/1‑150 (MQ=255)
                                                                                   cAGGCCGGATATGCGCTGGCGGGCGTAGGGGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGGATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATa                                                              >  1:733819/1‑150 (MQ=255)
                                                                                      gCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACt                                                          <  2:812959/151‑1 (MQ=255)
                                                                                         ggATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGGA‑ATTGGTGGGACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTg                                                       >  2:536478/1‑151 (MQ=255)
                                                                                             atGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATgg                                                   <  1:665713/151‑1 (MQ=255)
                                                                                             atGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGGACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATgg                                                   >  1:1006377/1‑151 (MQ=255)
                                                                                              tGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATgg                                                   >  2:860179/1‑150 (MQ=255)
                                                                                               gcgcTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCa                                                 <  2:518566/151‑1 (MQ=255)
                                                                                                cgcTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATgg                                                   >  1:222902/1‑148 (MQ=255)
                                                                                                     gcgggcgTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGa                                           <  1:64613/151‑1 (MQ=255)
                                                                                                       gggcgTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGt                                         <  2:1110252/151‑1 (MQ=255)
                                                                                                         gcgTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTc                                       <  2:1108245/151‑1 (MQ=255)
                                                                                                             aGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCtgg                                   >  1:405992/1‑151 (MQ=255)
                                                                                                                    ggAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACg                             <  1:505531/151‑1 (MQ=255)
                                                                                                                    ggAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACg                             <  2:1075381/151‑1 (MQ=255)
                                                                                                                      aaCTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGa                            >  2:1054482/1‑150 (MQ=255)
                                                                                                                       aCTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGCGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGaaa                          <  2:1160918/151‑1 (MQ=255)
                                                                                                                         tGGATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAAt                        >  2:260131/1‑151 (MQ=255)
                                                                                                                           gATGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCg                      <  1:961541/151‑1 (MQ=255)
                                                                                                                             tGCTATCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGcc                    >  2:921402/1‑151 (MQ=255)
                                                                                                                                  tCGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCgg               >  1:192950/1‑151 (MQ=255)
                                                                                                                                   cGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGt              <  1:1030597/151‑1 (MQ=255)
                                                                                                                                   cGCCTCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGt              >  1:603677/1‑151 (MQ=255)
                                                                                                                                       tCAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTAttt          >  2:796030/1‑151 (MQ=255)
                                                                                                                                        cAGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTAttt          <  2:836945/150‑1 (MQ=255)
                                                                                                                                         aGT‑GGTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTttg        >  2:374108/1‑151 (MQ=255)
                                                                                                                                             ggTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTttgttg     <  2:3270/151‑1 (MQ=255)
                                                                                                                                              gTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTtgttgt    >  1:951390/1‑151 (MQ=255)
                                                                                                                                              gTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTtgttgt    >  2:79178/1‑151 (MQ=255)
                                                                                                                                              gTA‑ATTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTtgttgt    >  2:428801/1‑151 (MQ=255)
                                                                                                                                                  aTTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGttt  >  2:147661/1‑150 (MQ=255)
                                                                                                                                                  aTTGGTGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGttt  >  2:811562/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                
CCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGT‑GGAACTGGATGCTATCGCCTCAGT‑GGTA‑ATTGGCGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTT  >  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: ‑