Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,337,626 T→C A159A (GCA→GCG adiA ← biodegradative arginine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,337,6260TC84.9% 687.0 / 115.1 298A159A (GCA→GCGadiAbiodegradative arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (26/19);  new base C (126/127);  total (152/146)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 304 total aligned reads displayed.

TCGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGTCATCGCGG  >  NC_000913/4337534‑4337714
                                                                                            |                                                                                        
  ggTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGcgccgc                                                                                        >  2:3954082/1‑93 (MQ=255)
  ggTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGcgccgc                                                                                        >  1:3605631/1‑93 (MQ=255)
   gTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCgccgcc                                                                                       >  1:5469717/1‑93 (MQ=255)
    tGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGccc                                                                                      >  1:1169182/1‑93 (MQ=255)
    tGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGccc                                                                                      >  2:3051804/1‑93 (MQ=255)
    tGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGccc                                                                                      <  2:1559880/93‑1 (MQ=255)
       ggAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAgg                                                                                   >  1:792512/1‑93 (MQ=255)
        gAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGa                                                                                  >  1:4123210/1‑93 (MQ=255)
        gAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGACGCTGCCCAGGa                                                                                  >  2:5419541/1‑93 (MQ=255)
          aCAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAat                                                                                >  2:7079402/1‑93 (MQ=255)
           cAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAata                                                                               <  2:5738941/93‑1 (MQ=255)
           cAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAata                                                                               <  1:729737/93‑1 (MQ=255)
           cAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAata                                                                               <  2:5278536/93‑1 (MQ=255)
            aGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAatat                                                                              <  1:6744911/93‑1 (MQ=255)
             gATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATAtt                                                                             <  2:2225154/93‑1 (MQ=255)
             gATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATAtt                                                                             >  1:4080065/1‑93 (MQ=255)
              aTTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTc                                                                            <  1:2983635/93‑1 (MQ=255)
               ttttCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCa                                                                           <  2:2066755/93‑1 (MQ=255)
               ttttCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCa                                                                           >  1:6668687/1‑93 (MQ=255)
                tttCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCAt                                                                          <  2:2619948/93‑1 (MQ=255)
                 ttCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATg                                                                         <  1:2099838/93‑1 (MQ=255)
                   cACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGa                                                                       <  2:5130383/93‑1 (MQ=255)
                   cACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGa                                                                       <  2:6792665/93‑1 (MQ=255)
                   cACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGa                                                                       <  2:431001/93‑1 (MQ=255)
                    aCCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGAt                                                                      <  2:281674/93‑1 (MQ=255)
                    aCCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGAt                                                                      <  2:3168560/93‑1 (MQ=255)
                    aCCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGAt                                                                      <  1:5537482/93‑1 (MQ=255)
                     ccATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATg                                                                     >  2:4466440/1‑93 (MQ=255)
                          gtagtCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACt                                                                <  1:3318547/93‑1 (MQ=255)
                          gtagtCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACt                                                                >  1:6245467/1‑93 (MQ=255)
                           tagtCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACta                                                               <  2:4620152/93‑1 (MQ=255)
                           tagtCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACta                                                               <  2:1677059/93‑1 (MQ=255)
                           tagtCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACta                                                               >  1:6436441/1‑93 (MQ=255)
                            agtCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACtat                                                              <  2:5879939/93‑1 (MQ=255)
                             gtCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACtata                                                             <  1:248524/93‑1 (MQ=255)
                             gtCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACtata                                                             >  2:2114902/1‑93 (MQ=255)
                             gtCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACtata                                                             <  2:2201087/93‑1 (MQ=255)
                               cATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATAtt                                                           >  2:6156107/1‑93 (MQ=255)
                               cATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATAtt                                                           >  2:2395109/1‑93 (MQ=255)
                                 tGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTtc                                                         >  1:6796663/1‑93 (MQ=255)
                                 tGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTtc                                                         >  2:6108473/1‑93 (MQ=255)
                                 tGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTtc                                                         >  2:1352580/1‑93 (MQ=255)
                                 tGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTtc                                                         <  2:4007147/93‑1 (MQ=255)
                                   gTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTtcat                                                       >  1:1633866/1‑93 (MQ=255)
                                   gTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTtcat                                                       >  2:4079363/1‑93 (MQ=255)
                                    tAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTtcatc                                                      <  1:1940125/93‑1 (MQ=255)
                                    tAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTtcatc                                                      <  1:883629/93‑1 (MQ=255)
                                    tAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTtcatc                                                      <  1:2682540/93‑1 (MQ=255)
                                       aaaCGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAgc                                                   >  2:5009547/1‑93 (MQ=255)
                                         aCGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAgcgc                                                 <  2:6200811/93‑1 (MQ=255)
                                         aCGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAgcgc                                                 >  1:4164543/1‑93 (MQ=255)
                                         aCGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAgcgc                                                 >  1:4880626/1‑93 (MQ=255)
                                           gTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAgcgcgc                                               >  2:3596514/1‑93 (MQ=255)
                                           gTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAgcgcgc                                               >  1:5726668/1‑93 (MQ=255)
                                            tCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCt                                              >  2:7348383/1‑93 (MQ=255)
                                              cggcggGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGa                                            >  2:5217396/1‑93 (MQ=255)
                                               ggcggGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGaa                                           <  2:3528641/93‑1 (MQ=255)
                                               ggcggGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGaa                                           <  2:6032984/93‑1 (MQ=255)
                                               ggcggGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGaa                                           <  2:5947932/93‑1 (MQ=255)
                                               ggcggGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGaa                                           >  1:5017974/1‑93 (MQ=255)
                                                   ggTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGt                                       <  1:2801701/93‑1 (MQ=255)
                                                   ggTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGt                                       >  1:3256276/1‑93 (MQ=255)
                                                   ggTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGt                                       <  1:5299342/93‑1 (MQ=255)
                                                   ggTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGt                                       <  1:1676952/93‑1 (MQ=255)
                                                    gtgtTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTg                                      <  2:4674491/93‑1 (MQ=255)
                                                    gtgtTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTg                                      <  1:2352510/93‑1 (MQ=255)
                                                      gtTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTggc                                    >  2:1877560/1‑93 (MQ=255)
                                                      gtTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTggc                                    <  2:6697531/93‑1 (MQ=255)
                                                       ttttGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTggcg                                   >  1:5595816/1‑93 (MQ=255)
                                                         ttGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTggcggc                                 <  2:3411392/93‑1 (MQ=255)
                                                           ggTAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCaa                               >  1:3154203/1‑93 (MQ=255)
                                                             tAAAACCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACa                             <  2:2044414/93‑1 (MQ=255)
                                                                aaCCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAgct                          >  2:4731970/1‑93 (MQ=255)
                                                                 aCCAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAgctg                         >  1:1583682/1‑93 (MQ=255)
                                                                  ccAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAgctgc                        >  2:6496013/1‑93 (MQ=255)
                                                                   cAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAgctgct                       >  2:5650566/1‑93 (MQ=255)
                                                                   cAACGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAgctgct                       >  1:3516312/1‑93 (MQ=255)
                                                                    aaCGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGctgctg                      >  1:3406287/1‑93 (MQ=255)
                                                                    aaCGCCGCCCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGctgctg                      >  1:5142676/1‑93 (MQ=255)
                                                                      cgccgcCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGc                    <  1:6558534/93‑1 (MQ=255)
                                                                       gccgccCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCg                   >  2:7517200/1‑93 (MQ=255)
                                                                        ccgccCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCgg                  <  1:4319171/93‑1 (MQ=255)
                                                                        ccgccCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCgg                  <  1:6965743/93‑1 (MQ=255)
                                                                        ccgccCTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCgg                  <  1:4796710/93‑1 (MQ=255)
                                                                             cTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCg             >  2:2547520/1‑93 (MQ=255)
                                                                             cTGGTGGCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCACCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCCGTGGCg             >  1:181223/1‑93 (MQ=255)
                                                                              tggtggCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCgg            >  1:2067417/1‑93 (MQ=255)
                                                                              tggtggCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCgg            >  1:3706735/1‑93 (MQ=255)
                                                                              tggtggCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCgg            >  1:4562512/1‑93 (MQ=255)
                                                                               ggtggCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCggg           >  1:709140/1‑93 (MQ=255)
                                                                               ggtggCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCggg           <  1:6038743/93‑1 (MQ=255)
                                                                                gtggCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGt          >  2:5817219/1‑93 (MQ=255)
                                                                                gtggCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGt          >  1:6368484/1‑93 (MQ=255)
                                                                                gtggCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGt          >  2:5507212/1‑93 (MQ=255)
                                                                                 tggCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGTc         >  2:6348806/1‑93 (MQ=255)
                                                                                 tggCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGTc         <  1:5852799/93‑1 (MQ=255)
                                                                                  ggCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGTCa        >  2:2585257/1‑93 (MQ=255)
                                                                                   gCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATCTTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGTCAt       >  2:2376647/1‑93 (MQ=255)
                                                                                   gCCTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGTCAt       >  2:3690003/1‑93 (MQ=255)
                                                                                     cTGGCGCCGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGTCATcg     >  1:1081286/1‑93 (MQ=255)
                                                                                            |                                                                                        
TCGGTGCGGAACAGATTTTCACCATAGTAGTCATGGTAGAAACGTCCGGCGGGTGTTTTGGTAAAACCAACGCCGCCCTGGTGGCCTGGCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGTCATCGCGG  >  NC_000913/4337534‑4337714

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