Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 391,002 A→C 8.8% S425R (AGC→CGC)  soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947391,0020AC 8.8% 187.3 / 10.3 68S425R (AGC→CGC) soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (33/29);  new base C (3/3);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

GCATACTCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATGTCATTGCCGAT  >  NC_002947/390853‑391151
                                                                                                                                                      |                                                                                                                                                     
gCATACTCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACa                                                                                                                                                       <  2:1369238/150‑1 (MQ=255)
 cATACTCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAg                                                                                                                                                      >  1:489554/1‑150 (MQ=255)
      tCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGa                                                                                                                                                 >  2:735214/1‑150 (MQ=255)
        tcttCCAGCGCCGTGCTCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATc                                                                                                                                               <  2:684196/147‑1 (MQ=255)
        aCCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGcccatga                                                                                                                                               >  2:1218065/1‑146 (MQ=255)
         ccTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCTGTAACCGCCCGagga                                                                                                                                              <  2:796973/150‑4 (MQ=255)
            ccAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGGGGCTGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCgaga                                                                                                                                           >  1:684196/1‑147 (MQ=255)
            ccAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCCACTTATCGCTGGCGTGCGATGTGCAGGCGACCAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGc                                                                                                                                           <  1:1028903/150‑1 (MQ=255)
               gcgcCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCAcc                                                                                                                                        >  1:702448/1‑150 (MQ=255)
                     tGATCGAGGCCAAAGGCAGCAAGCACGTCACAGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCtt                                                                                                                                  >  1:1177059/1‑150 (MQ=255)
                      gATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTc                                                                                                                                 >  2:918722/1‑150 (MQ=255)
                        tCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGc                                                                                                                               <  1:262187/150‑1 (MQ=255)
                         cGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCa                                                                                                                              >  1:635185/1‑150 (MQ=255)
                          gAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTg                                                                                                                                            >  2:508445/1‑135 (MQ=255)
                          gAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTg                                                                                                                                            <  1:508445/135‑1 (MQ=255)
                                aaaGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGc                                                                                                                       <  1:1254890/150‑1 (MQ=255)
                                 aaGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGCTGGATGGCCATGCGCACAAAGTCACAAGCCAAGGCGAAGTCATTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTAAAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCg                                                                                                                      >  2:984823/1‑150 (MQ=255)
                                        caagcaCGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGccc                                                                                                               <  2:117826/150‑1 (MQ=255)
                                           gcaCGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGt                                                                                                            <  2:489554/150‑1 (MQ=255)
                                                tCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGACCGcccaggt                                                                                                           >  2:1091645/1‑142 (MQ=255)
                                                ccACCGGCGCCCGTGTGGCGGCGACCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGACCGcccaggt                                                                                                           <  1:1091645/145‑5 (MQ=255)
                                                  aCCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGt                                                                                                     <  1:1323408/150‑1 (MQ=255)
                                                                     gcgATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTc                                                                                   >  1:496599/1‑150 (MQ=255)
                                                                           gatGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGc                                                                                                       <  1:624247/124‑1 (MQ=255)
                                                                           gatGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGc                                                                                                       >  2:624247/1‑124 (MQ=255)
                                                                           gatGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCg                                                                             <  2:1177059/150‑1 (MQ=255)
                                                                            atGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCgg                                                                            >  2:1036360/1‑150 (MQ=255)
                                                                                  aGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCTGTAACCGCCCGag                                                                                                                                                >  1:1371800/1‑75 (MQ=255)
                                                                                  aGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCTGTAACCGCCCGag                                                                                                                                                <  2:1371800/76‑2 (MQ=255)
                                                                                   ggCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATgc                                                                     <  1:735214/150‑1 (MQ=255)
                                                                                     cgcACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCGCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATgcgc                                                                   <  2:756054/150‑1 (MQ=255)
                                                                                      gcACAAAGTCACCAGCCCAGGCGAAATCCTTGAGTGCGACCTGATCGACTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGCTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGACGGGCGATGcacc                                                                  >  2:190862/1‑150 (MQ=255)
                                                                                          aaaGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAg                                                              >  2:358625/1‑150 (MQ=255)
                                                                                                aCCAGCCCAGGCGAAGTCCTTGAGTGTGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGt                                                        >  1:226453/1‑150 (MQ=255)
                                                                                                aCCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTgacga                                                                                                >  2:99135/1‑110 (MQ=255)
                                                                                                aCCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTgacga                                                                                                <  1:99135/110‑1 (MQ=255)
                                                                                                  cAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGa                                                      <  2:88670/150‑1 (MQ=255)
                                                                                                    gCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCTGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGt                                                    <  1:1159097/150‑1 (MQ=255)
                                                                                                          gCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGa                                                      >  2:471402/1‑142 (MQ=255)
                                                                                                          gCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGa                                                      <  1:471402/142‑1 (MQ=255)
                                                                                                          gCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCg                                              >  2:327717/1‑150 (MQ=255)
                                                                                                            gAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTGCAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCggt                                            >  1:1358991/1‑150 (MQ=255)
                                                                                                             ccgatcTTGAGTGCGACATGATCGACACCATCGACGGTTACATCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCggtg                                           <  2:866736/145‑1 (MQ=255)
                                                                                                                  cTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACAt                                                                                             >  2:20080/1‑95 (MQ=255)
                                                                                                                  cTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACAt                                                                                             <  1:20080/95‑1 (MQ=255)
                                                                                                                  cTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATc                                      >  2:1343598/1‑150 (MQ=255)
                                                                                                                   ttGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCggtgg                                          >  2:1399609/1‑145 (MQ=255)
                                                                                                                   ttGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCggtgg                                          <  1:1399609/145‑1 (MQ=255)
                                                                                                                   ttGAGTGCGACCTGATCGCCTCCTAGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGt                                                        >  1:930503/1‑131 (MQ=255)
                                                                                                                   ttGAGTGCGACCTGATCGCCTCCTAGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGt                                                        <  2:930503/131‑1 (MQ=255)
                                                                                                                     gAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCggtggagatcgg                                   >  1:866736/1‑143 (MQ=255)
                                                                                                                     gAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAAc                                   <  1:354637/150‑1 (MQ=255)
                                                                                                                        tGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGc                                <  1:984823/150‑1 (MQ=255)
                                                                                                                        tGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGc                                >  1:487130/1‑150 (MQ=255)
                                                                                                                           gACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCgtg                             >  2:1233111/1‑150 (MQ=255)
                                                                                                                              cTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCggt                                            <  2:833366/132‑1 (MQ=255)
                                                                                                                              cTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCggt                                            >  1:833366/1‑132 (MQ=255)
                                                                                                                                   tgcacgatcGGGCTGTAACCGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCtt                     >  1:1229524/8‑150 (MQ=255)
                                                                                                                                      ctcctcGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTggg                  >  2:286094/1‑150 (MQ=255)
                                                                                                                                        gatcGGGCTGTAACCGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCggtg                                           <  1:630357/121‑1 (MQ=255)
                                                                                                                                        gatcGGGCTGTAACCGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCggtg                                           >  2:630357/3‑123 (MQ=255)
                                                                                                                                          tcGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCggtg                                           >  2:359668/1‑121 (MQ=255)
                                                                                                                                          tcGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCggtg                                           <  1:359668/121‑1 (MQ=255)
                                                                                                                                            gggCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGc                                                  <  2:658676/112‑1 (MQ=255)
                                                                                                                                            gggCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGc                                                  >  1:658676/1‑112 (MQ=255)
                                                                                                                                                 gTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCg                                                                                  <  1:985326/75‑1 (MQ=255)
                                                                                                                                                 gTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCg                                                                                  >  2:985326/1‑75 (MQ=255)
                                                                                                                                                  ttACAGCCCGATCGTGGTCCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATGTCATTGc      <  2:516398/150‑1 (MQ=255)
                                                                                                                                                  ttACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATGTCATTGc      <  1:1395609/150‑1 (MQ=255)
                                                                                                                                                      aGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATgcgc                                                                   <  2:345087/85‑1 (MQ=255)
                                                                                                                                                      aGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATgcgc                                                                   >  1:345087/1‑85 (MQ=255)
                                                                                                                                                      aGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATGTCATTGCCGAt  <  1:1343598/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCATACTCACCTCCAGCGCCGTGATCGAGGCCAAAGGCAGCAAGCACGTCACCGGCGCCCGTGTGGCGGCGATCGATGTGCAGGCGCACAAAGTCACCAGCCCAGGCGAAGTCCTTGAGTGCGACCTGATCGCCTCCTCGGGCGGTTACAGCCCGATCGTGCACCTGGCTTCGCACCTGGGCGGTCGCCCGGTATGGCGTGACGACATCCTTGGCTTCGTGCCGGGCGATGCGCCGCAGAAGCGTGAGTGCGTCGGTGGTATCAACGGCGTGTATGCCTTGGGCGATGTCATTGCCGAT  >  NC_002947/390853‑391151

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