Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,724,099 T→G 10.0% 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,0990TG10.0% 281.1 / 24.0 90intergenic (+17/‑29)pssA/yfiMphosphatidylserine synthase; CDP‑diacylglycerol‑serine O‑phosphatidyltransferase/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (40/41);  new base G (4/5);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

GCGAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTC  >  NC_000913/2723952‑2724243
                                                                                                                                                   |                                                                                                                                                
gcgAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTa                                                                                                                                                 >  4:313963/1‑149 (MQ=255)
   aGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACgg                                                                                                                                              <  4:327458/149‑1 (MQ=255)
   aGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACg                                                                                                                                               <  1:739292/148‑1 (MQ=255)
      aTACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGgt                                                                                                                                           >  3:535719/1‑149 (MQ=255)
       tACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgtt                                                                                                                                          <  3:275494/149‑1 (MQ=255)
        accaccATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGgttt                                                                                                                                         <  3:316955/149‑1 (MQ=255)
        accaccATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgttg                                                                                                                                         <  3:896028/149‑2 (MQ=255)
         ccaccaTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgttgt                                                                                                                                        <  3:730855/149‑3 (MQ=255)
         ccaccaTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgttg                                                                                                                                         <  1:796336/148‑2 (MQ=255)
           accaTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGgtttgtt                                                                                                                                      <  4:241813/149‑1 (MQ=255)
            ccaTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGgtttgttt                                                                                                                                     <  4:535719/149‑1 (MQ=255)
             caTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGtttt                                                                                                                                    >  4:310841/1‑149 (MQ=255)
             caTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGtttt                                                                                                                                    >  4:89516/1‑149 (MQ=255)
              aTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGtttt                                                                                                                                    >  2:863663/1‑148 (MQ=255)
               tCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGttttt                                                                                                                                   >  2:655658/1‑148 (MQ=255)
               tCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGttttt                                                                                                                                   >  2:540942/1‑148 (MQ=255)
               tCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGttttt                                                                                                                                   <  1:643280/148‑1 (MQ=255)
                 gTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTgg                                                                                                                                >  3:857801/1‑149 (MQ=255)
                  ttAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGa                                                                                                                               >  4:903759/1‑149 (MQ=255)
                   tAAGCACTATCGCGATCTGCAAGGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGa                                                                                                                               <  1:709745/148‑1 (MQ=255)
                   tAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGa                                                                                                                               <  1:540942/148‑1 (MQ=255)
                   tAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAg                                                                                                                              <  4:857801/149‑1 (MQ=255)
                        aCTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCa                                                                                                                          >  2:45111/1‑148 (MQ=255)
                           aTCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGtt                                                                                                                      >  3:725816/1‑149 (MQ=255)
                            tCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGttt                                                                                                                     >  4:608311/1‑149 (MQ=255)
                             cgcgATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGttt                                                                                                                     >  2:852031/1‑148 (MQ=255)
                               cgATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTg                                                                                                                   <  2:484109/148‑1 (MQ=255)
                                 aTCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCg                                                                                                                 <  1:863663/148‑1 (MQ=255)
                                   cTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTa                                                                                                               <  2:774626/148‑1 (MQ=255)
                                    tGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTAtt                                                                                                             >  3:615598/1‑149 (MQ=255)
                                     gCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTAtt                                                                                                             >  1:798002/1‑148 (MQ=255)
                                      cAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCt                                                                                                           <  4:791793/149‑1 (MQ=255)
                                          gTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttttt                                                                                                       <  4:725816/149‑1 (MQ=255)
                                            aTTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGt                                                                                                     <  3:310841/149‑1 (MQ=255)
                                              tGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTc                                                                                                    <  1:771199/148‑1 (MQ=255)
                                                ccGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgtt                                                                                                                                          <  3:741626/108‑1 (MQ=255)
                                                ccGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgtt                                                                                                                                          >  4:741626/1‑108 (MQ=255)
                                                 cGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGtt                                                                                                <  4:32767/149‑1 (MQ=255)
                                                 cGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACggg                                                                                                                                             >  2:859142/1‑104 (MQ=255)
                                                    ttATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGc                                                                                             <  4:841/149‑1 (MQ=255)
                                                      aTCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCt                                                                                            <  1:249457/148‑1 (MQ=255)
                                                        ccGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTAtt                                                                                         <  3:863646/149‑1 (MQ=255)
                                                        ccGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTAt                                                                                          <  1:602380/148‑1 (MQ=255)
                                                                   tCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTg                                                                              <  3:817321/149‑1 (MQ=255)
                                                                     gTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTgg                                                                             <  1:682244/148‑1 (MQ=255)
                                                                      tAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGa                                                                            >  2:753160/1‑148 (MQ=255)
                                                                      tAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGAt                                                                           >  4:666398/1‑149 (MQ=255)
                                                                      tAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGAt                                                                           >  3:254686/1‑149 (MQ=255)
                                                                              tCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCat                                                                   >  3:322390/1‑149 (MQ=255)
                                                                              tCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCat                                                                   <  4:615598/149‑1 (MQ=255)
                                                                                cgccgTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCtttt                                                                                                        >  4:519343/1‑110 (MQ=255)
                                                                                cgccgTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCtttt                                                                                                        <  3:519343/110‑1 (MQ=255)
                                                                                cgccgTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCata                                                                  <  1:45111/148‑1 (MQ=255)
                                                                                      ttGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTa                                                            <  1:485059/148‑1 (MQ=255)
                                                                                             gTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGc                                                    <  3:903759/149‑1 (MQ=255)
                                                                                                      tCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCggg                                           <  4:333816/149‑1 (MQ=255)
                                                                                                           cGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAgg                                       <  1:892244/148‑1 (MQ=255)
                                                                                                            gATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGa                                      >  1:897934/1‑148 (MQ=255)
                                                                                                                 attaGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGc                                >  3:422758/1‑149 (MQ=255)
                                                                                                                   taGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCt                               >  2:506472/1‑148 (MQ=255)
                                                                                                                        gCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGa                                      <  3:805580/136‑1 (MQ=255)
                                                                                                                        gCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGa                                      >  4:805580/1‑136 (MQ=255)
                                                                                                                            ccTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTa                     >  4:816265/1‑149 (MQ=255)
                                                                                                                              tGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATc                   <  3:89516/149‑1 (MQ=255)
                                                                                                                              tGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATc                   >  4:311763/1‑149 (MQ=255)
                                                                                                                              tGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTGTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTAt                    >  2:477789/1‑148 (MQ=255)
                                                                                                                                tAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGc                 <  4:422758/149‑1 (MQ=255)
                                                                                                                                  aTCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCt               >  4:837329/1‑149 (MQ=255)
                                                                                                                                   tCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGtt                                                                                      <  4:435070/77‑1 (MQ=255)
                                                                                                                                   tCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGtt                                                                                      >  3:435070/1‑77 (MQ=255)
                                                                                                                                    cacaACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTc              >  1:860006/1‑148 (MQ=255)
                                                                                                                                    cacaACCCCGTCATGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCg             >  4:905197/1‑149 (MQ=255)
                                                                                                                                    acaaaCCCCGTACAGGACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCg             >  3:83444/4‑149 (MQ=255)
                                                                                                                                     acaACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGACTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCgg            <  3:220488/149‑1 (MQ=255)
                                                                                                                                       aaCCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCa                                                                    <  4:427631/91‑1 (MQ=255)
                                                                                                                                       aaCCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCa                                                                    >  3:427631/1‑91 (MQ=255)
                                                                                                                                        aCCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCg          >  2:227687/1‑148 (MQ=255)
                                                                                                                                        aCCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGa         >  3:663080/1‑149 (MQ=255)
                                                                                                                                        aCCCCGTACAGGACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCg          >  2:810759/1‑148 (MQ=255)
                                                                                                                                         ccccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGa         <  2:798002/148‑1 (MQ=255)
                                                                                                                                         ccccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGAt        >  4:814443/1‑149 (MQ=255)
                                                                                                                                         ccccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCACGATAAAGCTCAACACTTTATCGCCTCGGCGa         >  1:363291/1‑148 (MQ=255)
                                                                                                                                          cccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATg       <  3:311763/149‑1 (MQ=255)
                                                                                                                                           ccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCAGGTTTTGCGTATTCTTTTTGTTTGTTCGCTATAGTTACTTTCTGGATGTAGCCATATGGCTAACGATAGTTCGAGCGGGCAGGATAAGGCTCAACACTTTATCGCCGAGGCGAtaa      <  3:546162/149‑3 (MQ=255)
                                                                                                                                           ccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATg       >  1:188849/1‑148 (MQ=255)
                                                                                                                                           ccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTGTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATg       >  2:206360/1‑148 (MQ=255)
                                                                                                                                             gTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAg                                               <  2:809458/106‑1 (MQ=255)
                                                                                                                                             gTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAg                                               >  1:809458/1‑106 (MQ=255)
                                                                                                                                             gTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCt     >  2:659394/1‑148 (MQ=255)
                                                                                                                                               ccTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTc  <  4:768967/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCGAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTTC  >  NC_000913/2723952‑2724243

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: