Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,173,296 G→C 10.2% intergenic (‑289/+43) fadD‑I ← / ← fadD‑II long‑chain‑fatty acid‑‑CoA ligase/long‑chain‑fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,173,2960GC10.2% 289.2 / 22.8 99intergenic (‑289/+43)fadD‑I/fadD‑IIlong‑chain‑fatty acid‑‑CoA ligase/long‑chain‑fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (47/41);  new base C (6/4);  total (54/45)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

Only 100 of 106 total aligned reads displayed.

GGCCAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTGG  >  NC_002947/5173147‑5173445
                                                                                                                                                     |                                                                                                                                                     
ggCCAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTggg                                                                                                                                                       <  2:1521539/150‑1 (MQ=255)
   cAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATAAAGGTCATTTTGAGAGTAATTATCCAAAATGTCGCAGATAAGGATGTATAGAGCCTGTGACAGCGGTTTAAGGCATTCCCGGGTAAACCGGTTCCTACATCGGGCaa                                                                                                                                                    >  2:68893/1‑149 (MQ=255)
    aaaaCAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAgt                                                                                                                                                   >  2:773766/1‑150 (MQ=255)
         aGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAgg                                                                                                                                              >  1:461242/1‑150 (MQ=255)
          ggCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGa                                                                                                                                             <  1:1337401/150‑1 (MQ=255)
           gCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGATCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGTGCAGTGTa                                                                                                                                                >  2:401606/1‑146 (MQ=255)
           gCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTa                                                                                                                                                <  1:401606/146‑1 (MQ=255)
           gCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAg                                                                                                                                            >  2:1957354/1‑150 (MQ=39)
              ttATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCgg                                                                                                                                         <  1:773766/150‑1 (MQ=255)
              ttATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCgg                                                                                                                                         <  2:348364/150‑1 (MQ=255)
              ttATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCgg                                                                                                                                         <  2:2148510/150‑1 (MQ=255)
              ttATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCgg                                                                                                                                         <  2:2081149/150‑1 (MQ=38)
               tATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCggg                                                                                                                                        <  1:1957354/150‑1 (MQ=38)
                 tGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGtt                                                                                                                                      >  1:1155294/1‑150 (MQ=255)
                    tCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTACAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTAc                                                                                                                                   >  2:1190068/1‑150 (MQ=255)
                      acaAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGc                                                                                                                                           >  1:2314422/1‑140 (MQ=255)
                      acaAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGc                                                                                                                                           <  2:2314422/140‑1 (MQ=255)
                       caAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCcg                                                                                                                                <  2:153930/150‑1 (MQ=255)
                               aCTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTggg                                                                                                                                                       >  2:1218974/1‑119 (MQ=255)
                               aCTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTggg                                                                                                                                                       <  1:1218974/119‑1 (MQ=255)
                               aCTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCACGaagaag                                                                                                                        <  2:629876/150‑1 (MQ=255)
                                   aaTGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCAt                                                                                                                    >  2:1714486/1‑150 (MQ=255)
                                    aTGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCgg                                                                                                                                         <  1:1904334/128‑1 (MQ=38)
                                       aaCAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCg                                                                                                                >  1:269382/1‑150 (MQ=255)
                                            tACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGGCATCCCGGTATc                                                                                                           >  1:1156974/1‑150 (MQ=255)
                                                ggTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGc                                                                                                       <  2:269382/150‑1 (MQ=255)
                                                  tCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGGGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCta                                                                                                     <  1:1342022/150‑1 (MQ=255)
                                                           gagTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCa                                                                                            >  2:2032889/1‑150 (MQ=255)
                                                                       aaaaTGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCa                                                                                <  1:1714486/150‑1 (MQ=255)
                                                                                         cTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCgcgc                                                                                             >  1:915037/1‑119 (MQ=255)
                                                                                          tGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCACAAAATCTTGCCCACCGGCGt                                                             >  1:1068933/1‑150 (MQ=255)
                                                                                          tGTACAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGt                                                             <  1:1190068/150‑1 (MQ=255)
                                                                                            tataGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCcg                                                                                                                                >  1:2299421/1‑81 (MQ=255)
                                                                                            tataGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCcg                                                                                                                                <  2:2299421/81‑1 (MQ=255)
                                                                                             ataGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCAt                                                          <  1:940314/150‑1 (MQ=255)
                                                                                              taGAGCCTGTTAGAGCGGTTTAAGGCATTCTCGGGTAAACGCGCTCCTACCTTGGTCAGTGTAGTAGCGGGTTTACCCGCAAAAAAGCCATCACGCTATCACGCTATGTCGCGCGACTGCTGCCGCAAGAACTTGCCCAGTGGCGTCATg                                                         >  2:1205726/1‑150 (MQ=255)
                                                                                              taGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCggg                                                                                                                                        <  2:2245824/71‑1 (MQ=255)
                                                                                              taGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCggg                                                                                                                                        >  1:2245824/1‑71 (MQ=255)
                                                                                                gagCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATggg                                                       >  1:541502/1‑150 (MQ=255)
                                                                                                gagCCTGTTAGAGCGGTTTAAGGCATTCGCGGATAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATggg                                                       >  1:1489983/1‑150 (MQ=255)
                                                                                                   ccTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCgaagaa                                                                                                                         <  1:1254898/81‑1 (MQ=255)
                                                                                                   ccTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCgaagaa                                                                                                                         >  2:1254898/1‑81 (MQ=255)
                                                                                                   ccTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCaa                                                    >  1:1405084/1‑150 (MQ=255)
                                                                                                   ccTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCaa                                                    >  1:1395385/1‑150 (MQ=255)
                                                                                                    cTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAAc                                                   <  2:1155294/150‑1 (MQ=255)
                                                                                                      gTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGt                                                                                                                                       >  1:526803/1‑64 (MQ=255)
                                                                                                      gTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGt                                                                                                                                       <  2:526803/64‑1 (MQ=255)
                                                                                                      gTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATCGGGTCGTCTATGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGa                                                 <  1:2273039/150‑1 (MQ=39)
                                                                                                      gTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCACTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGa                                                 >  1:1033106/1‑150 (MQ=255)
                                                                                                       ttAGAGTGGTTTAAGGCATTCGCGGGTAACCCCGTTCCTACATTGGGTAGTGCGGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGAt                                                <  1:68893/150‑1 (MQ=255)
                                                                                                       ttAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATc                                                                           >  1:347557/1‑123 (MQ=255)
                                                                                                        tAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGAtt                                               <  1:2032889/150‑1 (MQ=255)
                                                                                                           agCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTcgc                                            >  2:393448/1‑150 (MQ=255)
                                                                                                             cGGTTTAAGGCATTCGCGGATAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCgcgc                                          <  2:1489983/150‑1 (MQ=255)
                                                                                                              ggattAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTcgcgaa                                         >  1:2334124/4‑148 (MQ=255)
                                                                                                              ggTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCgca                                         <  2:2124864/150‑1 (MQ=255)
                                                                                                               gTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCgcag                                        >  2:1922812/1‑150 (MQ=255)
                                                                                                                      gCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGc                                                               <  1:1342514/120‑1 (MQ=255)
                                                                                                                      gCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGc                                                               >  2:1342514/1‑120 (MQ=255)
                                                                                                                        aTTCGCGGGTAAACCCGCTTCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGAtgtg                               >  1:671174/1‑150 (MQ=255)
                                                                                                                        aTTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGAtgtg                               >  2:588954/1‑150 (MQ=255)
                                                                                                                           cgcgGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                            <  2:461242/150‑1 (MQ=255)
                                                                                                                           cgcgGAAAAACCCACTCCTACATTGGGCAGTGGAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAACCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                            <  1:1205726/150‑1 (MQ=255)
                                                                                                                            gcgGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCa                                                           >  2:58744/1‑118 (MQ=255)
                                                                                                                            gcgGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCa                                                           <  1:58744/118‑1 (MQ=255)
                                                                                                                              gggTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCgg                                                                                                               >  2:1740081/1‑64 (MQ=255)
                                                                                                                              gggTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCgg                                                                                                               <  1:1740081/64‑1 (MQ=255)
                                                                                                                              gggTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTggg                         >  2:516812/1‑150 (MQ=255)
                                                                                                                                   aaCCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACtt                    <  2:2158674/150‑1 (MQ=255)
                                                                                                                                   aaCCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACtt                    <  2:1405084/150‑1 (MQ=255)
                                                                                                                                    aCCCGCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGc                                                                       <  2:1543389/98‑1 (MQ=25)
                                                                                                                                     cccGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTg                  >  1:2008878/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTg                  >  2:1836033/1‑150 (MQ=38)
                                                                                                                                     cccGCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGAGGGGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTg                  >  2:81951/1‑150 (MQ=38)
                                                                                                                                      ccacTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATc                                                                           <  1:355425/89‑1 (MQ=255)
                                                                                                                                      ccacTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATc                                                                           >  2:355425/4‑92 (MQ=255)
                                                                                                                                      ccGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                            <  2:202630/139‑1 (MQ=255)
                                                                                                                                      ccGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                            >  1:202630/1‑139 (MQ=255)
                                                                                                                                       cGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCCACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTa                >  2:1134313/1‑150 (MQ=255)
                                                                                                                                       cGCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTa                >  1:617216/1‑150 (MQ=39)
                                                                                                                                         cTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGc              <  2:2008878/150‑1 (MQ=255)
                                                                                                                                          tCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATg                                                         <  1:405593/106‑1 (MQ=255)
                                                                                                                                          tCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATg                                                         >  2:405593/1‑106 (MQ=255)
                                                                                                                                           ccTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAGGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCAGCTTTACAGGAGTGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCg            <  2:2322123/150‑1 (MQ=255)
                                                                                                                                           ccTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCg            >  2:861722/1‑150 (MQ=255)
                                                                                                                                           ccTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCg            >  2:2099784/1‑150 (MQ=255)
                                                                                                                                            cTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGc                                                                       <  2:368513/90‑1 (MQ=255)
                                                                                                                                            cTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGc                                                                       >  1:368513/1‑90 (MQ=255)
                                                                                                                                            cTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCgg           >  1:1829178/1‑150 (MQ=255)
                                                                                                                                             tACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGt          >  1:1284875/1‑150 (MQ=255)
                                                                                                                                               cATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                            <  1:2047743/130‑1 (MQ=255)
                                                                                                                                               cATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                            >  2:2047743/1‑130 (MQ=255)
                                                                                                                                               cATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTg                           >  2:1327959/1‑131 (MQ=255)
                                                                                                                                               cATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTg                           <  1:1327959/131‑1 (MQ=255)
                                                                                                                                               cATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGa        >  2:195950/1‑150 (MQ=255)
                                                                                                                                                aTTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGaa       >  2:1837550/1‑150 (MQ=255)
                                                                                                                                                    ggCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTg   <  1:195950/150‑1 (MQ=255)
                                                                                                                                                     gCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTgg  >  2:1336342/1‑150 (MQ=255)
                                                                                                                                                     gCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTgg  >  1:1250453/1‑150 (MQ=255)
                                                                                                                                                     gCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTgg  <  2:1829178/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCCAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTGG  >  NC_002947/5173147‑5173445

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