Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,572,928 C→G 5.3% A759P (GCA→CCA)  pqqL ← putative periplasmic M16 family zinc metalloendopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,572,9280CG 5.3% 337.8 / 7.5 94A759P (GCA→CCA) pqqLputative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base C (37/52);  new base G (2/3);  total (39/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTA  >  NC_000913/1572780‑1573076
                                                                                                                                                    |                                                                                                                                                    
gaagaaACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGc                                                                                                                                                      >  2:280848/1‑149 (MQ=255)
 aagaaACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGcc                                                                                                                                                     >  4:143553/1‑149 (MQ=255)
     aaCGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAg                                                                                                                                                 <  3:312513/149‑1 (MQ=255)
       cGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGt                                                                                                                                               <  2:708622/149‑1 (MQ=255)
         cTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTc                                                                                                                                             <  4:880101/149‑1 (MQ=255)
         cTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTc                                                                                                                                             <  2:15851/149‑1 (MQ=255)
           gTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAt                                                                                                                                           <  4:139427/149‑1 (MQ=255)
            tATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAACCCTGTGCCACAGGTTCAtt                                                                                                                                          <  4:635911/149‑1 (MQ=255)
            tATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtt                                                                                                                                          <  3:428783/149‑1 (MQ=255)
            tATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtt                                                                                                                                          >  1:869116/1‑149 (MQ=255)
            tATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtt                                                                                                                                          <  1:707575/149‑1 (MQ=255)
            tATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGATAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtt                                                                                                                                          >  2:732055/1‑149 (MQ=255)
             aTGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcaca                                                                                                                                         <  3:529261/149‑6 (MQ=255)
             aTGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAttt                                                                                                                                         <  1:129957/149‑1 (MQ=255)
              tGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtttt                                                                                                                                        <  1:819787/149‑1 (MQ=255)
                cTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGt                                                                                                                                      <  3:708054/149‑1 (MQ=255)
                cTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGt                                                                                                                                      >  1:56435/1‑149 (MQ=255)
                 tCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGtt                                                                                                                                     >  3:78089/1‑149 (MQ=255)
                 tCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGtt                                                                                                                                     >  2:745709/1‑149 (MQ=255)
                    aGATGCCTGTTCACGAATATTAACACGTAGATCTTGTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGTGTCGGCAGATTAACAAGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCtt                                                                                                                                  <  2:176306/149‑1 (MQ=255)
                     gATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCttt                                                                                                                                 <  1:274761/149‑1 (MQ=255)
                       tGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTa                                                                                                                               <  3:663584/149‑1 (MQ=255)
                        gCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTAc                                                                                                                              <  4:761127/149‑1 (MQ=255)
                         ccTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACa                                                                                                                             <  2:254377/149‑1 (MQ=255)
                         ccTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACa                                                                                                                             <  1:218967/149‑1 (MQ=255)
                          cTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAg                                                                                                                            >  3:445097/1‑149 (MQ=255)
                           tGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGt                                                                                                                           >  1:333767/1‑149 (MQ=255)
                                aCGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCg                                                                                                                      >  4:39988/1‑149 (MQ=255)
                                 cGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGa                                                                                                                     <  3:818012/149‑1 (MQ=255)
                                 cGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGa                                                                                                                     >  1:134457/1‑149 (MQ=255)
                                 cGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGa                                                                                                                     >  1:311263/1‑149 (MQ=255)
                                  gAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGAt                                                                                                                    >  1:681102/1‑149 (MQ=255)
                                  gAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGAt                                                                                                                    >  1:207775/1‑149 (MQ=255)
                                   aaTATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATg                                                                                                                   <  2:681102/149‑1 (MQ=255)
                                    atatTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGc                                                                                                                  <  4:378363/149‑1 (MQ=255)
                                         aaCACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGt                                                                                                             >  4:129368/1‑149 (MQ=255)
                                           cacGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTcc                                                                                                           <  4:202401/149‑1 (MQ=255)
                                                  aTCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcgcgc                                                                                                    >  2:378973/1‑149 (MQ=255)
                                                         gCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTaa                                                                                             >  3:101775/1‑149 (MQ=255)
                                                          ccAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAt                                                                                            >  3:445302/1‑149 (MQ=255)
                                                          ccAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAt                                                                                            <  2:221992/149‑1 (MQ=255)
                                                           cAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATg                                                                                           <  1:387718/149‑1 (MQ=255)
                                                           cAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATg                                                                                           >  4:777488/1‑149 (MQ=255)
                                                           cAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATg                                                                                           >  1:241755/1‑149 (MQ=255)
                                                                  aCGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTAcc                                                                                    <  2:352424/149‑1 (MQ=255)
                                                                   cGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCt                                                                                   <  1:137792/149‑1 (MQ=255)
                                                                     ttAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGc                                                                                 <  4:445302/149‑1 (MQ=255)
                                                                       aaaaGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCgg                                                                               <  2:698116/149‑1 (MQ=255)
                                                                           gCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTaa                                                                           >  1:41761/1‑149 (MQ=255)
                                                                               cGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGc                                                                       <  4:125110/149‑1 (MQ=255)
                                                                                 aGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGa                                                                     <  4:445097/149‑1 (MQ=255)
                                                                                    gcCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATc                                                                  <  3:164696/149‑1 (MQ=255)
                                                                                       aTACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCaga                                                               <  2:333767/149‑1 (MQ=255)
                                                                                         aCGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAgt                                                             <  2:56435/149‑1 (MQ=255)
                                                                                               cGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttga                                                       <  2:102318/149‑1 (MQ=255)
                                                                                                ggCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgat                                                      <  2:669411/149‑1 (MQ=255)
                                                                                                 gCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTaa                                                                                             <  3:408214/109‑1 (MQ=255)
                                                                                                 gCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTaa                                                                                             >  4:408214/1‑109 (MQ=255)
                                                                                                 gCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgatt                                                     >  4:756301/1‑149 (MQ=255)
                                                                                                     aTTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATc                                                 <  3:129368/149‑1 (MQ=255)
                                                                                                      ttAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATcc                                                >  3:614839/1‑149 (MQ=255)
                                                                                                      ttAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATcc                                                >  4:73488/1‑149 (MQ=255)
                                                                                                       tAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGACACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCt                                               >  1:67361/1‑149 (MQ=255)
                                                                                                         aCAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTaa                                             <  4:428909/149‑1 (MQ=255)
                                                                                                         aCAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTaa                                             >  4:582530/1‑149 (MQ=255)
                                                                                                         aCAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTaa                                             <  2:241755/149‑1 (MQ=255)
                                                                                                          cAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATc                                                 <  4:169876/144‑1 (MQ=255)
                                                                                                          cAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATc                                                 >  3:169876/1‑144 (MQ=255)
                                                                                                           aGGTGTCCGGGAATCATAACGCTTACACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGAAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGt                                           >  4:168083/1‑149 (MQ=255)
                                                                                                              tgtCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAc                                        >  4:355768/1‑149 (MQ=255)
                                                                                                                tCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAcgc                                      >  4:675259/1‑149 (MQ=255)
                                                                                                                   gggAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTa                                   <  3:769877/149‑1 (MQ=255)
                                                                                                                    ggAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcgc                                                                                                      >  4:128522/1‑81 (MQ=255)
                                                                                                                    ggAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcgc                                                                                                      <  3:128522/81‑1 (MQ=255)
                                                                                                                    ggAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTaa                                  <  4:412906/149‑1 (MQ=255)
                                                                                                                        tCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTaa                              <  1:280848/149‑1 (MQ=255)
                                                                                                                         cATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAAc                             <  1:745709/149‑1 (MQ=255)
                                                                                                                           tAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAc                                        >  1:527776/1‑136 (MQ=255)
                                                                                                                           tAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAc                                        <  2:527776/136‑1 (MQ=255)
                                                                                                                           tAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGc                           <  2:381161/149‑1 (MQ=255)
                                                                                                                           tAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGc                           <  3:538212/149‑1 (MQ=255)
                                                                                                                           tAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGc                           >  3:161930/1‑149 (MQ=255)
                                                                                                                                  tCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGt                    >  2:435352/1‑149 (MQ=255)
                                                                                                                                   ccACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGtt                   <  3:620198/149‑1 (MQ=255)
                                                                                                                                    cACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACATTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGttt                  <  1:732055/149‑1 (MQ=255)
                                                                                                                                    cACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGTGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGttt                  >  3:761365/1‑149 (MQ=255)
                                                                                                                                    cACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGttt                  >  3:393062/1‑149 (MQ=255)
                                                                                                                                      cTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGt                <  2:41761/149‑1 (MQ=255)
                                                                                                                                      cTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGt                >  2:421238/1‑149 (MQ=255)
                                                                                                                                      cTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGt                <  2:869116/149‑1 (MQ=255)
                                                                                                                                       aaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtc               >  4:463230/6‑149 (MQ=255)
                                                                                                                                        gtgAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtct              >  3:305243/1‑149 (MQ=255)
                                                                                                                                        aatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCaga                                                               >  3:373254/5‑100 (MQ=255)
                                                                                                                                        aatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCaga                                                               <  4:373254/96‑1 (MQ=255)
                                                                                                                                               cTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGa       <  3:582530/149‑1 (MQ=255)
                                                                                                                                                tgtgGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGAc      <  3:463230/149‑1 (MQ=255)
                                                                                                                                                    ccACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTa  >  4:117465/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTA  >  NC_000913/1572780‑1573076

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