Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,632,259 T→C A207A (GCA→GCG rob ← right origin‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,632,2590TC86.3% 682.9 / 99.0 285A207A (GCA→GCGrobright origin‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (18/21);  new base C (135/111);  total (153/132)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

Only 100 of 286 total aligned reads displayed.

CTGCACGCCGGTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCAC  >  NC_000913/4632169‑4632351
                                                                                          |                                                                                            
cTGCACGCCGGTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGc                                                                                            >  2:5417435/1‑93 (MQ=255)
cTGCACGCCGGTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGc                                                                                            <  1:2792351/93‑1 (MQ=255)
   cACGCCGGTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTg                                                                                         >  1:4178349/1‑93 (MQ=255)
   cACGCCGGTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTg                                                                                         >  2:7462569/1‑93 (MQ=255)
   cACGCCGGTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTg                                                                                         <  2:2978322/93‑1 (MQ=255)
     cGCCGGTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGAt                                                                                       <  2:7339433/93‑1 (MQ=255)
      gCCGGTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATc                                                                                      >  2:7007302/1‑93 (MQ=255)
      gCCGGTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATc                                                                                      >  1:4952191/1‑93 (MQ=255)
         ggTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTg                                                                                   <  1:5956963/93‑1 (MQ=255)
         ggTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTg                                                                                   >  2:1715744/1‑93 (MQ=255)
         ggTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTg                                                                                   >  1:337152/1‑93 (MQ=255)
         ggTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTg                                                                                   >  2:5758149/1‑93 (MQ=255)
            tCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGc                                                                                >  2:1605180/1‑93 (MQ=255)
              ccAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTa                                                                              >  1:5290342/1‑93 (MQ=255)
               cAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTaa                                                                             >  1:2966081/1‑93 (MQ=255)
               cAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTaa                                                                             >  2:211384/1‑93 (MQ=255)
               cAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTaa                                                                             >  2:1210825/1‑93 (MQ=255)
                 gACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAAcg                                                                           >  2:3816948/1‑93 (MQ=255)
                 gACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAAcg                                                                           >  2:4138313/1‑93 (MQ=255)
                 gACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAAcg                                                                           >  2:5656117/1‑93 (MQ=255)
                  aCCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAAcgc                                                                          <  1:1943050/93‑1 (MQ=255)
                     ttCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCggt                                                                       <  2:695608/93‑1 (MQ=255)
                     ttCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCggt                                                                       <  1:7396547/93‑1 (MQ=255)
                     ttCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCggt                                                                       <  1:7093860/93‑1 (MQ=255)
                      tCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCggtg                                                                      <  1:3106186/93‑1 (MQ=255)
                        aTAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCggtggt                                                                    >  1:4600502/1‑93 (MQ=255)
                        aTAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCggtggt                                                                    >  1:718161/1‑93 (MQ=255)
                         tAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGta                                                                   >  1:5759775/1‑93 (MQ=255)
                         tAGGTAAACATCACATATTCGCCACCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGta                                                                   >  2:243455/1‑93 (MQ=255)
                          aGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGtat                                                                  >  2:1110768/1‑93 (MQ=255)
                          aGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGtat                                                                  >  2:2275157/1‑93 (MQ=255)
                          aGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGtat                                                                  >  1:3727101/1‑93 (MQ=255)
                           ggTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGtata                                                                 <  1:6753675/93‑1 (MQ=255)
                            gTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAg                                                                >  1:6378954/1‑93 (MQ=255)
                             tAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGa                                                               <  1:244745/93‑1 (MQ=255)
                              aaaCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGaa                                                              >  1:5804214/1‑93 (MQ=255)
                              aaaCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGaa                                                              >  2:5869895/1‑93 (MQ=255)
                               aaCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAAt                                                             >  2:3605108/1‑93 (MQ=255)
                                aCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATa                                                            <  1:5480007/93‑1 (MQ=255)
                                aCATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATa                                                            <  2:5102961/93‑1 (MQ=255)
                                 cATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATAc                                                           >  1:291636/1‑93 (MQ=255)
                                 cATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATAc                                                           >  1:2715911/1‑93 (MQ=255)
                                 cATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATAc                                                           <  1:2486318/93‑1 (MQ=255)
                                  aTCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTATCCTGGGCTAACGCGGTGGTATAGAATACct                                                         <  2:1472098/93‑1 (MQ=255)
                                  aTCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATAcc                                                          <  1:983962/93‑1 (MQ=255)
                                  aTCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATAcc                                                          <  1:6905851/93‑1 (MQ=255)
                                  aTCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATACTGGGCTAACGCGGTGGTATAGAATAcc                                                          <  2:5893722/93‑1 (MQ=255)
                                   tCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACct                                                         >  2:4812740/1‑93 (MQ=255)
                                    cacaTATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACctc                                                        <  1:4269062/93‑1 (MQ=255)
                                    cacaTATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACctc                                                        <  2:4908041/93‑1 (MQ=255)
                                     acaTATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACctct                                                       <  1:860735/93‑1 (MQ=255)
                                      caTATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCtt                                                      >  2:6592868/1‑93 (MQ=255)
                                        tatTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGc                                                    >  1:7367846/1‑93 (MQ=255)
                                        tatTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGc                                                    <  2:3948319/93‑1 (MQ=255)
                                        tatTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGc                                                    >  1:1553709/1‑93 (MQ=255)
                                         atTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCt                                                   <  2:1968478/93‑1 (MQ=255)
                                          ttCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCtc                                                  >  1:2579251/1‑93 (MQ=255)
                                             gccgccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCtcgtc                                               >  2:5362266/1‑93 (MQ=255)
                                             gccgccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCtcgtc                                               >  2:3986310/1‑93 (MQ=255)
                                             gccgccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCtcgtc                                               <  1:1098086/93‑1 (MQ=255)
                                              ccgccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCtcgtcg                                              <  2:7263227/93‑1 (MQ=255)
                                               cgccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                             >  2:3356590/1‑93 (MQ=255)
                                                gccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCgtcgtc                                            <  1:624622/93‑1 (MQ=255)
                                                gccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCgtcgtc                                            >  2:2827278/1‑93 (MQ=255)
                                                 cccTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCt                                           <  1:157498/93‑1 (MQ=255)
                                                  ccTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCtt                                          >  2:790982/1‑93 (MQ=255)
                                                  ccTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCtt                                          >  2:6782841/1‑93 (MQ=255)
                                                  ccTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCtt                                          >  1:4589298/1‑93 (MQ=255)
                                                   cTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCttt                                         >  1:3596020/1‑93 (MQ=255)
                                                   cTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCttt                                         >  1:6996255/1‑93 (MQ=255)
                                                   cTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCttt                                         <  1:6098855/93‑1 (MQ=255)
                                                       agcaTCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATcc                                     <  2:3207413/93‑1 (MQ=255)
                                                        gcaTCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCt                                    >  1:4395613/1‑93 (MQ=255)
                                                           tCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGAc                                 <  1:1904618/93‑1 (MQ=255)
                                                           tCACCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGAc                                 >  2:6136008/1‑93 (MQ=255)
                                                             aCCGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTc                               <  2:4323796/93‑1 (MQ=255)
                                                              ccGGATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCg                              >  1:4347216/1‑93 (MQ=255)
                                                                ggATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGa                            <  1:1767008/93‑1 (MQ=255)
                                                                ggATGCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGa                            <  2:5523802/93‑1 (MQ=255)
                                                                    gCCCCGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGAcgcg                        >  1:6821634/1‑93 (MQ=255)
                                                                     ccccGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGt                       >  2:1768451/1‑93 (MQ=255)
                                                                      cccGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGtt                      <  1:2317631/93‑1 (MQ=255)
                                                                      cccGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGtt                      >  1:5142105/1‑93 (MQ=255)
                                                                       ccGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGttt                     >  2:466978/1‑93 (MQ=255)
                                                                       ccGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGttt                     <  2:1801767/93‑1 (MQ=255)
                                                                        cGTCAGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTttc                    <  1:1235751/93‑1 (MQ=255)
                                                                          tCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTttcat                  <  2:269724/93‑1 (MQ=255)
                                                                            aGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTttcattc                >  1:1102945/1‑93 (MQ=255)
                                                                            aGTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTttcattc                >  1:3612675/1‑93 (MQ=255)
                                                                             gTACATAGCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTttcattca               <  2:2226942/93‑1 (MQ=255)
                                                                                    gCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTa        <  1:5588780/93‑1 (MQ=255)
                                                                                    gCCATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTa        <  1:2514714/93‑1 (MQ=255)
                                                                                     ccATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTag       >  1:2637613/1‑93 (MQ=255)
                                                                                     ccATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTag       <  2:6339479/93‑1 (MQ=255)
                                                                                     ccATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTag       <  2:4251484/93‑1 (MQ=255)
                                                                                     ccATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTag       >  2:7153946/1‑93 (MQ=255)
                                                                                      cATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTaga      >  1:4132567/1‑93 (MQ=255)
                                                                                      cATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTaga      >  1:6582091/1‑93 (MQ=255)
                                                                                      cATCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTaga      >  2:847340/1‑93 (MQ=255)
                                                                                       aTCCGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTagag     >  2:5248209/1‑93 (MQ=255)
                                                                                          |                                                                                            
CTGCACGCCGGTTCCCAGACCTTCATAGGTAAACATCACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCAC  >  NC_000913/4632169‑4632351

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