Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,551,460 T→C A72A (GCT→GCC hemF → coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,551,4600TC88.5% 523.1 / 57.7 209A72A (GCT→GCChemFcoproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base T (12/12);  new base C (96/89);  total (108/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

Only 100 of 212 total aligned reads displayed.

CGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTT  >  NC_000913/2551370‑2551550
                                                                                          |                                                                                          
cggcggGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTc                                                                                          <  1:2397970/93‑1 (MQ=255)
 ggcgggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTcc                                                                                         >  2:3268431/1‑93 (MQ=255)
     ggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCa                                                                                     >  1:6900096/1‑93 (MQ=255)
     ggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCa                                                                                     >  1:97449/1‑93 (MQ=255)
      gcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCAc                                                                                    <  1:2646144/93‑1 (MQ=255)
      gcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCAc                                                                                    <  1:4049495/93‑1 (MQ=255)
        gTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCg                                                                                  >  2:1508775/1‑93 (MQ=255)
         tAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGc                                                                                 >  2:7339979/1‑93 (MQ=255)
          aGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCt                                                                                >  2:3950365/1‑93 (MQ=255)
            tCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCa                                                                              >  1:2545105/1‑93 (MQ=255)
              gggTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATc                                                                            <  1:7158250/93‑1 (MQ=255)
              gggTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATc                                                                            >  2:5264805/1‑93 (MQ=255)
               ggTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCg                                                                           <  2:2137711/93‑1 (MQ=255)
                 tgtTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGcc                                                                         >  2:4078833/1‑93 (MQ=255)
                 tgtTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGcc                                                                         <  2:3053546/93‑1 (MQ=255)
                 tgtTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGcc                                                                         >  2:2938043/1‑93 (MQ=255)
                 tgtTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGcc                                                                         >  1:3523228/1‑93 (MQ=255)
                    tGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCgg                                                                      >  1:4727875/1‑93 (MQ=255)
                       gTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAAc                                                                   >  1:6029805/1‑93 (MQ=255)
                         aaTGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACtt                                                                 >  1:1728273/1‑93 (MQ=255)
                           tggtggTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGc                                                               >  2:840749/1‑93 (MQ=255)
                             gtggtgTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCg                                                             <  2:137537/93‑1 (MQ=255)
                              tggtgTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCgg                                                            <  2:5209518/93‑1 (MQ=255)
                              tggtgTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCgg                                                            >  2:1575957/1‑93 (MQ=255)
                                gtgtTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgc                                                          <  2:5022547/93‑1 (MQ=255)
                                gtgtTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgc                                                          <  2:4384295/93‑1 (MQ=255)
                                  gtTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgcgc                                                        <  1:6100716/93‑1 (MQ=255)
                                  gtTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgcgc                                                        <  2:2068135/93‑1 (MQ=255)
                                   ttttCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCa                                                       >  1:5536631/1‑93 (MQ=255)
                                      tCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGtt                                                    <  1:3991550/93‑1 (MQ=255)
                                      tCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGtt                                                    <  1:6050329/93‑1 (MQ=255)
                                       cGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGttt                                                   <  1:1623942/93‑1 (MQ=255)
                                        gAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTc                                                  <  1:2696119/93‑1 (MQ=255)
                                        gAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTc                                                  <  2:4658008/93‑1 (MQ=255)
                                         aaCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCg                                                 >  2:7246883/1‑93 (MQ=255)
                                         aaCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCg                                                 >  1:1614556/1‑93 (MQ=255)
                                          aCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGa                                                >  1:5397641/1‑93 (MQ=255)
                                          aCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCACCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGa                                                <  2:2118838/93‑1 (MQ=255)
                                           caggcaggCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAg                                               <  2:4708985/93‑1 (MQ=255)
                                            aggcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAgg                                              <  2:4909677/93‑1 (MQ=255)
                                             ggcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGc                                             >  2:4742702/1‑93 (MQ=255)
                                             ggcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGc                                             >  1:2171770/1‑93 (MQ=255)
                                             ggcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCAGGAACTTGCCGGGCGCAGTTTCGAGGc                                             <  1:7256790/93‑1 (MQ=255)
                                              gcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGACCGGAACTTGCCGGGCGCAGTTTCGAGGCg                                            >  1:2743731/1‑93 (MQ=255)
                                              gcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCg                                            >  2:1947675/1‑93 (MQ=255)
                                              gcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCg                                            >  1:4235724/1‑93 (MQ=255)
                                               caggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGa                                           <  2:4778553/93‑1 (MQ=255)
                                               caggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGa                                           <  2:7498900/93‑1 (MQ=255)
                                                aggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGAt                                          <  1:1899320/93‑1 (MQ=255)
                                                 ggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATg                                         >  1:5575965/1‑93 (MQ=255)
                                                 ggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATg                                         <  2:1465593/93‑1 (MQ=255)
                                                  gcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATgg                                        <  2:1520433/93‑1 (MQ=255)
                                                  gcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATgg                                        <  2:6541734/93‑1 (MQ=255)
                                                  gcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATgg                                        >  1:1346497/1‑93 (MQ=255)
                                                  gcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATgg                                        <  1:4558067/93‑1 (MQ=255)
                                                   cGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATggg                                       <  1:2464792/93‑1 (MQ=255)
                                                   cGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATggg                                       <  2:3766293/93‑1 (MQ=255)
                                                   cGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATggg                                       <  1:4953023/93‑1 (MQ=255)
                                                    gTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGc                                      <  2:2528798/93‑1 (MQ=255)
                                                     tCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCg                                     >  2:5851938/1‑93 (MQ=255)
                                                       aaCTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGtt                                   >  2:1355245/1‑93 (MQ=255)
                                                        aCTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGttt                                  >  2:1786209/1‑93 (MQ=255)
                                                         cTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTc                                 <  2:2056127/93‑1 (MQ=255)
                                                          ttttCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCa                                >  2:2351237/1‑93 (MQ=255)
                                                          ttttCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCa                                <  1:6994491/93‑1 (MQ=255)
                                                           tttCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCAc                               <  1:6584569/93‑1 (MQ=255)
                                                           tttCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCAc                               >  2:6374377/1‑93 (MQ=255)
                                                           tttCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCAc                               <  1:6311963/93‑1 (MQ=255)
                                                           tttCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCAc                               >  1:3178410/1‑93 (MQ=255)
                                                              cGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTgg                            <  1:6201346/93‑1 (MQ=255)
                                                              cGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTgg                            <  1:7087847/93‑1 (MQ=255)
                                                                cATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTa                          >  2:1218536/1‑93 (MQ=255)
                                                                 aTGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAg                         >  1:7319976/1‑93 (MQ=255)
                                                                  tGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGt                        <  2:5792396/93‑1 (MQ=255)
                                                                   gTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTg                       >  2:3192141/1‑93 (MQ=255)
                                                                    tCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGc                      >  1:4469223/1‑93 (MQ=255)
                                                                    tCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGc                      >  1:5721664/1‑93 (MQ=255)
                                                                       aCGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATc                   >  1:850259/1‑93 (MQ=255)
                                                                        cGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATcc                  <  1:6670207/93‑1 (MQ=255)
                                                                         ggTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCg                 >  1:2897725/1‑93 (MQ=255)
                                                                         ggTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCg                 >  2:4229832/1‑93 (MQ=255)
                                                                          gTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGc                <  2:4406364/93‑1 (MQ=255)
                                                                          gTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGc                >  2:6935568/1‑93 (MQ=255)
                                                                             aGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATa             <  1:1505979/93‑1 (MQ=255)
                                                                              ggCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATaa            >  1:4003016/1‑93 (MQ=255)
                                                                                cGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAAcc          >  1:7522605/1‑93 (MQ=255)
                                                                                 gATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAAccc         >  1:3761401/1‑93 (MQ=255)
                                                                                 gATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAAccc         >  1:5289080/1‑93 (MQ=255)
                                                                                 gATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAAccc         <  1:3481218/93‑1 (MQ=255)
                                                                                  aTGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCg        >  2:5744111/1‑93 (MQ=255)
                                                                                  aTGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCg        >  1:3686760/1‑93 (MQ=255)
                                                                                  aTGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCg        >  1:4826350/1‑93 (MQ=255)
                                                                                  aTGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCg        >  2:4482884/1‑93 (MQ=255)
                                                                                  aTGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCg        >  2:3289664/1‑93 (MQ=255)
                                                                                   tgcctgccTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGt       >  1:3974032/1‑93 (MQ=255)
                                                                                     ccTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTAt     <  1:2806277/93‑1 (MQ=255)
                                                                                      cTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATg    <  2:5056859/93‑1 (MQ=255)
                                                                                      cTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATg    <  2:6174835/93‑1 (MQ=255)
                                                                                       tGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGt   >  2:2730593/1‑93 (MQ=255)
                                                                                        gCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGtt  <  2:5453079/93‑1 (MQ=255)
                                                                                          |                                                                                          
CGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTT  >  NC_000913/2551370‑2551550

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