Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,256,450 C→T 7.9% V68I (GTT→ATT)  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,4500CT 7.9% 301.8 / 10.1 99V68I (GTT→ATT) yhaMputative L‑serine dehydratase alpha chain
Reads supporting (aligned to +/- strand):  ref base C (47/44);  new base T (3/5);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTACTTCCTC  >  NC_000913/3256305‑3256597
                                                                                                                                                  |                                                                                                                                                   
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:667887/149‑1 (MQ=255)
                          ggCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGAt                                                                                                                        <  2:442428/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                                                                                                                    >  2:722403/1‑149 (MQ=255)
                              tCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTgg                                                                                                                    >  2:63485/1‑149 (MQ=255)
                              tCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTgg                                                                                                                    >  1:163656/1‑149 (MQ=255)
                              tCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTgg                                                                                                                    <  2:108228/149‑1 (MQ=255)
                                ggCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGACg                                                                                                                  >  3:702757/1‑149 (MQ=255)
                                 gCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGa                                                                                                                 >  2:811100/1‑149 (MQ=255)
                                 gCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGa                                                                                                                 <  2:676645/149‑1 (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                                                                                            <  4:457982/149‑1 (MQ=255)
                                                      tCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTc                                                                                            <  3:843163/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                                                                                                                                                 <  1:429887/90‑1 (MQ=255)
                                                            aGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCG‑GCACGGGAATGGTg                                                                                                                                                 >  2:429887/1‑90 (MQ=255)
                                                                  tCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTc                                                                                <  1:886776/149‑1 (MQ=255)
                                                                     aGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGGTTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAg                                                                             >  2:607253/1‑149 (MQ=255)
                                                                            cGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGcc                                                                      <  3:858470/149‑1 (MQ=255)
                                                                            cGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGcc                                                                      <  2:374502/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                                                                           <  2:536534/136‑1 (MQ=255)
                                                                                    tttCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGtt                                                                           >  1:536534/1‑136 (MQ=255)
                                                                                    tttCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAAc                                                              >  4:141919/1‑149 (MQ=255)
                                                                                    tttCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAAc                                                              >  4:51835/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                                                                                                   >  3:813314/1‑96 (MQ=255)
                                                                                                    ccAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTa                                                                                                   <  4:813314/96‑1 (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                                    <  1:584809/149‑1 (MQ=255)
                                                                                                              cgcAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCg                                    >  1:630752/1‑149 (MQ=255)
                                                                                                              cgcAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCg                                    >  2:506427/1‑149 (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                     >  3:594469/1‑149 (MQ=255)
                                                                                                                             caccaTTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCa                     >  2:68351/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)
                                                                                                                                          c‑cGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTAc        <  3:85276/149‑1 (MQ=255)
                                                                                                                                             gggAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTACtt      <  1:218046/149‑1 (MQ=255)
                                                                                                                                               gAACAGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTAAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTACTtcc    <  1:68351/149‑1 (MQ=255)
                                                                                                                                                 aCGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTACTtcctc  <  2:448138/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGC‑CGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTTACTTCCTC  >  NC_000913/3256305‑3256597

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