Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,316,831 C→T 6.1% E130K (GAG→AAG)  phnM ← ribophosphonate triphosphate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,316,8310CT 6.1% 212.6 / 3.8 66E130K (GAG→AAG) phnMribophosphonate triphosphate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (36/26);  new base T (2/2);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

ATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGC  >  NC_000913/4316685‑4316975
                                                                                                                                                  |                                                                                                                                                
aTACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCgg                                                                                                                                                >  4:457237/1‑149 (MQ=255)
aTACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGAGGTCCATCAGCGACACCAGCGGCACCGGCTCGCCCTGCACCAGTTTTTCAAACTGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCgt                                                                                                                                                >  2:224146/1‑148 (MQ=255)
   cTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGgcgc                                                                                                                                             <  1:597412/149‑1 (MQ=255)
    ttCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGcgcg                                                                                                                                            >  1:162207/1‑149 (MQ=255)
      ctcGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGt                                                                                                                                          >  4:70698/1‑149 (MQ=255)
             ttGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGAcgcc                                                                                                                                   <  3:136803/149‑1 (MQ=255)
                        gcgcTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCg                                                                                                                        >  3:314151/1‑149 (MQ=255)
                          gcTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGtc                                                                                                                      >  1:462879/1‑149 (MQ=255)
                               cccGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACCGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTc                                                                                                                 <  1:610809/149‑1 (MQ=255)
                               cccGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTc                                                                                                                 <  2:432889/149‑1 (MQ=255)
                                 cGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGa                                                                                                               <  4:219517/149‑1 (MQ=255)
                                        tGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCACACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGtt                                                                                                        <  2:462879/149‑1 (MQ=255)
                                          gTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGa                                                                                                      <  1:587342/149‑1 (MQ=255)
                                           tCCATCAGCGAGACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGat                                                                                                     >  3:536335/1‑149 (MQ=255)
                                           tCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGat                                                                                                     <  1:231640/149‑1 (MQ=255)
                                             cATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTTGGCGCGGTTGACGCCGCGTTTCTGCGTctct                                                                                                                   <  1:305372/133‑1 (MQ=255)
                                             cATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTTGGCGCGGTTGACGCCGCGTTTCTGCGTctct                                                                                                                   >  2:305372/1‑133 (MQ=255)
                                             cATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTTGGCGCGGTTGACGCCGCGTTTCTGCGTctct                                                                                                                   >  4:525643/1‑133 (MQ=255)
                                             cATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTTGGCGCGGTTGACGCCGCGTTTCTGCGTctct                                                                                                                   <  3:525643/133‑1 (MQ=255)
                                             cATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGatca                                                                                                   >  1:425458/1‑149 (MQ=255)
                                             cATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGatca                                                                                                   >  3:155386/1‑149 (MQ=255)
                                              aTCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGatcat                                                                                                  >  3:143567/1‑149 (MQ=255)
                                               tCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGatcatc                                                                                                 >  3:521096/1‑149 (MQ=255)
                                                       acCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAg                                                                                         >  4:108442/1‑149 (MQ=255)
                                                                aCCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAg                                                                                >  3:523494/1‑149 (MQ=255)
                                                                      tCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGa                                                                                                      >  3:354273/1‑121 (MQ=255)
                                                                      tCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGa                                                                                                      <  4:354273/121‑1 (MQ=255)
                                                                      tCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATc                                                                          <  4:274594/149‑1 (MQ=255)
                                                                        gcgcTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATcgc                                                                        <  3:108442/149‑1 (MQ=255)
                                                                         cgcTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATcgcc                                                                       <  3:70698/149‑1 (MQ=255)
                                                                          gcTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGa                                                                                                      >  1:555586/1‑117 (MQ=255)
                                                                          gcTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGa                                                                                                      <  2:555586/117‑1 (MQ=255)
                                                                                  aGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTcgcg                                                              >  2:17383/1‑149 (MQ=255)
                                                                                      tttCAAACAGCGGCGGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACg                                                          >  3:43514/1‑149 (MQ=255)
                                                                                          aaaCAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGc                                                      >  3:313987/1‑149 (MQ=255)
                                                                                            aCAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCa                                                    >  2:64867/1‑149 (MQ=255)
                                                                                            aCAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCa                                                    >  2:552534/1‑149 (MQ=255)
                                                                                              aGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCAtct                                                                                                <  3:197344/103‑1 (MQ=255)
                                                                                              aGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCAtct                                                                                                >  4:197344/1‑103 (MQ=255)
                                                                                              aGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAAt                                                  >  3:201540/1‑149 (MQ=255)
                                                                                              aGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAAt                                                  >  1:37141/1‑149 (MQ=255)
                                                                                                   cAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCa                                             <  2:425458/149‑1 (MQ=255)
                                                                                                           tgtgATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATcc                                     >  2:502929/1‑149 (MQ=255)
                                                                                                             tgATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAg                                   >  3:213819/1‑149 (MQ=255)
                                                                                                               aTGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTa                                 >  4:315358/1‑149 (MQ=255)
                                                                                                               aTGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTa                                 >  4:404766/1‑149 (MQ=255)
                                                                                                                    gCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACggtg                            <  2:77194/149‑1 (MQ=255)
                                                                                                                      aGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACggtggt                          <  2:162207/149‑1 (MQ=255)
                                                                                                                        ttCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGa                        >  3:417672/1‑149 (MQ=255)
                                                                                                                         tCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGAt                       >  2:250141/1‑149 (MQ=255)
                                                                                                                           gcagcgcagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGc                     >  3:32199/1‑149 (MQ=255)
                                                                                                                           acagcgcagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGc                     <  3:53506/148‑1 (MQ=255)
                                                                                                                                gcagATGCAGACGGTGCTCGGCGCGTTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTc                <  3:404766/149‑1 (MQ=255)
                                                                                                                                gcagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTc                <  4:155386/149‑1 (MQ=255)
                                                                                                                                gcagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTc                >  1:351998/1‑149 (MQ=255)
                                                                                                                                 cagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCaa                                                   >  4:235777/1‑113 (MQ=255)
                                                                                                                                 cagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCaa                                                   <  3:235777/113‑1 (MQ=255)
                                                                                                                                 cagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTAc                                <  4:165593/132‑1 (MQ=255)
                                                                                                                                 cagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTAc                                >  1:345445/1‑132 (MQ=255)
                                                                                                                                 cagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTAc                                >  3:165593/1‑132 (MQ=255)
                                                                                                                                 cagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTAc                                <  2:345445/132‑1 (MQ=255)
                                                                                                                                 cagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCg               >  3:173225/1‑149 (MQ=255)
                                                                                                                                   gATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCGcc             <  3:412082/149‑1 (MQ=255)
                                                                                                                                   gATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCGcc             <  1:39556/149‑1 (MQ=255)
                                                                                                                                        aGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATCCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCAt        >  1:321746/1‑149 (MQ=255)
                                                                                                                                              tGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAgcgc  <  4:232005/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGC  >  NC_000913/4316685‑4316975

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