Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,724,099 T→G 8.7% intergenic (+17/‑29) pssA → / → yfiM phosphatidylserine synthase; CDP‑diacylglycerol‑serine O‑phosphatidyltransferase/putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,724,0990TG 8.7% 411.4 / 26.0 126intergenic (+17/‑29)pssA/yfiMphosphatidylserine synthase; CDP‑diacylglycerol‑serine O‑phosphatidyltransferase/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (60/55);  new base G (5/6);  total (65/61)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 126 total aligned reads displayed.

GCGAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCG  >  NC_000913/2723952‑2724245
                                                                                                                                                   |                                                                                                                                                  
gcgAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTa                                                                                                                                                   <  3:764766/149‑1 (MQ=255)
 cgAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTAc                                                                                                                                                  <  3:331423/149‑1 (MQ=255)
    gCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACggg                                                                                                                                               <  2:53222/149‑1 (MQ=255)
       tACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGgtt                                                                                                                                            <  4:285007/149‑1 (MQ=255)
       tACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgtt                                                                                                                                            <  3:708799/149‑1 (MQ=255)
       tACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgtt                                                                                                                                            <  4:513269/149‑1 (MQ=255)
         ccaccaTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgttgt                                                                                                                                          <  2:594025/149‑3 (MQ=255)
          caccaTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTATACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgttgtg                                                                                                                                         <  4:733088/149‑4 (MQ=255)
          caccaTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGgtttgt                                                                                                                                         >  2:428516/1‑149 (MQ=255)
          caccaTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGgtttgt                                                                                                                                         >  3:652792/1‑149 (MQ=255)
          caccaTCATTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGgtttgt                                                                                                                                         >  3:627272/1‑149 (MQ=255)
             caTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGtttt                                                                                                                                      >  2:446014/1‑149 (MQ=255)
             caTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGtttt                                                                                                                                      >  3:444246/1‑149 (MQ=255)
              aTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGttttt                                                                                                                                     <  1:374489/149‑1 (MQ=255)
              aTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGttttt                                                                                                                                     <  2:289362/149‑1 (MQ=255)
               tCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGtttttt                                                                                                                                    <  4:409122/149‑1 (MQ=255)
               tCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGtttttt                                                                                                                                    <  3:554341/149‑1 (MQ=255)
                cGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTg                                                                                                                                   <  4:316581/149‑1 (MQ=255)
                cGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTg                                                                                                                                   <  2:776765/149‑1 (MQ=255)
                  ttAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGa                                                                                                                                 <  3:874347/149‑1 (MQ=255)
                    aaGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAgg                                                                                                                               >  2:137048/1‑149 (MQ=255)
                     aGCACTATCGCGATCTGCAAAGTATTGGCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGTCCGATTAATTAGCCACATCCTGTAATCACAACCCCGTCCTGTACggg                                                                                                                                               <  2:238650/132‑1 (MQ=255)
                     aGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACggg                                                                                                                                               >  1:238650/1‑132 (MQ=255)
                      gCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGcc                                                                                                                             >  2:862320/1‑149 (MQ=255)
                       cACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCa                                                                                                                            >  3:70913/1‑149 (MQ=255)
                           aTCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGtt                                                                                                                        >  4:848069/1‑149 (MQ=255)
                             cgcgATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGtttt                                                                                                                      >  3:187759/1‑149 (MQ=255)
                              gcgATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTg                                                                                                                     >  1:442169/1‑149 (MQ=255)
                               cgATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGc                                                                                                                    <  1:426766/149‑1 (MQ=255)
                               cgATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGc                                                                                                                    <  3:85756/149‑1 (MQ=255)
                                        aaGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttt                                                                                                           >  4:530268/1‑149 (MQ=255)
                                        aaGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttt                                                                                                           <  3:10475/149‑1 (MQ=255)
                                          gTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttttt                                                                                                         <  4:488844/149‑1 (MQ=255)
                                          gTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttttt                                                                                                         <  1:423547/149‑1 (MQ=255)
                                          gTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttttt                                                                                                         <  3:607950/149‑1 (MQ=255)
                                          gTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttttt                                                                                                         <  1:572171/149‑1 (MQ=255)
                                             ttGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTc                                                                                                      <  3:11285/149‑1 (MQ=255)
                                                ccGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTc                                                                                                              >  2:244836/1‑138 (MQ=255)
                                                ccGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGt                                                                                                   >  2:190759/1‑149 (MQ=255)
                                                 cGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACggg                                                                                                                                               >  4:405395/1‑104 (MQ=255)
                                                 cGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACggg                                                                                                                                               <  3:405395/104‑1 (MQ=255)
                                                  taTTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACTTTTCGCGTATTCTTTTTGTCTGTTc                                                                                                 >  2:106099/2‑149 (MQ=255)
                                                  gATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAg                                                                                                                                >  2:677456/1‑118 (MQ=255)
                                                  gATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTc                                                                                                 >  2:553280/1‑149 (MQ=255)
                                                         cGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTg                                                                                          >  3:673271/1‑149 (MQ=255)
                                                          ggTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGt                                                                                         <  1:259934/149‑1 (MQ=255)
                                                           gTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGtt                                                                                        >  1:824558/1‑149 (MQ=255)
                                                              aaGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACt                                                                                     >  4:549057/1‑149 (MQ=255)
                                                                  ttCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCt                                                                                 >  4:619341/1‑149 (MQ=255)
                                                                   tCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTg                                                                                <  4:673271/149‑1 (MQ=255)
                                                                       aaaCTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATg                                                                            >  1:465361/1‑149 (MQ=255)
                                                                         aCTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTg                                                                                <  3:769233/143‑1 (MQ=255)
                                                                         aCTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTg                                                                                >  4:769233/1‑143 (MQ=255)
                                                                         aCTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCa                                                                          >  4:293985/1‑149 (MQ=255)
                                                                           tCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTAtt                                                                                           >  1:847831/1‑130 (MQ=255)
                                                                           tCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTAtt                                                                                           <  2:847831/130‑1 (MQ=255)
                                                                              tCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGGTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCat                                                                     >  1:51814/1‑149 (MQ=255)
                                                                                 gccgTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATg                                                                  >  4:494263/1‑149 (MQ=255)
                                                                                 gccgTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATg                                                                  >  2:557089/1‑149 (MQ=255)
                                                                                 gccgTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATg                                                                  <  1:862320/149‑1 (MQ=255)
                                                                                          gcCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACg                                                                                                                                                 <  4:168575/61‑1 (MQ=255)
                                                                                           ccGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATa                                                        >  3:563752/1‑149 (MQ=255)
                                                                                            cGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAg                                                       <  4:447644/149‑1 (MQ=255)
                                                                                            cGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAg                                                       >  1:218503/1‑149 (MQ=255)
                                                                                               aTCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTg                                                    <  3:293985/149‑1 (MQ=255)
                                                                                                 ccGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGa                                                  >  4:379463/1‑149 (MQ=255)
                                                                                                         aCCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAg                                          <  3:379463/149‑1 (MQ=255)
                                                                                                           cGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGa                                        <  1:774621/149‑1 (MQ=255)
                                                                                                             attaattaGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATa                                      >  3:309458/1‑149 (MQ=255)
                                                                                                             attaattaGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATa                                      >  3:165438/1‑149 (MQ=255)
                                                                                                                  ttaGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCt                                 >  1:697573/1‑149 (MQ=255)
                                                                                                                      ccGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAac                             >  3:635879/1‑149 (MQ=255)
                                                                                                                      ccGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAac                             >  2:671144/1‑149 (MQ=255)
                                                                                                                        gCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAacac                           >  2:717699/1‑149 (MQ=255)
                                                                                                                          aTCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTGTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGCGCGGGCAGGATAAAGCTTAACACtt                         >  4:145179/1‑149 (MQ=255)
                                                                                                                             cTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCATGTTTTGCGTATTGTTTTTGTCTGTTCGTTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTAt                      <  2:299684/149‑1 (MQ=255)
                                                                                                                              tGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATc                     >  1:255096/1‑149 (MQ=255)
                                                                                                                                tAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGc                   >  1:63305/1‑149 (MQ=255)
                                                                                                                                aaaaacaaaCCCCGTACAGGACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGc                   >  4:678716/8‑149 (MQ=255)
                                                                                                                                 aaTCACAATCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGCTGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGcc                  >  1:639136/1‑149 (MQ=255)
                                                                                                                                  aTCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCt                 >  1:464665/1‑149 (MQ=255)
                                                                                                                                   tCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTc                >  1:174538/1‑149 (MQ=255)
                                                                                                                                   tCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTc                <  4:187410/149‑1 (MQ=255)
                                                                                                                                   tCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTc                >  2:656504/1‑149 (MQ=255)
                                                                                                                                    acaaaCCCCGTACAGGACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCg               >  1:520078/4‑149 (MQ=255)
                                                                                                                                     acaACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCgg              >  1:166184/1‑149 (MQ=255)
                                                                                                                                       aaCCCCGTACAGGACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAg                                          >  4:840078/1‑119 (MQ=255)
                                                                                                                                       aaCCCCGTACAGGACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAg                                          <  3:840078/119‑1 (MQ=255)
                                                                                                                                       aaCCCCGTACAGGACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCg            >  1:336535/1‑149 (MQ=255)
                                                                                                                                         ccccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAg                                                                         <  4:5828/86‑1 (MQ=255)
                                                                                                                                         ccccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAg                                                                         >  3:5828/1‑86 (MQ=255)
                                                                                                                                         ccccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGAt          <  1:656504/149‑1 (MQ=255)
                                                                                                                                            cGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCt       >  3:438264/1‑149 (MQ=255)
                                                                                                                                            cGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCGAGGCGAtaaa       <  3:566918/149‑4 (MQ=255)
                                                                                                                                             gTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCtt      <  2:824558/149‑1 (MQ=255)
                                                                                                                                             gTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCtt      <  2:464665/149‑1 (MQ=255)
                                                                                                                                               ccTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCGAGGCGAtaaagtg    <  1:428516/149‑7 (MQ=255)
                                                                                                                                                cTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTcc   <  3:526956/149‑1 (MQ=255)
                                                                                                                                                cTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTcc   <  4:652792/149‑1 (MQ=255)
                                                                                                                                                 tGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTccg  >  3:801599/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTCCG  >  NC_000913/2723952‑2724245

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