Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,725,886 G→T 6.1% Q774K (CAG→AAG)  plsB ← glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,725,8860GT 6.1% 210.6 / 2.2 65Q774K (CAG→AAG) plsBglycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (26/35);  new base T (2/2);  total (28/37)
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.68e-01

AGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGTTGATC  >  NC_002947/1725738‑1726034
                                                                                                                                                    |                                                                                                                                                    
aGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGa                                                                                                                                                     >  2:73235/1‑150 (MQ=255)
aGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTCTTCTTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTTCTGa                                                                                                                                                     <  1:101604/150‑1 (MQ=255)
 gCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACGGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGaa                                                                                                                                                    <  1:424741/150‑1 (MQ=255)
     gtcgtcACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCaa                                                                                                                                                <  1:1103966/150‑1 (MQ=255)
     gtcgtcACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCaa                                                                                                                                                <  2:429791/150‑1 (MQ=255)
            cGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTTCAGCAAggcgtgc                                                                                                                                         <  1:656568/150‑6 (MQ=255)
             ggTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtgaa                                                                                                                                        >  2:611953/1‑150 (MQ=255)
              gTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtgaat                                                                                                                                       <  2:773166/150‑1 (MQ=255)
                 cAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtgaatgaa                                                                                                                                    <  1:469842/150‑1 (MQ=255)
                 cAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtgaatgaa                                                                                                                                    <  2:1391572/150‑1 (MQ=255)
                   gggCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGt                                                                                                                                  <  1:792418/150‑1 (MQ=255)
                          cTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGaa                                                                                                                           <  2:652032/150‑1 (MQ=255)
                                gATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGt                                                                                                                     <  1:73235/150‑1 (MQ=255)
                                         gcgcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCgaa                                                                                                            <  2:1397189/150‑1 (MQ=255)
                                           gcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGt                                                                                                                >  1:161763/1‑144 (MQ=255)
                                           gcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGt                                                                                                                <  2:161763/144‑1 (MQ=255)
                                             aGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAAc                                                                                                        >  2:1321256/1‑150 (MQ=255)
                                                   gCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGt                                                                                                  <  2:323983/150‑1 (MQ=255)
                                                   gCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGt                                                                                                  <  1:578577/150‑1 (MQ=255)
                                                           tACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGc                                                                                                                       <  2:1334937/121‑1 (MQ=255)
                                                           tACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGc                                                                                                                       >  1:1334937/1‑121 (MQ=255)
                                                              tCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGGTAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGt                                                                                                                <  1:798872/125‑1 (MQ=255)
                                                              tCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGt                                                                                                                >  2:798872/1‑125 (MQ=255)
                                                              tCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAg                                                                                          <  2:1312135/147‑1 (MQ=255)
                                                              tCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAg                                                                                          >  1:1312135/1‑147 (MQ=255)
                                                              tCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGcc                                                                                       >  2:552189/1‑150 (MQ=255)
                                                              tCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGcc                                                                                       >  1:1100651/1‑150 (MQ=255)
                                                                gCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCAt                                                                                     >  2:1332596/1‑150 (MQ=255)
                                                                gCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCAt                                                                                     >  1:967135/1‑150 (MQ=255)
                                                                gCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCAt                                                                                     >  2:835609/1‑150 (MQ=255)
                                                                    ggCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAg                                                                                 >  1:336840/1‑150 (MQ=255)
                                                                        aCGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATc                                                                             <  1:978816/150‑1 (MQ=255)
                                                                         cGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGc                                                                                                                               <  2:802412/99‑1 (MQ=255)
                                                                         cGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGc                                                                                                                               >  1:802412/1‑99 (MQ=255)
                                                                               cGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGgc                                                                      <  1:611953/150‑1 (MQ=255)
                                                                                    aGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTgg                                                                 <  1:835609/150‑1 (MQ=255)
                                                                                    aGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTgg                                                                 <  2:35829/150‑1 (MQ=255)
                                                                                           cGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAgg                                                                                         >  2:1120203/1‑119 (MQ=255)
                                                                                           cGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAgg                                                                                         <  1:1120203/119‑1 (MQ=255)
                                                                                              gCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTACACAGGAATGCGTGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCAtcatca                                                       <  2:1416631/150‑1 (MQ=255)
                                                                                                   tCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAg                                                                                          >  2:1347753/1‑110 (MQ=255)
                                                                                                   tCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAg                                                                                          <  1:1347753/110‑1 (MQ=255)
                                                                                                   tCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGcac                                                  >  1:776232/1‑150 (MQ=255)
                                                                                                    cGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGGCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGAGCTGGGTCATCATCACGcaca                                                 <  2:490367/150‑1 (MQ=255)
                                                                                                     gTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAg                                                                                          <  2:1142650/108‑1 (MQ=255)
                                                                                                     gTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAg                                                                                          >  1:1142650/1‑108 (MQ=255)
                                                                                                             tCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCa                                                                                  >  1:925476/1‑108 (MQ=255)
                                                                                                             tCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCa                                                                                  <  2:925476/108‑1 (MQ=255)
                                                                                                                    gccttgcGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCAc                                                      <  2:1282868/129‑1 (MQ=255)
                                                                                                                    gccttgcGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCAc                                                      >  1:1282868/1‑129 (MQ=255)
                                                                                                                      cttgcGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCgg                               >  2:1177942/1‑150 (MQ=255)
                                                                                                                       ttgcGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGTTCTTCCAGCGCTTCGgc                              >  1:1092916/1‑150 (MQ=255)
                                                                                                                         gcgcGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGgcgc                            >  1:868541/1‑150 (MQ=255)
                                                                                                                            cgTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGCCTGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCa                         <  1:1144242/150‑1 (MQ=255)
                                                                                                                               cGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGc                      <  2:336840/150‑1 (MQ=255)
                                                                                                                                   gTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCTATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGtt                  <  1:1332596/150‑1 (MQ=255)
                                                                                                                                        gcaCGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGcc             >  2:795668/1‑150 (MQ=255)
                                                                                                                                           gaccTTGCTTCAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGc              >  1:698004/3‑146 (MQ=255)
                                                                                                                                           gaccTTGCTTCAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGc              <  2:698004/144‑1 (MQ=255)
                                                                                                                                           cGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGc          <  1:709851/150‑1 (MQ=255)
                                                                                                                                            accTTGCTTCAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCt         >  1:1256662/2‑150 (MQ=255)
                                                                                                                                              cTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGt       >  1:451119/1‑150 (MQ=255)
                                                                                                                                                 gCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGTTGa    <  1:1321256/150‑1 (MQ=255)
                                                                                                                                                   tGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGTTGATc  >  1:654223/1‑150 (MQ=255)
                                                                                                                                                   tGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGTTGATc  <  2:1100651/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGTTGATC  >  NC_002947/1725738‑1726034

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