Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,256,442 G→C 7.3% G70G (GGC→GGG) ‡ yhaM ← putative L‑serine dehydratase alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,256,4420GC 7.3% 308.7 / 10.6 96G70G (GGC→GGG) ‡yhaMputative L‑serine dehydratase alpha chain
Reads supporting (aligned to +/- strand):  ref base G (46/42);  new base C (3/4);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTA  >  NC_000913/3256294‑3256590
                                                                                                                                                    |                                                                                                                                                     
tCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTg                                                                                                                                                       >  1:476596/1‑149 (MQ=255)
  ttAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCACCGCAATCGGCAGCCCCACCATTCCCGTGc‑c                                                                                                                                                    <  1:833873/149‑1 (MQ=255)
         agaCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGGCA‑CGGGaatg                                                                                                                                             <  2:827614/149‑3 (MQ=255)
           aCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCGCTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCGTTCCCGGCA‑CGGGAATGGt                                                                                                                                           <  2:35919/149‑1 (MQ=255)
           aCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCGTTCCCGGCA‑CGGGAATGGt                                                                                                                                           <  4:851150/149‑1 (MQ=255)
           aCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGt                                                                                                                                           <  2:752219/149‑1 (MQ=255)
           aCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGt                                                                                                                                           <  3:595660/149‑1 (MQ=255)
            cTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGGCA‑CGGGAATGGTg                                                                                                                                          <  3:61863/149‑1 (MQ=255)
            cTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCG‑GCACGGGAATGGTg                                                                                                                                          <  1:115706/149‑1 (MQ=255)
             tttCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGa                                                                                                                                         >  4:683247/1‑149 (MQ=255)
              ttCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGAc                                                                                                                                        <  3:683247/149‑1 (MQ=255)
                  ccgccgCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGTCCGGCGTTGGCATTGCCACCTAACGCCCCCAGCGCCGCCGCAATCGGTAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGccc                                                                                                                                    >  1:558574/1‑149 (MQ=255)
                  ccgccgCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGccc                                                                                                                                    >  1:369576/1‑149 (MQ=255)
                   cgccgcCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCa                                                                                                                                   >  3:507213/1‑149 (MQ=255)
                   cgccgcCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCa                                                                                                                                   <  3:195140/149‑1 (MQ=255)
                          agcagTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGtt                                                                                                                            >  4:208002/1‑149 (MQ=255)
                           gcagTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGttc                                                                                                                           <  1:871059/149‑1 (MQ=255)
                                 cTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTtcatc                                                                                                                     <  3:208002/149‑1 (MQ=255)
                                 cTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTtcatc                                                                                                                     <  2:558574/149‑1 (MQ=255)
                                  tttGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTtcatca                                                                                                                    >  4:858470/1‑149 (MQ=255)
                                     ggCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGAt                                                                                                                 <  2:442428/149‑1 (MQ=255)
                                     ggCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGAt                                                                                                                 <  2:667887/149‑1 (MQ=255)
                                      gCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGAtt                                                                                                                <  1:722561/149‑1 (MQ=255)
                                      gCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGAtt                                                                                                                >  1:165903/1‑149 (MQ=255)
                                        aTCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTg                                                                                                              >  1:436513/1‑149 (MQ=255)
                                        aTCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTg                                                                                                              <  4:132075/149‑1 (MQ=255)
                                         tCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTgg                                                                                                             >  1:163656/1‑149 (MQ=255)
                                         tCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTgg                                                                                                             <  2:108228/149‑1 (MQ=255)
                                         tCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTgg                                                                                                             >  2:63485/1‑149 (MQ=255)
                                         tCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTgg                                                                                                             >  2:722403/1‑149 (MQ=255)
                                           ggCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGACg                                                                                                           >  3:702757/1‑149 (MQ=255)
                                            gCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGa                                                                                                          <  2:676645/149‑1 (MQ=255)
                                            gCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGa                                                                                                          >  2:811100/1‑149 (MQ=255)
                                                 tGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAAccc                                                                                                     <  2:584529/149‑1 (MQ=255)
                                                  gCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCa                                                                                                    <  3:677853/149‑1 (MQ=255)
                                                        gcTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTc                                                                                              <  3:545469/149‑1 (MQ=255)
                                                          tGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTa                                                                                            <  4:686499/149‑1 (MQ=255)
                                                            ttGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTaca                                                                                          <  3:108326/149‑1 (MQ=255)
                                                              gCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATGTGGCGAAACCCAGGCTTCTACACg                                                                                        >  2:584809/1‑149 (MQ=255)
                                                                 tCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTc                                                                                     <  3:843163/149‑1 (MQ=255)
                                                                 tCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTc                                                                                     <  4:457982/149‑1 (MQ=255)
                                                                 tCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGAAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCTCCGACCGTTCTTCATCAGATTTGGCGAAACCTAGGCTTCTACACGTTc                                                                                     >  2:313489/1‑149 (MQ=255)
                                                                  cTTTCAGCACTTCCAGCCCCGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCa                                                                                    >  4:852825/1‑149 (MQ=255)
                                                                   tttCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCaa                                                                                   <  2:243068/149‑1 (MQ=255)
                                                                       aGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCG‑GCACGGGAATGGTg                                                                                                                                          >  2:429887/1‑90 (MQ=255)
                                                                       aGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCG‑GCACGGGAATGGTg                                                                                                                                          <  1:429887/90‑1 (MQ=255)
                                                                             tCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTc                                                                         <  1:886776/149‑1 (MQ=255)
                                                                                aGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGGTTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAg                                                                      >  2:607253/1‑149 (MQ=255)
                                                                                       cGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGcc                                                               <  2:374502/149‑1 (MQ=255)
                                                                                       cGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGcc                                                               <  3:858470/149‑1 (MQ=255)
                                                                                         ttGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAAcc                                                                                 >  3:812128/1‑129 (MQ=255)
                                                                                         ttGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAAcc                                                                                 <  4:812128/129‑1 (MQ=255)
                                                                                         ttGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCgc                                                             >  3:271169/1‑149 (MQ=255)
                                                                                          tGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCgca                                                            <  2:163656/149‑1 (MQ=255)
                                                                                            gCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCgcagc                                                          >  2:884291/1‑149 (MQ=255)
                                                                                              aTTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCaa                                                        >  2:25208/1‑149 (MQ=255)
                                                                                               tttCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGtt                                                                    >  1:536534/1‑136 (MQ=255)
                                                                                               tttCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGtt                                                                    <  2:536534/136‑1 (MQ=255)
                                                                                               tttCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAAc                                                       >  4:51835/1‑149 (MQ=255)
                                                                                               tttCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAAc                                                       >  4:141919/1‑149 (MQ=255)
                                                                                                    aCCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGTTTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAAccgccg                                                  >  3:234895/1‑149 (MQ=255)
                                                                                                      cTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAAccgccgcc                                                >  2:467048/1‑149 (MQ=255)
                                                                                                          cGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgcca                                            >  1:817376/1‑149 (MQ=255)
                                                                                                           gCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccag                                           <  3:85384/149‑1 (MQ=255)
                                                                                                              cccAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagcgc                                        >  1:219430/1‑149 (MQ=255)
                                                                                                              cccAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagcgc                                        >  2:701681/1‑149 (MQ=255)
                                                                                                               ccAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTa                                                                                            <  4:813314/96‑1 (MQ=255)
                                                                                                               ccAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTa                                                                                            >  3:813314/1‑96 (MQ=255)
                                                                                                               ccAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagcgcc                                       >  4:388683/1‑149 (MQ=255)
                                                                                                                 aGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagcgccag                                     >  2:874780/1‑149 (MQ=255)
                                                                                                                 aGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagcgccag                                     <  2:473990/149‑1 (MQ=255)
                                                                                                                       gccgcAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAAt                               >  1:90048/1‑149 (MQ=255)
                                                                                                                        tcgcAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATc                              <  1:313489/148‑1 (MQ=255)
                                                                                                                        ccgcAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATc                              >  4:37275/1‑149 (MQ=255)
                                                                                                                         cgcAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCg                             >  2:506427/1‑149 (MQ=255)
                                                                                                                         cgcAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCg                             >  1:630752/1‑149 (MQ=255)
                                                                                                                         cgcAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCg                             <  1:584809/149‑1 (MQ=255)
                                                                                                                            aaTCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGtt                          <  4:820158/149‑1 (MQ=255)
                                                                                                                               cGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGtt                          <  1:242175/146‑1 (MQ=255)
                                                                                                                               cGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGtt                          >  2:242175/1‑146 (MQ=255)
                                                                                                                                     ccccACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATc                 >  3:280568/1‑149 (MQ=255)
                                                                                                                                      cccACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATcc                >  4:539472/1‑149 (MQ=255)
                                                                                                                                      cccACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATcc                <  2:219430/149‑1 (MQ=255)
                                                                                                                                       ccaccaTTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCACCCGCCGCCGCCAGCGCCAGTGAAATCGGTTCGGTACATccc               <  2:369576/149‑1 (MQ=255)
                                                                                                                                       ccaccaTTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagcgcca                                      <  3:480859/126‑1 (MQ=255)
                                                                                                                                       ccaccaTTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagcgcca                                      >  4:480859/1‑126 (MQ=255)
                                                                                                                                        caccgTTCCCGGCA‑CGGGAATGGTGGCGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCa              >  4:675050/1‑149 (MQ=255)
                                                                                                                                        caccaTTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCa              >  2:68351/1‑149 (MQ=255)
                                                                                                                                        caccaTTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCa              >  3:594469/1‑149 (MQ=255)
                                                                                                                                         accgTTCCCGGCA‑CGGGAATGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAg             >  2:872542/1‑149 (MQ=255)
                                                                                                                                         accgTTCCCGGCA‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAg             >  4:357709/1‑149 (MQ=255)
                                                                                                                                         accaTTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGATCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAg             >  3:549921/1‑149 (MQ=255)
                                                                                                                                         accaTTCCCGTAC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAg             >  1:467551/1‑149 (MQ=255)
                                                                                                                                               cccGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCgg       <  3:51835/149‑1 (MQ=255)
                                                                                                                                                ccGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGt      >  4:669261/1‑149 (MQ=255)
                                                                                                                                                ccGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGt      <  2:165903/149‑1 (MQ=255)
                                                                                                                                                ccGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGt      >  3:340158/1‑149 (MQ=255)
                                                                                                                                                    gC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTa  <  1:63485/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTA  >  NC_000913/3256294‑3256590

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑